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PREFACE 


GENERAL 

The Fractional Watt Vuilleumier Cryogenic Refrigerator Program was ad- 
ministered by the Goddard Space Flight Center, Greenbelt, Maryland, under 
Contract No. NAS 5-21715* The NASA Technical Monitor for this contract was 
Mr. Max Gasser. 

The Task 1 preliminary design activities were completed by AiResearch 
Manufacturing Company, a division of The Garrett Corporation, in Torrance, 
California. The Program Manager was V. L. Potter and the Program Engineer 
was W. S. Miller. Major contributors to this Task I effort were V. L. Potter, 
mechanical design; M. Morimoto, stress analysis; M. Soliman, heat pipes; 

J. Denk, motor; and W. S. Miller, thermodynamics and heat transfer. In 
addition, C. W. Browning served as technical consultant. 

OBJECTIVE 

An engineering model long life VM refrigerator has been fabricated and 
delivered to NASA/GSFC under contract NAS 5~ 21 096 . The emphasis in the design 
of the engineering model was the attainment of long life. A natural extension 
to that work is application of the technology developed to flight type hard- 
ware. Thus, the objective of this program is to design and build a flight 
prototype VM refrigerator. 

SCOPE OF WORK 

The VM cryogenic refrigerator program consists of five major tasks: 

(1) Task 1, Preliminary Design; (2) Task I I , Anal yt i ca 1 and Test Programs; 

(3) Task III, Model Specification; (4) Task IV, Model Fabrication and Test, 
and (5) Task V, Program Management and Periodic Documentation. 

This technical teport describes the Task I (Preliminary Design) phase of 
the VM cryogenic refrigerator design, development, and fabrication program. 

The general requirement of the refrigerator is to produce 0.25 watts of cool- 
ing at 65° K after operating 2 to 5 years in a space environment. These 
Task I accomplishments provide the basis for the further design, analysis, 
development testing, and evaluation to be conducted during the next phase 
(Task 1 1 ) ' o f this program. 
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SECTION I 


INTRODUCTION 


The progressing development of meteorological and communication spacecraft 
requires increasingly sophisticated sensing devices. A coexisting requirement 
is the development of cryogenic refrigeration; the sensitivity and performance 
of certain advanced sensing devices is enhanced when cooled to cryogenic tem- 
peratures, while others must be cooled to cryogenic temperatures in order to 
f unct i on. 

Most NASA space applications for cryogenic refrigeration require opera- 
tional lifetimes of I to 5 years. State-of-the-art production VM refrigerators, 
however, have life expectancies far less than the minimum acceptable for space 
applications. For this reason, NASA undertook the development of a long life 
cryogenic refrigeration system compatible with spacecraft requirements. The 
heat driven Vuilleumier (VM) cycle was the refrigeration system selected for 
development because of the following characteristics that are potentially 
advantageous for space applications: (l) long operating life due to low work- 

ing stresses placed on bearings and internal seals, (2) use of thermal energy 
as the source of energy to drive the system, (3) small size and weight, (4) 
moderate efficiencies for low cooling capacity systems, and (5) low noise and 
vibration levels. 

^Government funded efforts in this area have shown significant progress in 
demonstrating the characteristics listed above, including long life. An 
engineering model long life VM refrigerator has been fabricated and delivered 
to NASA/GSFC under Contract NAS 5-21096. To demonstrate the feasibility of 
long life VM refrigerators for NASA space applications, Contract NAS 5-21715 
was entered into in March 1972 between AiResearch and the NASA Goddard Space 
Flight Center. The contract covers design, development and fabrication of a 
flight prototype fractional watt VM refrigerator for integration into a radi- 
ometer device. The primary design requirements of the unit being developed 
are long operating life, with a 2 year operating life specified as a minimum 
requirement, and the use of design techniques and assembly methods compatible with 
unmanned satellite environmental conditions. The refrigerator shall be capable 
of cooling a 1/4 watt thermal load to 65°K. 

The program under Contract No. NAS 5-21715 consists of the following 
major tasks: 

Task I Preliminary Design 

Task II Analytical and Test Program 

Task III Model Specification 

Task IV Model Fabrication and Test 

Task V Program Management and Periodic Documentation 
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The technical specification for this program is given by NASA/GSFC Statement 
of Work and Specification APD-763-28, dated September 15, 1971. 

This report describes the Task I effort, which consisted of definition of 
the Preliminary Design of the Flight Prototype Fractional Watt Vuilleumier 
C ryogeni c Ref r i gerator . 

TASK I PRELIMINARY DESIGN 

The principal goal of the Task I effort was to establish and recommend a 
preliminary design VM refrigerator which would meet the operational require- 
ments for the intended application. The general design requirements are that 
the refrigerator cool a 1/4 watt load at 65°K. A summary of the preliminary 
design VM refrigerator characteristics is presented in Section 2 of this report. 
The detailed requirements and the technical discussion of the procedures uti- 
lized to arrive at the preliminary design configuration are presented in Section 
3. 


The Task I effort was delineated into several major subtasks, which are 
described briefly in the following paragraphs. The experience gained in the 
design and fabrication of the engineering model, GSFC 5 watt VM refrigerator, 
was utilized to the greatest possible extent during the present program. 

Subtask 1-1 --Prel imi nary Engine Sizing and Thermal Analysis 

The sizing of the refrigerator and its thermal analysis paralleled the 
mechanical design, which was a concurrent effort. The thermal analysis forms 
a major part of the design iteration loop, and thus was integrated into the 
overall preliminary design effort. 

The initial design parameters used as the starting point for the thermal 
analysis were those derived for the conceptual design evolved at the beginning 
of the program. After establishing revised design parameters, including inputs 
from the mechanical design analysis and other associated tasks, the thermo- 
dynamic cycle was reevaluated. By repeating this process while holding the 
basic operating parameters constant — cycle pressure, hot-end temperature, ambient 
temperature, refrigeration temperature, and machine speed — the optimum physical 
characteristics of a refrigerator (based on the initial design parameters) were 
defined. Then, changing those basic parameters which are subject to variation 
for purposes of optimization, the characteristics based on a new set of param- 
eters were defined. By comparing these successive iterations, the optimum 
preliminary design was selected. 

Definition of the operating cycle described above allows determination of 
the required thermal characteristics of the various critical components of the 
machine. As these devices were analytically defined, modifications to the 
operating characteristics were dictated. These changes were fed back into the 
cycle analysis and the resulting changes in thermal requirements determined. 

The individual components were then modified, and the iterative procedure 
repeated. The entire VM refrigerator was thus iterated to its preliminary 
design, or present, configuration. 
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Subtask I-2--Prel iminary Mechanical Design and Stress Analysis 


The mechanical design considerations associated with VM refrigerator 
design are many and varied. Special attention was given to areas involving 
the drive mechanism, bearings, seals, refrigerant leakage, and contamination. 
Each of these were considered independently, as well as their interaction on 
one another and their impact on the thermal design. 

The basic mechanical design approach developed for the previous engineer- 
ing model machine was utilized for the fractional watt VM refrigerator. This 
approach makes use of solid compact bearings and non-contact seals to achieve 
long-life capability. Based on test results on the GSFC 5 watt VM, there is 
every reason to believe that the fractional watt VM refrigerator will meet all 
design requirements. 

Two areas that required careful consideration in the mechanical design of 
the fractional watt machine are: (l) scaling down the drive mechanism and 

bearing approach used in the 5-watt machine to reduce weight and size and (2) 
the devising of a method of closure of the working volume of the engine that 
allows take-down of the engine for inspection during the run-in period and 
final close-out of the system by welding to prevent external leakage upon final 
assembly. 

The preliminary mechanical design and stress analysis task also covered 
the structural design of the pressure containment members of the refrigerator 
and the dynamics of the system including the vibration loading of the key ele- 
ments of the design. These analyses will be refined in Task 2. 

Subtask I-3--Heat Pipe Interface 


A preliminary heat pipe design was evolved to evaluate the removable heat 
pipe interface assembly. The heat pipe assembly design was only carried to the 
depth necessary to define the interface requirements. This, in turn, required 
definition of the number of heat pipes, their size, the working fluid, the wick 
configuration, and the materials of construction. 

The heat pipe assembly design makes use of four separate heat pipes (two 
on each side of the refrigerator crankcase). This approach allows direct use 
of the heat pipe interface clamp design for the water-cooled interface. This 
offers a considerable advantage in that complete definition of the interface 
clamp requires extensive stress analysis to ensure good thermal contact under 
all operating conditions. With direct use of the identical clamp, duplicate 
effort was avoided. 

It should be noted that the originally specified I80°F sump temperature 
was reduced to I40°F by program Technical Direction No. I, thus allowing use 
of ammonia as the heat pipe working fluid. 
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Sublask l-4--Motor Design and Specification 


A preliminary specification for the drive motor was prepared and submitted 
to the AiResearch motor group for review. It was determined that lamination 
dies were available in the size range and required pole pair- number at AiResearch, 
the products of previous production fabrication orders. Based on this avail- 
ability, further evaluation of other potential sources was eliminated. 

The EMI specifications submitted by Honeywell for review were also subjected 
to analysis by the motor EMI Analysts. It is the consensus that the 
AC motor depicted in the VM layout drawing is a device that is free of radiated 
and conducted EMI. This is, however, an opinion since no one with the exception 
of Honeywell has tested at the levels required by their specification. 

Subtask I-5--Layout Drawings 

Layout drawings were developed as required in support of the above four 
Task I subtasks. A full system layout drawing has been constructed depicting 
the preliminary design. This drawing, together with the supporting analyses, 
form the basis for the Task I design review as well as the Task I Summary 
Report . 

Subtask I-6--Summary Report 

This Task I Summary Report describes the preliminary system design and 
summarizes the analyses conducted in Task I. 
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SECTION 2 


PRELIMINARY DESIGN SUMMARY 


INTRODUCTION 

The GSFC fractional watt Vuilleumier cryogenic refrigerator is a gas-cycle, 
reciprocating machine, intended to operate continuously for five years at a 
design speed of 400 rpm. Continuous operation for such a long time period pre- 
sents technical problems, namely long-life bearing and seal design. This 
technology has been developed, and is applied to a flight prototype machine 
for the first time in this design. As an aid in understanding the analysis 
described in subsequent sections of this Task I report, a brief physical and 
functional description and summary of the VM refrigerator is given in this sec- 
tion. The preliminary design layout drawing of this refrigerator is presented 
in Figure 2-1 . 

PHYSICAL DESCRIPTION 

Although the primary energy input to the refrigerator is thermal energy, 
moving components (two displacers) are required so that the VM cycle will 
function properly. The two displacers (hot and cold) which make the refrigera- 
tor a reciprocating machine are driven by connecting rods attached to a crank- 
shaft (Figure 2-1). Crankshaft throws are arranged 90-de^rees apart so that 
the hot displacer is leading the cold displacer when the crankshaft is driven 
by an electric motor located at one end of the shaft. The displacers travel 
inside cylinders surrounded by packed-bed regenerators which in turn, are 
enclosed in pressure shells joined to the crankshaft housing, making the entire 
assembly pressure tight. 

The displacers are arranged in a horizontally opposed configuration, with 
the cold displacer positioned at an angle of 180 degrees from the hot displacer, 
thus simplifying the crankcase and sump heat exchanger mechanical design and 
minimizing the interaction between the hottest and coldest parts of the system. 

As indicated earlier, the primary energy for driving the refrigerator is 
introduced at the hot end in the form of heat. The interface at the hot end 
of the machine is shown in Figure 2-1. For this design, an electric heater is 
used to supply the heat; future spacecraft systems may use radioisotope heat 
sources or a solar collector to conserve solar cell produced energy. 

Proper operation of the VM cycle also depends on rejecting heat in the 
crankcase-sump region. The GSFC VM refrigerator rejects heat to water cooling 
coils that interface with the crankcase and sump heat exchanger. These water 
cooling coils replace the ammonia heat pipes and a radiator which would be 
used in a flight-type system. 

The refrigeration heat load is absorbed at the cold end via the cold-end 
heat exchanger. For test purposes, the refrigeration load is generated by a 
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Figure 2-1. GSFC Flight Prototype 
Fractional Watt 
Vuilleumier Cryogenic 
Refrigerator Prelimi- 
nary Design 
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small resistance-type heater bonded to the exterior surface of the cold-end 
heat exchanger. In a spacecraft system, the refrigerator load will be mounted 
directly to the cold end. 

SYSTEM OPERATION 

Externally, the VM refrigerator is viewed as a machine which absorbs heat 
at both a high and low temperature and then rejects this heat at some inter- 
mediate temperature (heat sink temperature). Heat input at the hot end of the 
machine provides the energy required to produce cooling at the cold end, while 
heat is being rejected from the sump or crankcase region of the system. The 
internal processes leading to the production of refrigeration are complex and 
difficult to describe in their real, rather than idealized form. Only a sim- 
plified description of system operation is provided here. 

The total internal volume of the VM refrigerator is constant and consists 
of two parts: an active volume and a void volume. The active volume divides 

into three regions at different temperatures and of varying size due to recip- 
rocating displacer motion. One of the active volumes is located at the hot end 
of the machine, one at the cold end, and the third in the crankcase-sump region. 
The remainder of the internal volume, or void volume, consists of fluid passages, 
component clearances and voids in the packed bed regenerators. 

The arrangement of the two regenerators relative to the active volumes is 
such that gas passing from one volume to the next must first pass through a 
regenerator.. In operation, gas passing through a regenerator is heated or 
cooled depending on flow direction. The regenerators are energy storage devices 
whose heat capacity greatly exceeds that of the gas. For example, gas flowing 
from the sump volume toward the hot volume is heated by the hot regenerator. 

The energy added to the gas was stored in the regenerator packing, or matrix, 
by gas flow in the reverse direction during a previous part of the cycle. The 
gas is heated to a temperature approaching that of the hot end by the time it 
exits the hot regenerator and enters the hot volume. The energy storage capa- 
bility of the regenerator is generally such that its temperature profile is 
only slightly altered by providing energy to the gas. Flow in the reverse 
direction results in cooling of the gas by the regenerator until it enters the 
sump at a temperature approaching that of the heat sink. 

The net effect of the operation of the regenerators can be viewed as isola- 
ting the three active volumes of the refrigerator in a thermal sense. Though 
this is not strictly true, the regenerators tend to maintain each active volume 
at nearly constant, but different temperature levels. These temperature levels 
are almost independent of the size of the volume, which is a function of crank- 
shaft position. Thus with the hot displacer at the bottom of its stroke, the 
major part of the gas would be in the hot volume at an elevated temperature. 

With the hot displacer at the top of its stroke, the majority of the gas would 
be in the sump volume at a temperature approaching that of the heat sink. 
Remembering that the total internal volume is constant, as is the total mass of 
gas contained, rotation of the crankshaft brings about a change in the average 
temperature of the gas and an associated change in pressure. The pressure drop 
in the regenerators is very small and thus the pressure throughout the internal 



AIRESEARCH MANUFACTURING COMPANY 


72-881*6 
Page 2-3 




volume of the machine is nearly the same at a given crankshaft position. This 
pressure depends on the amount of gas at a given temperature and therefore 
varies cyclically as the crankshaft is rotated. The only function of the drive 
motor is to provide the reciprocating motion of the displacers, not to provide 
mechanical work to the working fluid. 

With the above processes described, the VM refrigerator can be viewed as 
a Stirling-type expansion machine (cold end) driven by a thermal compressor 
(hot end). Hot displaced volume is generally large compared to the cold dis- 
placed volume. Thus, the shuttling of gas between the hot and sump active 
volumes via the hot regenerator is the major factor causing the internal pres- 
sure to vary in a cyclic manner. These pressure variations (decreased slightly 
by pressure drop) are imposed on the cold end of the machine. By proper phas- 
ing of the cold displacer motion relative to the hot displacer, the cold dis- 
placed volume can be operated as an expansion volume to produce the cooling 
effect. A 90 degree phase angle, with the hot displacer leading the cold dis- 
placer, has been shown to provide near optimum performance. 

The production of the cooling effect in the cold end of the machine can 
be viewed as follows: The gas entering the cold volume is cooled to near the 

refrigeration temperature by the cold regenerator. This gas then undergoes an 
expansion thereby further reducing its temperature and providing the capab i 1 i ty 
to absorb heat from the refrigerator load. The expansion or reduction in pres- 
sure is brought about primarily by the processes taking place in the hot and 
sump regions of the machine rather than movement of the cold displacer. The 
phasing of the cold displacer movement is set to take' maximum advantage of the 
pressure variations established within- the unit. Subsequent to the expansion 
process, the gas is then returned to the regenerator at a temperature slightly 
below the original incoming temperature. This serves to reestablish the tem- 
perature profile in the regenerator matrix by absorption of heat from the pack- 
ing and makes up for regenerator inefficiencies. The regenerator is thus 
placed in a thermal status to cool incoming gas during the next cycle process. 

PERFORMANCE AND DESIGN SUMMARY 

A summary of major performance and design parameters for the GSFC frac- 
tional watt Vuilleumier cryogenic refrigerator is given in Table 2-1. 
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GSFC FLIGHT PROTOTYPE FRACTIONAL WATT VUILLEUMI ER CRYOGENIC REFRIGERATOR 
SUMMARY OF NOMINAL PERFORMANCE AND DESIGN VALUES 
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SECTION 3 


TECHNICAL DISCUSSION 


INTRODUCTION 

This section presents a detailed discussion of all of the analyses per- 
formed leading to definition of the recommended preliminary design fractional 
watt Vuilleumier refrigerator. The characteristics of the design have been 
summarized in Section 2 of this report. This section provides the reasoning 
and supporting data leading to the selection of the particular configuration 
of each major component of the machine. 

The design requirements are presented in detail in order to completely 
define the overall objectives. Next the design approach is discussed, with 
particular emphasis on the basis of selection of each of the critical operating 
parameters and mechanical features of the machine. Under thermal design, both 
the individual components and development of optimization tradeoff parameters 
are presented. The heat pipe interface is discussed, including the heat pipes 
themselves, the interface clamp, and the sump water cooling coils. The long 
life aspects of the refrigerator, as well as stress analysis and interfacing 
requirement details, are presented under the subsection entitled mechanical 
design. Finally, the specialized features of the refrigerator drive motor are 
discussed in detail. •• 

DESIGN REQUIREMENTS 

The flight prototype fractional watt Vuilleumier cryogenic refrigerator 
has been designed to meet the requirements of NASA/GSFC Statement of Work 
APB- 763 - 28 , dated September 15. 1971* The detail requirements are summarized 
in Table 3"1 • In addition to the detail requirements specified, the refriger- 
ator is to interface with a radiometer device. Considerable attention has 
been paid to the interfacing requirements, both through close coordination 
with NASA and the radiometer manufacturer, and by detail design of various 
candidate interface component concepts. Thus not only the detail requirements 
have been considered, but also the intent of the program. This summary is 
presented in order to assure clear understanding of all program requirements. 

DESIGN APPROACH 

The general approach taken in the design of the VM refrigerator has been 
one of achieving all of the specified requirements. The power specified has 
dictated minimization of all losses. Considerable effort has been expended in 
developing methods of joining the various components of the refrigerator with 
a minimum of interface weight. The primary objectives, however, have been 
absolute assurance of the required refrigeration capacity and long life. 
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TABLE 3-1 


DESIGN REQUIREMENTS 
FLIGHT PROTOTYPE 

FRACTIONAL WATT VU I LLEUMI ER CRYOGENIC REFRIGERATOR 


Parameter 

Requ i remen t 

Thermal 


Refrigeration Heat Load 

0.25 watt 

Ref r i gera t i on Temperature 

65°K 

Thermal Input 

80 watts maximum 

Heat Rejection Temperature 

333°K maximum 

Operating Life 

Two years minimum at 100% duty cycle, 
unattended in a spacecraft environ- 
ment 

Mechanical 


Configuration 

In-line 

Weight, Including Motor 

18 lb maximum 

Heat Rejection Mechanism 

Water cooled sump interface, compatible 
wi th heat pi pes 

Heat Supply Mechanism 

Electrical heater conductively inter- 
faced with the hot cylinder 

Size 

Not specified 

Electric Motor Power 

10 watts maximum 

Envi ronmental * 


Shock 

Not specified 

Humid i ty 

90% ±3% at 30°C ±2°C for 2 hours 

V i b ra t i on 

3 axis, sinusoidal and random motion 

(a) S inusoi dal -Sweep Frequency 

See Table 3“2 

(b) Random Motion 

See Table 3~3 

P rel aunch 

24°C ±2.8°C, 50% relative humidity 
maximum, for 3 weeks 

Sand and Dust 

Not speci f ied 

Design Load Factors 

30 g i n any d i rect ion 

Storage Temperature 

-17.8°C to +54.5°C up to 3 years 

Storage Pressure 

1 atm to 10 10 torr up to 3 years 

Att i tude 

Operational or transportable in any 
att i tude 

Acce 1 e rat i on 

Operable under 5 g perpendicular to 
ini ine axis 

Tempe rature 

Operable in a vacuum at 0 to 1 20°F 
steady state or cycling between 0 
and 120°F with a minimum cyclic 
period of 6.2 hours. 

Operating Pressure 

Operable in ambient pressure from 

1 atmosphere to 10 10 torr for five 
years 


•'Testing to confirm compliance with these requirements is not required. 
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TABLE 3-2 


DELTA SINUSOIDAL VIBRATION SCHEDULE 
SWEEP APPLIED FREQUENCY ONCE THROUGH EACH RANGE 


THRUST AXIS Z-Z 

Frequency 

Leve 1 

Sweep Rate 

Range (Hz) 

(O-to-peak) 

(Octaves/mi nute) 

5-11 

12.2 mm. d.a. 

2.0 ! 


(0.48 in. d.a.) 


1 1-17 

±2.3 g 

2.0 

17-23 

±6.8 g 

1.5 

23-200 

±2.3 g 

2.0 

LATERAL AXES X-X AND Y-Y 

5-9.9 

19.0 mm d.a. 

(0. 75 in. d.a.) 

2.0 

7.1-14 

±2.0 g 

2.0 

14-200 

±1.5 g 

2.0 


TABLE 3-3 

DELTA RANDOM VIBRATION 
APPLIED IN ALL THREE AXES 


Frequency 
Range (Hz) 


A cce lerat ion 
(g - rms ) 

Duration 

20-300 

0.0029 to 0.045 
increasing from 
20 Hz at rate of 
3dB/octave 

9.1 

4 mjnutes 
ax i s 

each 

300-2000 

0.045 

9. 1 

4 minutes 
axi s 

each 
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In the process of synthesizing the design of the refrigerator, many 
tradeoffs were made in arriving at the preliminary design recommended in this 
report. These tradeoffs were made based on quantitative evaluations as well 
as qualitative considerations. Much of the experience gained during the GSFC 
5 watt refrigerator development program has been applied directly in the 
evolution of this design. 

Details of each design area and its design selection are discussed in 
another portion of this section; however, an overall summary of major design 
considerations is presented in the following paragraphs. 

Temperatures 

The selection of the design temperatures was based on system requirements 
of the overall refrigerator system. The cold end cooling temperature of 65 K 
(117 R) was specified as a design requirement; by calculating the internal gas 
heat transfer coefficients, the average gas temperature for the cold expansion 
volume was found. 


The gas temperature of the sump is primarily controlled by the heat rejec 
tion system. For this system, the maximum coolant temperature at the sump is 
specified as 333 K (600 R) . An average sump gas temperature of 344 K (620 R) 
was selected by adding the temperature drops across the sump heat exchanger, 
the interface between the sump and the sump clamp, the clamp itself, the inter 
face between the clamp and the water cooling coil, and the film temperature 
drop between the coil and the cooling water. 


The selection of the hot end temperature was based primarily on material 
considerations. With the use of electrical heating, the power source tempera- 
ture does not play an important part in the selection of the hot end tempera- 
ture. Heaters are available which are capable of operation at temperatures 
considerably above the maximum operating temperature of the hot end pressure 
containment dome. From a thermodynamic standpoint, the higher the hot end 
temgerature^ the higher the thermal efficiency. The selected temperature of 
867 K (1100 F) was based on the strength characteristics of Inconel 718 and 
also on the bearing characteristics of Boeing Compact 6-84-1. Figure 3“1 
shows the sudden drop in strength of Inconel 718 at temperatures above 992°K 
(1200°F) . A 55 .6°K (100°F) design margin of safety was allowed below this 
value. The hot end is designed structurally for pressure containment of 992 K 
(1200 F) ; thus the use of 867 K (1100 F) hot end temperature represents a margin 
of safety on performance, and not stress level. Based on heat transfer analysis 
of the hot end heat exchanger, the average gas temperature in the hot expansion 
volume is 853°K ( 1 535°R) - 


Weight 

As previously indicated, the weight of the refrigerator was considered in 
detail in the design. Minimum weight materials have been utilized wherever 
possible. An example of this is the sacrifice of thermal performance incurred 
by the use of aluminum instead of copper for the sump interface clamp. The 
weight of the moving elements of the refrigerator has been minimized in order 
to reduce both the induced dynamic loads and total weight. The assembly con- 
cepts have been developed with minimum weight a primary consideration. 
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STRESS, KSI 



Figure 3-1 • Ultimate Strength of Inconel 718 



AIRESEARCH MANUFACTURING COMPANY 


72-8846 
Page 3~5 





Power 


Power was also a major concern in the design of the refrigerator. In the 
process of synthesizing the design, considerable attention has been focused on 
losses. The thermal power input was minimized by requiring very low thermal 
losses. The motor power requirement has been kept low by maintaining a very 
low pressure drop across the hot displacer. 

Pressure 


The cycle pressure was selected in a somewhat qualitative manner. Thermo- 
dynamically, the higher the operating pressure, the higher the refrigeration 
capacity (or the smaller the size of the refrigerator). A maximum operating 


6 2 

pressure of 10.33 x 10 N/m (1500 psia) was selected in trying to minimize the 
effect of axial conduction losses for the cold end. For low pressures, the 
wall thickness of the cold end is determined by fabrication limitations. How- 
ever, as the pressure increases, the wall thickness requirement remains con- 
stant until pressure level, .rather than fabrication requirements, becomes the 
design consideration. By thermally designing the refrigerator to operate at 


6 2 

6.895 x 10 N/m (1000 psia) and by structurally designing the refrigerator for 

6 2 

10.33 x 10 N/m (1500 ps.i ) , the pressure of the working gas can later be 

6 6 ? 

increased by 50 percent 6.895 x 10 (1000) to 10.33 x 10 N/m (1500 psia) , 

if necessary, to increase the capacity of the engine. This provides a margin 
of almost 50 percent in the thermal refrigeration capacity. The maximum operat- 
ing pressure was multiplied by 1.67 to determine the proof pressure, and by 
2.25 to obtain the design burst pressure. 


Cycle Speed 

The speed of the refrigerator was selected based on the extensive investi- 
gations performed in support of the GSFC 5 watt VM refrigerator. At very low 
speeds, the residence times of the gas in the heat transfer devices are such 
that near isothermal operation can be achieved. This would, however, require 
a refrigerator of relatively large size, and the other thermal losses, such as 
axial conduction, would eliminate any gains in thermal performance. If, on the 
other hand, the refrigerator was operated at a very high speed, the temperature 
response of the gas and heat exchanger devices would be such that the thermal 
performance of the refrigerator would be greatly decreased. 

A speed of 400 rpm was selected as a good compromise. This speed also 
results in a very low bearing speed and low dynamic loading for the bearings 
to support. To further increase life, it may be possible to reduce the speed, 
with an accompanying reduction in capacity. 

Mechanical Assembly 

The approach taken here has been the use of mechanical joints which allow 
the refrigerator to be disassembled after checkout tests. Upon verification of 
adequate performance, all joints may be permanently closed by welding without 
disturbing the mechanical integrity. Temporary non-organic seals are utilized 
at all assembly interfaces. Provisions have been incorporated for destruction 
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In this manner, the attainment 
assured. Many mechanical closure 
on attainment of the goals 


THERMAL OES I GN 
I n t roduc t ion 


of these seals prior to final joint closure, 
of absolute leak tight final closure welds is 
concepts were investigated, with the emphasis 
described above with minimum weight. 


The thermal design of a VM refrigerator is a lengthy i terat i ve process. 
Initially, rough cut design calculations are made in order to establish the 
design approach and basic sizes of the machine's elements. After the basic 
design is defined, effort must be concentrated on matching the thermodynamic 
design with that of the heat transfer devices (heat exchangers and regenera- 
tors). Typically, the configurations and volumes of the heat transfer devices 
are adjusted to improve their heat transfer and pressure drop character i s t i cs . 
These adjustments imply that changes be made to the active displaced volumes 
to compensate for the influence of the heat transfer devices on the thermo- 
dynamic processes of the working fluid. In turn, once the active volumes are 
changed, the heat transfer devices require adjustment to account for the vari- 
ations in flows, pressure levels, and heat loads. This iterative process is 
carried out until the thermodynamic cycle parameters match the designs of the 
heat transfer devices. By examining several matched designs, a near-optimum 
refrigerator can be selected. 

In this program, emphasis has not been placed on complete optimization of 
the refrigerator with respect to thermal performance; long operational life 
and reliability were considered the most important. It is believed, 
however, that due to the care taken in the detail thermal design, a near-optimum 
thermal design has resulted. 

This Task I report presents the preliminary design of the fractional watt 
VM . As detail design progresses, mod i f i cat ions to the preliminary design com- 
ponents will be made as necessary. Where possible, the experience gained in 
the design and fabrication of the GSFC 5 watt VM refrigerator has been applied 
to the configurations presented herein. Based on this background, the design 
changes necessary as detail design progresses should be minimized. This section 
summarizes the design analyses; the topics covered are: 

• Thermal Design Summary 

• Cold-end regenerator 

• Hot-end regenerator 

• Cold-end heat exchanger 

• Hot-end heat exchanger 

• Ambient heat exchanger 

• Cold-end insulation 
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Hot-end insulation 


• 

• Pressure drop and dead volume tradeoff factors 

• Cold-end flow distributor 

• Design of cold end seal 

• Design of hot end seal 

The general method of analysis applied to both the hot and the cold 
regenerators, and the regenerator matrix characterization is presented in 
Appendix A of this Task I report. 

Thermal Design Summary 

The thermal design of the VM refrigerator is a lengthy iterative process; 
therefore, the cycle and thermal design are summarized before the detailed 
discussions. Many of the thermal features of the preliminary design VM have 
been presented in Section 2 of this report. This sub-section gives details 
of the thermodynamic cycle' and the breakdown of the hot end power and cold end 
losses. 

1 . Thermodynamic Cycle 

The thermodynamic cycle is analyzed by use of the ideal cycle analysis 
computer program. The information obtained includes the cyclic pressure pro- 
file, the mass flows in various portions of the machine, the ideal refrigera- 
tion, and the heat input required. Figure 3-2 presents the program output for 
the operating conditions of the preliminary design fractional watt VM refrigera- 
tor. The effects of pressure drop on the hot and cold-end have been accounted 
for in the calculations. Table 3“ 4 gives the nomenclature utilized in Figure 3-2. 

It should be noted that the output from the ideal cycle analysis program 
forms a major portion of the input for the regenerator analysis. This analysis 
procedure is presented in Appendix A. 

2. Hot End Power 


The hot end power input is specified at 80 watts. This power input is 
divided into several portions: (1) the ideal required input, from Figure 3 - 2 , 

(2) the losses to ambient, and (3) losses to the sump region. The losses have 
been calculated based on the preliminary design configuration of the VM 
refrigerator. The calculated total is slightly less than the maximum allow- 
able, which provides a margin of safety for changes incurred as detail design 
progresses. The calculated hot end power distribution is presented in 
Table 3~5- 



AIRESEARCH MANUFACTURING COMPANY 


72-8846 
Page 3-8 



FRACTIONAL WATT VM-PRELtMt NAR y DESIGN UNIT WITH 6STIMATE0 DELTA P 
OPERATING PARAMETERS 



COLD VOLUME TEMP, » 

112.00 


SUMP VOLUME TEMP, • 

620.00 


HOT VOLUME TEMP, ■ 

1533.00 


COLO REGEN, TEMP, • 

366.00 


HOT REGEN, TEMP, • 

1077,30 


COLO DJSPLACEO 

VOU. » 

.05530 


HOT DISPLACED 

you. ■ 

1,88500 


COLD OEAD VQL, 

■ 

,10000 


SUMP OEAD VOL. 

9 

3.00000 


HOT DEAD VOL, 

u 

.33000 


COUP REGEN, yOU. ■ 

1,32800 


HOT REGEN, VOU 

. « 

3,33300 


GAS CONSTANT 

* 

6634.40 


SPEED 

■ 

400,00 


CHARGE PRESSURE » 

753,41 


CHARGE TEMPERATURE ■ 

535,00 


MASS OF FLUID 

• 

,0029 


total volume 

» 

10.05130 

PRtSSUKE 

- MASS - FLOW 

profile 


angle 

phess ye 

VA 

VH 

DEG 

PS I A CU-IN 

CiMN 

CU-IN 

21) . 

969.93 .1017 

3.6738 

1,6146 

40, 

986.03 .1063 

3.3836 

1,6982 

60. 

996.63 .1138 

3.1678 

2,1086 

80. 

999.98 ,1228 

3,0468 

2.2206 

100. 

995.56 ,1324 

3,0371 

2.2207 

120. 

984.08 ,1415 

3.139? 

2,1069 

140. 

967.36 .I486 

3.3428 

1,8967 

160. 

947.82 ,1336 

3,6214 

1,6153 

180. 

928.03 .1553 

3,9419 

1.2931 

200. 

910.36 ,1336 

4.265^ 

.9707 

220. 

896.38 .1468 

4,5543 

,6872 

240. 

888.01 .1415 

4,7722 

,4766 

260. 

885,33 ,1323 

4,8934 

.3645 

280. 

888.84 .1229 

4,9033 

,3642 

iOO. 

898,18 ,U38 

4,8006 

,4738 

->20. 

912.58 • 1063 

4,597? 

.6859 

340. 

930.67 ,1017 

4, 3193 

.9692 

360. 

950.56 .1000 

3,9989 

1,2914 

THE CYCLE PRESSURE ABOV 

£ IS IDEAL. THE P 

hot maximum op* 1,0000 

psj cold maxi 

IDEAL REFRIGERATION ANO 

HEAT INPUT 

REFRIGERATION * 

3.3313 

WATTS 

thermal heat ■ 

27,0030 

WATTS 

MAX. 

PRESSURE » 1000,0000 

PSIA AT 


• • '.i 

HN 

J-IH 


’IN 


MOOTC 

LB/SfC 

.00102 
.00194 
.001?* 
•00193 
•00172 
•00126 
• 000(6 
•00003 
••00099 
••00110 
• i 001(7 
••00166 
••00171 
••00197 
••00127 
••00061 
••00024 
.00039 


MDOTA 

LB/SEC 

••00999 

■,00848 

-.00614 

-.00286 

.00082 

.00430 

.00706 

•00877 

.00932 

.00879 

.00734 

.00919 

,00299 

••00034 

-.00326 

-•00994 

*.00810 

••00941 


MOOTH 

LB/SEC 

.00931 

. 0 Q 462 

.00304 

•00099 

.00120 

.00317 

.00463 

.00942 

.00992 

•00498 

.00394 

.00293 

.00089 

.00064 

.00292 

.00399 

• 009 Q 9 

.00969 


HDOTRCA MOOTRHA 
LB/SEC LB/SEC 


.00264 

.00277 

.00292 

.00190 

.00101 

.00000 

. 000*7 

.00176 

.00230 

.00294 

.00290 

.00220 

,00169 

.00100 

.00020 

.00069 

.00146 

.00218 


,00695 
,00971 
,00362 
,00096 
-.00183 
» , 00430 
-.00609 
•,00701 
•,00703 
• , 00629 
-.00484 
.,00299 
-, 00087 
,00134 
,00346 
,00929 
,00662 
,00723 


PSI 


ANGLE • 78,59 DEGREES 


Figure 3-2. Data Output From Idealized Cycle Program 
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TABLE 3-k 


IDEAL VUILLEUMIER CYCLE ANALYSIS 
NOMENCLATURE KEY 


Symbol 

Press 

Angle 

VC 

VA 

VH 

MDOTC 

MDOTA 

MDOTH 

MDOTRCA 

MDOTRHA 


Def i ni tion 


Cycle Pressure 


Crankshaft angle referenced to 
cold displacer top dead center 


Cold displaced volume 


Ambient displaced volume 


Hot displaced volume 


Flow rate into cold volume 


Flow rate into ambient volume 


Flow rate into hot volume 


Flow rate into cold regenerator at 
the end toward the sump 

Flow rate into hot regenerator at 
end toward the sump 
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TABLE 3-5 


HOT END POWER REQUIREMENTS 



Heat Input, 
Watts 

Regenerator Loss (Jhc dT) 

5.67 

Regenerator Inner Wall Conduction Loss 

6.34 

Regenerator Bed Conduction Loss 

3.71 

Regenerator Outer Wall Conduction Loss 

12.82 

Displacer Conduction Loss 

15.33 

Insulation Loss 

5.90 

Seal Leakage Loss 

0.54 

Ideal Heat Input (Includes AP Correction) 

27.00 

Total Calculated Heat Input 

77.31 watts 


3. Cold End Losses 


The cold end losses are calculated in a similar manner to the hot end. In 
this case, however, the calculated losses represent heat introduced to the cold 
end in addition to the refrigeration heat load. Thus the sum of the losses and 
the refrigeration requirement represents the total heat that must be carried to 
the sump region and rejected to the heat sink. These losses are summarized 
in Table 3-6. It should be noted that the displacer gas. leakage loss is based 
on very conservative assumptions and bearing clearances representative of two 
years of operation. The bearing wear rate and thus the clearances after two 
years are also considered conservative. Thus it is expected that the refrig- 
erator will be capable of somewhat greater net cooling than indicated in Table 
3-6 after two years of operation. 

Cold Regenerator Design 

1 . I n t roduct i on 


The cold regenerator is one of the most important components of the VM 
refrigerator, it not the most important. The regenerator cools gas as it flows 
from the sump region of the machine to the cold expansion volume. The gas 
expansion process further reduces the temperature, providing cooling at the 
cold temperature. After heat absorption from the refrigeration load, the gas 
is returned to the sump through the regenerator. This reverse flow process 
heats the gas by removal of energy stored in the regenerator matrix, and 
reestablishes the matrix temperature profile for cooling the gas during the 
next cycle. The periodic storage and removal of energy from the matrix allows 
the working fluid to pass from one essentially constant temperature portion of 
the machine to another. An ideal regenerator may thus be viewed as a tempera- 
ture isolator. 



AIRESEARCH MANUFACTURING COMPANY 

I iii Mm. v. Uilit m m.i 


72-8846 
Page 3-11 









TABLE 3-6 

SUMMARY OF COLD END LOSSES 



Heat Loss , 
Watts 

Regenerator Loss C/flc dT) 

0.929 

Regenerator Inner Wall Conduction Loss 

0.200 

Regenerator Bed Conduction Loss 

0.65^ 

Regenerator Outer Wall Conduction Loss 

0.629 

Displacer Conduction Loss 

0.256 

Insulation Loss 

0 .0^7 

Displacer Leakage Loss (2 year worn bearings) 

0.569 

Total losses 

3.284 

Ideal Refrigeration (Includes AP Correction) 

3.551 

Ideal Q Minus Losses = Net Refrigeration 
(after 2 years operation) 

0.267 (0.25 req'd) 


The design requirements for an efficient regenerator demand that it must 
(1) absorb heat from the gas stream while at nearly the same temperature as 
the gas, (2) store this energy without significant temperature change in a 
given matrix locality, and (3) resupply the energy to the gas stream when the 
flow reverses direction, again while at a temperature very closely approaching 
that of the gas. These design requirements dictate that the regenerator packin 
(1) must have a very large heat capacity relative to that of the gas, (2) must 
have large values of heat transfer coefficient, heat transfer area, and thermal 
diffusivity, and (3) must be configured to limit axial conduction of heat from 
one end to the other. In addition, for the VM application, the void volume in 
the regenerator must be minimized. Obtaining these desirable regenerator 
characteristics from a specific packing generally leads to increased pressure 
drop as the thermal characteristics improve. Since pressure drop as well as 
void volume degrades the overall refrigerator performance, the basic tradeoff 
in the design of a regenerator is between heat transfer potential and the 
detrimental factors accompanying this potential. 

2 . Method of Analysis 

The analysis of the regenerators for a VM refrigerator cannot be based on 
the classical effectiveness parameters that make use of end point temperatures. 
The system pressure fluctuates and a considerable amount of gas is stored in 
the regenerator void volume during various parts of the flow cycle. This 
periodic mass storage characteristic, coupled with the basic transient nature 
of regenerators, dictates that a finite difference technique be used in the 
anal ys i s . 
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AiResearch has developed a computer program which utilizes finite difference 
techniques to analyze regenerators. The required program inputs are the regene- 
rator physical characteristics, the matrix heat transfer and pressure drop 
characteristics, and the initial and boundary conditions. The computer program 
and the general characteristics of both spherical shot and screen regenerator 
packing materials are discussed in detail in Appendix A of this report. 

The required boundary conditions are the time dependent pressures and flow 
rate profiles and the gas temperature at the sump end of the regenerator. The 
cold end gas temperature under reverse flow conditions is also required. These 
parameters are obtained from the design requirements of the VM, and the idealized 
cycle analysis computer program. 

3. Initial Iterations 


The regenerator design concept studies performed on the 5 watt VM program 
(Ref. 3-1 ) clearly indicated the advantages of the annular regenerator, as 
compared to those located inside the displacers or those located in parallel 
with the cold displacer. Therefore, the annular regenerator was the only type 
considered for the present machine. 

The cold regenerator design was evolved along with the rest of the VM 
through numerous design iterations. However, certain initial ground rules were 
utilized which simplified the process. The cold displacer diameter was fixed 
early in the program, thus defining the inner diameter of the annulus. It was 
decided that, if at all possible, the packing would be composed entirely of 
spherical shot, in order to provide the maximum matrix heat capacity. The 
ground rule was also established that the annular height would be 30 shot 
diameters. This is necessary in order to prevent channeling of the gas along 
the regenerator wall. 

The first regenerator analyzed utilized a packing material composed of 
0.000178m (0;007 in.) diameter monel shot. The annular height was 0.00566m 
(0.21 in,), and the regenerator length was 0.1132m (b.S in.). The gas and metal 
temperatures at the cold end are presented in Figure 3-3 as a function of crank- 
shaft angular position. Figure 3"^ presents the pressure drop characteristic 
of this regenerator. 

Figure 3 - 3 indicates that the heat transfer coefficient between the gas 
and the spherical shot matrix is sufficiently high. The maximum temperature 
deviation between the gas and Q the metal, except at the instant of flow reversal, 
is approximately 0.167°K (0.3 R) . However, the heat capacity of the matrix is 
insufficient. This is evidenced by the gas temperature deviation of 1 . 333°k 
(2.A°R) , which is considered excessive. The maximum pressure drop of 

4 2 

2.76 x 10 N/m (*4.0 ps i ) , as shown in Figure 3 _ *t, is also considered too high. 

The regenerator deficiencies pointed out above led to the next design 
iteration. The first (warmer) half of the regenerator was filled with 0.00025**m 
(0.010 in.) diameter monel shot, and the annular height was increased to 0.00762m 
(0.3 in.). Both of these changes reduce pressure drop and increase matrix heat 
capacity. The second (colder) half of the regenerator was packed with 0.000178m 
(0.007 in.) diameter shot, in order to maintain the high heat transfer coeffi- 
cient. 



AIRESEARCH MANUFACTURING COMPANY 

Imr.init; ( .tlr'.iff... 


72-88*46 
Page 3-13 





c 

o 


03 

3 — 

<D 


t/) 

l— 


"D 

C 

UJ 


o 

o 


c 

o 


u 

o 


(D <D 
C 

3 <U 
■M CD 
03 0) 

s- OT 

cu 

E" ^ 
<D O 
H* O 


rr\ 

I 

ro 

0) 

L_ 


3 

cr> 


y Q ‘ 3dfll\/a3dW31 



AIRESEARCH MANUFACTURING COMPANY 


72 - 88^6 
Page 3-14 



o 



■J* CN o CN ' kO 

-r | III 


isd ‘doaa 3dnss3dd 


lT\ 

r-- 

r-» 


CO 

UJ 

UJ 

CC 

CD 

UJ 

Q 


CD 

2 

< 


Ll. 

< 


< 

CC 

<_> 


O 

4 — < 

«J 

U 

0) 

c 

<u 

cn 

Q) 

oc 

“U 

o 

o 

c 


03 

<u 


i- 


Q_ 

O 


<D 

u 

Z3 

</> 

c/1 

Q) 

c. 

Q_ 


I 

c+\ 

<D 


L_ 

D 

cn 



AiRESEARCH MANUFACTURING COMPANY 


72-8846 
Page 3“ 1 5 






The metal and gas temperatures at the cold end of this second iteration 
regenerator are presented in Figure 3 - 5. Figure 3-6 presents the gas pressure 
drop characteristic. The maximum temperature deviation is reduced to 0.855°K 

o 4 2 

(1.54 R) , and maximum pressure drop is 1.048 x 10 N/m (1.52 psi). These 

values represent considerable improvement over the first design iteration. 
However, when all losses were considered, it was decided that the remaining 
net refrigeration with this configuration was too marginal for this early stage 
in the program. The regenerator was accordingly lengthened, to reduce conduc- 
tion losses, thus leading to the selected preliminary design. 

4. Selected Design 

The cold regenerator selected for the preliminary design fractional watt 
VM is of an annular configuration, packed with two sizes of spherical monel 
shot. Monel was chosen as the material because of its superior heat capacity 
in the temperature range of interest (see Appendix A) . The total length is 
0 . 127 m (5 in.), and the annular height is 0 . 00762 m ( 0.3 in.), 30 diameters of 
the largest shot, in order to prevent gas channeling along the walls. The 
first 0 . 0762 m (3 in.) from the sump end are packed with 0.000354m ( 0.010 in.) 
diameter shot; the remaining . 0508 m (2 in.) are packed with 0 . 000178 m ( 0.007 in.) 

2 2 

diameter shot. The frontal area is 0. 00044m (0.682 in. ). 

The first section of the regenerator, 0.000254m (0.010 in.) dia shot, has a 

2 3 2 3 

porosity of 39 percent, an area to volume ratio of 14,400m /nr (366 in. /in. ), 
and a hydraulic diameter of 0 . 0001086. (0.004274 in.). For the second section 
0 . 000178 m ( 0.007 in.) dia. shot, the porosity is 39 percent, the area to volume 
2 3 2 3 

ratio is 20 , 600 m /m (523 in. /in. ), and the hydraulic diameter is 0 . 0000759 m 
( 0.002988 in.) . 

The larger diameter shot is used in the warm end of the cold regenerator 
to reduce the pressure drop, but still maintain a large thermal capacity. The 
smaller diameter shot at the cold end yields the required high heat transfer 
coef f i c i ent . 

Table 3 _ 7 presents the detailed output from the regenerator analysis com- 
puter program for this preliminary design cold regenerator. The nodewise pres- 
sures and temperatures of the gas are listed as a function of the angular 
position of the crankshaft, or time. Angular displacement is referenced to. the 
top dead center position of the cold displacer. The output parameters are 
matrix temperature, gas temperature, gas density, gas pressure, and gas flow 
rate. Positive gas flow rate denotes mass flow from the sump end of the 
regenerator toward the cold end. Node 0 represents the sump end of the regene- 
rator, and Node 13 the outlet face at the cold end. Nodes 1 through 7 are the 
portion of the core packed with 0.000254m (0.010 in.) dia shot, and Nodes 8 
through 12 represent the 0 . 000178 m ( 0.007 in.) dia shot. 

Figures 3 - 7 through 3"9 present plots of key parameters from the computer 
program printout. The matrix and gas temperatures at the cold end of the 
regenerator are plotted in Figure 3~7- The small difference between the two 
temperatures is indicative of excellent heat transfer, and the 0.889 K (1.6 R) 
temperature swing of the gas is indicative of adequate heat capacity. 
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TABLE 3-7 


PRELIMINARY DESIGN 
COLD REGENERATOR 
PERFORMANCE CHARACTERISTICS 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 

CD 

II 

°o 


REGULAR 

PRINTOUTS 






DATE 

■ 27 SEP 72 TIME ■ 17119133 

ME( SEC , 

) * 4.49989-01 


4722 

« COUNT (NO, OF 

CALCULATIONS) 

M 

TM(N) 

TG(N> 

RG(N) 

PG(N> 

WG ( N ) 


OEG.R 

OEG.R 

LBM/CF 

PS! A 

LBM/SeC 

-o 

S. 19488.02 

6.20000*02 

3,33439-01 

9,50334*02 

?, 17912-03 

1 

5 . 96764*02 

5.97083*02 

3,73384-01 

9,50408*02 

2.13282-03 

2* 

5.56808*02 

5.57110*02 

6,14103-01 

9.50185*02 

2.03389-63 

3 

5 . 14275*02 

5.14913*02 

6,62202-01 

9.49997*0? 

1.92747-03 

4 

4.71956*02 

4.72609*02 

7.19777-01 

9.49841*02 

1.81214-63 

5 

4 . 29641*02 

4.30309*02 

7.86209-01 

9.49716*02 

1.68630-63 

6 

3.87336*02 

3.88019*02 

8,68513-01 

9.49617*02 

1.54818-63 

7 

3.45036*02 

3.45736*02 

9.69237-01 

9.49542*0? 

1,39484-63 

8 

3.02796*02 

3.03211*02 

1,09818*00 

9,49481*02 

1.22221-63 

9 

2.60416*02 

2.60828*02 

1.26766*00 

9.49419*0? 

1.03780-63 

10 

2.18133*02 

2.18563*02 

1.49931*00 

9.49380*02 

8.22114-64 

11 

1.76017*02 

1.76472*02 

1.83978*00 

9.49360*02 

5.62211-64 

12 

1,33948*02 

1.36472*02 

2.33832*00 

9.49332*0? 

2,38966-04 

13 

1.12002*02 

1.12119*02 

2.80969*00 

9.49331*02 

4.88867-05 


o 

o 

ca 

II 

CD 



DATE 

* 27 SEP 72 TIMF « {7il9l5l 

IMECSEC.) * 4,62479- 

31 

4868 

■ COUNT (NO. OF 

CALCULATIONS) 

N 

n<N> 

TGCN) 

RG(N) 

PG(N) 

WG( N ) 


DEG . R 

OEG.R 

LBM/CF 

PSIA 

LBM/SEC 

-0 

6,19639*02 

6.20000*02 

5,69082-01 

9,78303*02 

2,73966-03 

t 

5.96980*02 

3.97310*02 

5.89314-01 

9.78136*02 

2.7Q238-63 

2 

5.36900*02 

5.37316*02 

6.30*66-01 

9.77834*02 

2.62275-63 

3 

5.14683*02 

3.13339*02 

6,80132-01 

9,77573*0? 

2.53711-63 

4 

4.72363*02 

4.73039*02 

7,39064-01 

9,77350*0? 

2.44433-63 

5 

4.30046*02 

4.30740*02 

8.09166-01 

9.77163*0? 

2.34311-03 

6 

3.87737*02 

3.88430*02 

8.91340-01 

9.77010*02 

2.23208-63 

7 

3.43437*02 

3.46170*02 

9,94387-01 

9.76887*02 

2.10*83-63 

6 

3.03227*02 

3.03664*02 

1.12625*00 

9,76780*02 

1. 97016-63 

9 

2. 60*41*02 

2.61280*02 

1.29933*00 

9.76658*02 

1,82213-63 

10 

2.18580*02 

2.19038*02 

1,53568*00 

9.76571*02 

1,64918-03 

11 

1.76502*02 

1.76979*02 

1.87815*00 

9.76514*02 

1,44108-03 

12 

1 . 36511*02 

1.36988*02 

2.38827*00 

9.76479*02 

1.18285-63 

13 

1 . 12340*02 

1.12594*02 

2.86844*00 

9.76469*0? 

1.03120-63 


NOTES: 


1. TEMP ,°K m °R/1 . 8 

2. DENSITY, kg/m 3 - 16.02 X LB/FT 3 

3. PRESSURE, N/m 2 - 6.895 X 1 0 3 X PS I 

4. FLOW, kg/sec = 0.4536 X LB/SEC 
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TABLE 3-7 (Continued) 


NODE 

MATRIX 

GAS 

GAS 


GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 


PRESSURE 

RATE 

9 = 60° 


REGULAR 

PRINTOUTS 







OATE 

• 

27 SEP 72 TIME < (7128111 

I SEC. 

) * 4,74931-01 


5013 

• ■ 

COUNT <N0, OP 

CALCULATIONS) 

N 

TM(N> 

TCI N) 

RSI N) 


PCI N) 

WCIN) 


oeg.r 

OEG.R 

LBM/CF 


PSti 

L0M/SEC 

-0 

6.19917* 0 2 

6.20000*02 

5.79333-01 


9.96553*02 

2.52277-83 

1 

5.97(83*82 

5.97460*02 

5,99985-01 


9.96406*02 

2.50146*53 

2 

5.57769*02 

5.57781*02 

6.41923-01 


9.96137*02 

2.45595-63 

3 

5.15069*02 

5.15622*02 

6,92012-01 


9.95900*02 

2.40700-63 

4 

4.72755*02 

4.73325*02 

7,51861-01 


9.95694*02 

2.35400-63 

5 

4.30438*02 

4.31029*02 

8,23126-01 


9,95518*02 

2.29619-63 

6 

3. 88{3i*02 

3.86742*02 

9.06569-01 


9.95370*02 

2.23278-03 

7 

3.45837*02 

3.46466*02 

1.01118*00 


9.95247*02 

2,16244-83 

8 

3.03664*02 

3.04042*02 

1.14477*00 


9,95136*02 

2.08333-63 

9 

2,61285*02 

2.61673*02 

1,32017*00 


9,95004*02 

(, 99893-53 

10 

2.19059*02 

2.19467*02 

1.55929*00 


9.94905*02 

1.9oo38-63 

11 

1.77641*02 

1.77471*02 

1,90494*00 


9.94833*02 

1.78199-63 

12 

1.37(54*02 

1.37588*02 

2.41715*00 


9,94785*02 

1, 63545-63 

13 

1.12798*02 

1.13047*02 

2.90462*00 


9.94769*02 

(.54938-63 

e = 90° 



OATS 

s 

27 SEP 72 TIME « i 7 i 2h 1 29 

[ ME ( SEC . : 

1 a 4.87426-01 


5147 

• 

COUNT InO, of 

calculations) 

N 

TMCN) 

TG(N) 

RGIN) 


PGINj 

toGI V) 


oec.R 

OEC.R 

LBM/CF 


PSlA 

L0M/SEC 

-o 

6.19956*02 

6.20000*02 

5 , 8q38(*01 


9.98419*02 

(.47802-63 

1 

5.975o’*02 

5.97433*02 

6.01148-01 


9,98348*02 

1.48822-q3 

2 

5.57494*02 

5.57712*02 

6,43289-01 


9.98213*02 

1.50999-63 

3 

5.15312*02 

5.15560*02 

6.93549-01 


9.98088*02 

1.53340-63 

4 

4.73004*02 

4.73268*02 

7.53617-01 


9.97974*02 

(, 55876-63 

5 

4.30695*02 

4.30976*02 

8,25074-01 


9.97871*02 

(.58641-63 

6 

3.88395*02 

3.86697*02 

*.0875(-01 


9.97779*02 

1,61675-63 

7 

3.46(12*02 

3.46436*02 

1.0l36i*00 


9.97698*02 

i. 65041-63 

8 

3.03971*02 

3.04172*02 

1.14699*00 


9.97616*02 

1.68824-53 

9 

2.61414*02 

2.61835*02 

1.32256*00 


9.97515*0? 

1,72857-83 

10 

2.19426*02 

2.19671*02 

1.56l37*oo 


9.97430*02 

{,77565-53 

11 

1.77474*02 

1.77747*02 

1.90649*00 


9.97361*02 

1.83718-63 

12 

1.37499*02 

1.37983*02 

2.41563*00 


9,97308*02 

1.90209-83 

13 

1.13239*02 

1.13432*02 

2.90239*00 


9.97287*02 

1.94313-83 




NOTES: 

1 . 

TEMP ,°K = °R/1. 

,8 





2. 

OENSITY, kg/m 3 

- 16.02 X LB/FT 3 





3* 

PRESSURE, N/m 2 

= 6.895 X 10 3 X PSI 





4* 

FLOW, kg/sec ■ 

0.4536 X LB/SEC 
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TABLE 3-7 (Continued) 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 

9 = 120° 


REGULAR 

1 PRINTOUTS 






OATE 

a 27 SEP 72 TIME » l7i2fil42 

t i me ( sec * 

) a 4. 99999-01 


5250 

■ COUNT (NO. OF 

CALCULATIONS) 

N 

TM(M) 

TG(N> 

RC ( N ) 

PG(N) 

W6( N) 


DEC, R 

OEC.R 

L9M/CF 

P8I A 

LBM/SEC 


6.19972*02 

6.20000*02 

5.72329-01 

9.84081*0? 

8,44207-08 

1 

5.97322*02 

5.97155*02 

5.93124-01 

9.84Q8l*0? 

2.94475-05 

2 

5.57510*02 

5.57231*02 

6,35054-01 

9,84079*02 

9.21793-05 

3 

5.15340*02 

5.15103*02 

6,84771-01 

9.84075*02 

1,59634-04 

4 

4.73039*02 

4.72835*02 

7.44244-01 

9.84069*02 

2.32711-04 

5 

4.30738*02 

4.30563*02 

8.14926-01 

9.84060*0? 

3.12428-04 

6 

3.88448*02 

3.88303*02 

8,97767-01 

9.84048*02 

3.99871-04 

7 

3. 46(77*02 

3.46061*02 

1.00154*00 

9.84034*02 

4.96893-04 

8 

3.04053*02 

3.04003*02 

1 , 13291*00 

9.84016*0? 

6,05914-04 

9 

2.61723*02 

2.61704*02 

1,30644*00 

9,83988*0? 

7.22239-04 

10 

2.19566*02 

2.19572*02 

1.54284*00 

9.83959*0? 

8,58020-04 

11 

1.77664*02 

1.77701*02 

1.88375*00 

9.83931*02 

1.02111*03 

12 

1.379'67*02 

1.38034*02 

2.38610*00 

9,83905*02 

1.22279-63 

13 

1.13318*02 

. 1.13605*0.2 

2,86499*00 

9,83893*02 

1. 34116-63 


CD 

11 

U3 

o 

o 



OATE 

• 27 SEP 72 TIME a f 7 j 2n 1 55 

T I ME< SEC , ) 

« 5.12496-01 


5353 

a COUNT (NO. OF 

CALCULATIONS) 

N 

TM( N ) 

TC ( N ) 

Rf.(N) 

P C ( M } 

WG< N) 


DEC. R 

OEC.R 

LBM/CF 

PSU 

LBM/SEC 

-o 

6.19870*02 

6.19606*02 

5,57844-01 

9.57720*0? 

-1,38458-63 

1 

5.97133*02 

5.96639*02 

5,78315-01 

9.57787*02 

-1.33932-63 

2 

5.57325*02 

5.56811*02 

6.i9235-01 

9,57904*0? 

-1.24262-63 

3 

5.15(62*02 

5.14641*02 

6.67894-01 

9.57998*02 

-1. 13860-63 

4 

4.72869*02 

4.72340*02 

7.26166-01 

9.58072*02 

-1. 02568-63 

3 

4.30378*02 

4.30042*02 

7.95339-01 

9,58128*02 

-9,02898-64 

6 

3.88304*02 

3.87759*02 

8.76519-01 

9.58168*0? 

-7,67927-g4 

7 

3,46048*02 

3.45494*02 

9.78252-01 

9.58196*02 

-6,18082-04 

S 

3.03928*02 

3.03616*02 

1.10831*00 

9.58213*02 

-4,49737-64 

9 

2.61627*02 

2.61299*02 

1.27641*00 

9.58227*0? 

-2.69993-64 

10 

2.19477*02 

2.19395*02 

1.50683*00 

9,58232*0? 

-6.01729-65 

11 

(.77566*02 

1.77285*02 

1.84326*00 

9,58231*02 

1.92377-64 

12 

1.37896*02 

1.37681*02 

2.33695*00 

9.58225*0? 

5.05(90-34 

13 

1.13363*02 

1.13512*02 

2.80223*00 

9.58221*02 

6.88703-64 


NOTES: 1. TEMP ,°K » °R/1 . 8 

2. DENSITY, kg/m 3 » 16.02 X LB/FT 3 

3. PRESSURE, N/m 2 = 6.895 X 10 3 X PS I 

4. FLOW, kg/ sec « 0.4536 X LB/SEC 
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TABLE 3-7 (Continued) 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 

<D 

11 

CO 

o 

o 


REGULAR 

PRINTOUTS 






OATS i 

27 SEP 72 TINS ■ t7 i 24 113 

TIMEtSEC.) * 5.2493? 

-01 

5491 ■ 

COUNT (NO, OF 

' CALCUL AT I ONS ) 

N 

TM(N> 

TGI N J 

RG(N> 

PGINj 

WC(N) 


0EG.R 

0EG.R 

L9M/CF 

PSIA 

LgN/Sgc 

-o 

6 . 19688*02 

6.19368*02 

5,41404-01 

9.28209*0? 

-2.29608*03 

1 

5.96803*02 

5.96218*02 

5.61522-01 

9.28347*02 

-2.25073-03 

2 

5.56990*02 

5.56375*02 

6.01417-01 

9.28592*02 

-2.15383*03 

3 

5, 14*32*02 

5.14205*02 

6.48825-01 

9,28800*02 

-2.04956*03 

4 

4.72544*02 

4.71906*02 

7,05648-01 

9.28973*02 

-1.93657-03 

5 

4.30?6i*02 

4.29613*02 

7,73051-01 

9.29114*02 

-1. 81324-03 

6 

3.87995*02 

3.87336*02 

8.52224-01 

9.29226*0? 

-1 . 47784*03 

7 

3.45749*02 

3.45080*02 

9.81436-01 

9.29312*02 

•1.52746*03 

8 

3.03427*02 

3.03249*02 

1.07618*00 

9.29382*02 

-1.35845-03 

9 

2.61339*02 

2.60951*02 

1.24218*00 

9.29456*0? 

-1 .17788-03 

10 . 

2.19204*02 

2.18817*02 

1.46874*00 

9.29504*02 

-9,66729-04 

11 

1.77351*02 

1,76967*02 

1,79615*00 

9,29531*0? 

-7,12516-04 

12 

1.37716*02 

1.37507*02 

2.27769*00 

9.29543*02 

-3,97321-fl* 

13 

l.l3{3i*02 

1.12000*02 

2.76124*00 

9.29545*0? 

-2,10433-04 

e = 2 io° 



OATS i 

27 SEP 72 Tins i 1 7 1 2? 1 33 

‘IMEfSEC.) s 5.37500- 

01 

5645 * 

COUNT (NO, OF 

CALCULATIONS) 

N 

TM(M) 

TGI N) 

RG(N) 

PG( N i 

WC( N) 


0EG.R 

DEG.R 

LSM/CF 

PSIA 

LBM/SEC 

-o 

6.19479*02 

6.19174*02 

5.27374-01 

9.03100*02 

-2,54000-03 

1 

5.96429*02 

5.95680*02 

5,47148-01 

9,03263*02 

-2.50796-03 

2 

5.56604*02 

5.56022*02 

5,86124-01 

9.03557*0? 

-2,43950-63 

3 

5.14449*02 

5.13852*02 

6,32418-01 

9.03812*02 

-2.36582-63 

4 

4.72162*02 

4.71552*02 

6.87967-01 

9.04030*02 

-2,28595-03 

3 

4.29483*02 

4.29258*02 

7,53818-01 

9.04214*02 

-2.19875-03 

6 

3.87421*02 

3.86982*02 

8.31234-01 

9.04365*02 

-2.10300-63 

7 

3.45376*02 

3.44723*02 

9.28266-01 

9.04467*0? 

-1,99663-63 

8 

3.03242*02 

3.02864*02 

1.05036*00 

9.0*594*0? 

-i, 87700-03 

9 

2.60<»45*02 

2.60558*02 

1.21294*00 

9,04716*0? 

' -1.7*911-63 

10 

2.18403*02 

2.18405*02 

1.43517*00 

9,04805*0? 

-1.599*1-03 

11 

1.76934*02 

1.76543*02 

1.75673*00 

9.0*864*02 

-1. 41897-63 

12 

i . 37450*02 

1.37218*02 

2,22856*00 

9.0*901*0? 

•i . 19506*03 

13 

1 . 12254*02 

1.12000*02 

2,69740*00 

9.0*912*0? 

-1,06237-03 




NOTES : 1 . 

TEMP , °K - °R/1. 

8 




2. 

DENSITY, kg/m 3 

- 16.02 X LB/FT 3 




3. 

PRESSURE, N/m 2 

- 6.895 X 10 3 X PS- 




4. 

FLOW, kg/sec * 

0.4536 X LB/SEC 
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TABLE 3-7 (Continued) 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

temperature 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 

CD 

II 

°o 


REGULAR 

PRINTOUTS 






OATS ■ 

27 SEP 72 TIME ■ 1 7 i 2? 1 52 

TI ME ( SEC • ) 

« 5,49977-01 


5860 « 

COUNT (NO. OF CALCULATIONS! 

N 

TM(N> 

TGINi 

RG(N) 

PGiNi 

WC(N) 


deg.r 

OEG.R 

Lbm/cf 

PS! A 

L8M/SSC 

-o 

4,19298*02 

6.19048*02 

3.18945-01 

8.88030*02 

-2,20096-03 

1 

5.96109*02 

5.95668*02 

5.38487-01 

8.88167*02 

-2,18364-63 

2 

5.56269*02 

5.55796*02 

5.76865-01 

8.88416*02 

-2,14662-63 

3 

5, 14112*02 

5.13624*02 

6, 22445-01 

8.88636*02 

-2,10679-63 

4 

4.7iS24*02 

4.71322*02 

6.77177-01 

8.88826*02 

-2,06359-63 

5 

4.29544*02 

4.29027*02 

7,4204(1-01 

8,8899Q*02 

-2.01643-63 

A 

3 . 87?8l*02 

3.86748*02 

8,18357-01 

8.89128*02 

-1,96463-03 

7 

3.45033*02 

3.44483*02 

9,i4037-01 

8,89242*02 

-1,90706-63 

• A 

3.02875*02 

3.02551*02 

1,03463*00 

8.89346*02 

-1,84231-63 

9 

2.60564*02 

2.60230*02 

1.19522*00 

8.89469*02 

-1.77305-63 

10 

2.18398*02 

2.16047*02 

1.41493*00 

8.89563*02 

-1.69191-03 

11 

1.76504*02 

1.76154*02 

1.93338*00 

8,89630*02 

-1,59401-63 

12 

1.37135*02 

1.36916*02 

2,19997*00 

8,89675*02 

-1,47243-63 

13 

1.12043*02 

1.12000*02 

2,65770*00 

8,89690*02 

-1 ,. 40045-03 

6 = 270° 



DATE 

i 27 SEP 72 TIME * i 7 i 2i»t 12 

TIMEtSEC. > • 5.62500-fll 

5946 

« COUNT (NO, OF 

CALCULATIONS ) 

N 

TM( N ) 

TGtNJ 

RG(N) 

PGIN) 

WC ( VI) 


DEC . R 

OEG.R 

LBM/CF 

PSI4 

LBM/SEC 

-o 

6.19164*02 

6.19062*02 

5,18130-01 

8.86610*02 

-1.35500-63 

1 

5.95912*02 

5.95706*02 

5.37579-01 

8.86661*02 

-1.36314-03 

2 

5.56060*02 

5.55829*02 

3,75827-01 

8.8*815*02 

-1.38052-63 

3 

5.13898*02 

5.13654*02 

6.21*63-01 

8.86938*02 

-1,39923-63 

4 

4.71604*02 

4.71345*02 

6,75844-01 

8.87050*02 

-1 .41952-03 

3 

4.29320*02 

4.29044*02 

7,40542-01 

8,87l5l*02 

-1 , 44167-63 

6 

3.87051*02 

3.86757*02 

8.16682-01 

8.87241*02 

-1.46601*63 

7 

3.44793*02 

3.44478*02 

9.12175-01 

8.87321*02 

-1,49306-63 

a 

3.02607*02 

3.02410*02 

1,03293*00 

8.67398*02 

-1.52350-63 

9 

2.60779*02 

2.60063*02 

1.19343*00 

8.87496*02 

-1.55606-03 

10 

2.18077*02 

2.17841*02 

1.41328*00 

8.87576*02 

-1.59421-63 

11 

1. 761 48*02 

1.75902*02 

1.73221**00 

8.87643*02 

-1 .64028-63 

12 

1.36843*02 

1.36664*02 

2.19966*00 

8.87693*02 

-i. 69754-03 

13 

1.12006*02 

1.12000*02 

2.65253*00 

8.877i2*02 

-1 . 73138-63 




NOTES: 1. 

TEMP ,°K = °R/1.8 





2. 

DENSITY, kg/m 3 = 

16.02 X LB/FT 3 




3. 

PRESSURE, N/m 2 - 

6.895 X 1 0 3 X PS 1 




4. 

FLOW, kg/sec - 0.4536 X LB/SEC 
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TABLE 3-7 (Continued) 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO, 

TEMPERATURE 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 

0 = 300° 


REGULAR PRINTOUTS 






DATE 

a 27 SEP 72 

TIME c 17 i 2? 1 27 

TIME(SEC 

I.) » S, 74999' 

-01 

6063 

» COUNT (NO. OF CALCULATIONS) 

N 

TM(N) 

TGI N) 

R6(N) 

PGIN) 

WG( SI) 


□EC. R 

OEG.R 

L0H/CF 

PSIA 

L0M/SEC 

-o 

4.19162*02 

6.19229*02 

5.24543-01 

8.98179*02 

-2.00107-04 

1 

5.95581*02 

5.96017*02 

5,44044-01 

8.98184*02 

-2.24369-04 

2 

5,56019*02 

5.56132*02 

5.82654-01 

8.98197*0? 

-2,76223*04 

3 

5.13851*02 

5.13946*02 

6,28533-01 

8.98212*02 

-3,32020-04 

4 

4.71550*02 

4.71626*02 

6,83626-01 

6.98229*02 

-3,92322-04 

ft 

4.29257*02 

4.29313*02 

7,48963-01 

8.98248*0? 

-4,58576-04 

6 

3 • 86979*02 

3.87012*02 

8,25828-01 

6.98268*02 

-5.31126*04 

7 

3.44710*02 

3.44716*02 

9.22252-01 

6.98290*02 

-6,11750-04 

ft 

5.02500*02 

3.02480*02 

1.04471*00 

8.96314*0? 

-7,02505-04 

9 

2. 6oi5i*02 

2.60111*02 

1.20696*00 

8.98351*0? 

-7.99583-04 

10 

2.17917*02 

2.17851*02 

1.42914*00 

6,98385*02 

-9.13294-64 

11 

1.75952*02 

1.75865*02 

1.75172*00 

6.98418*02 

-1, 05057-63 

12 

1.36644*02 

1.36549*02 

2.22531*00 

6,98446*02 

-1. 22117-63 

13 

1.12001*02 

1.12000*02 

2.68059*00, 

8.98459*0? 

-1,32185-63 

e = 330° 



DATE ' 

i 27 SEP 72 TIME » I7i2?l39 

T I ME ( SEC 

.) 8 5.87485- 

•01 

6161 ■ 

COUNT (NO. of 

CALCULATIONS) 

M 

TM(SJ) 

TGI N > 

Rfi ( N ) 

PGIN) 

WC( VI ) 


OEG.R 

OEG.R 

LBM/CF 

PS I A 

LBM/S5C 

-o 

6.19422*02 

6.20000*02 

5.37013-01 

9.21362*02 

1.06847-03 

1 

5.95971*02 

5.96219*02 

3.57406-01 

9.21313*0? 

1 .02641-63 

2 

5,56169*02 

5.56639*02 

3.96539-01 

9.21230*02 

9,36590-04 

3 

5.13994*02 

5.14489*02 

6.43354-01 

9.21164*02 

8.39951-04 

4 

4,71488*02 

4.72193*02 

6.99497-01 

9.21114*0? 

7,35199-04 

5 

4.29384*02 

4.299q6*02 

7,66163-01 

9.21077*0? 

6 , 2o87l-64 

6 

3.87q96*02 

3.87635*02 

8,44430-01 

9.21053*02 

4.95363-64 

7 

3.44814*02 

3.45384*02 

9,42602-01 

9.21038*0? 

3,55955-04 

ft 

3.02387*02 

3.02954*02 

1.06804*00 

9.21030*02 

1 , 96983-04 

9 

2.60?12*02 

2.60725*02 

1.23265*00 

9.21027*02 

3.12932-65 

10 

2.17954*02 

2.18183*02 

1.46033*00 

9.21 028*0? 

-1,65193-04 

11 

1.75963*02 

1 . 76iq8*02 

1,78940*00 

9.21033*02 

-4,02253-64 

12 

1.36581*02 

1.36594*02 

2.27432*00 

9.21044*02 

-6,96765-04 

13 

1.12000*02 

1.12000*02 

2,73928*00 

9.21050*02 

-8,70*46-04 




NOTES: 1. 

TEMP ,°K - °R/1 . 8 



2. DENSITY, kg/m 3 - 16.02 X LB/FT 3 

3. PRESSURE, N/m 2 - 6.895 X 10 3 X PS I 


4. FLOW, kg/sec => 0.4536 X LB/SEC 
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TABLE 3-7 (Continued) 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

OENS ITY 

PRESSURE 

RATE 

0 = 360° 

TIMECSEC. > 

■ 3,99988-01 

regular 

PRINtOUTS 

DATE 
6296 ■ 

• 27 SEP 72 TIME « 3 7 i 2? 1 S7 
COUNT (no. OF CALCULATIONS) 

N 

tmcm) 

TG(N) 

RG(N) 

PGINj 

WG ( N ) 


DEG, R 

DEG. R 

LBH/CF 

PSIA 

LSM/SEC 

-o 

6, 19682*02 

6.20000*02 

3.33438-01 

9.56332*02 

2.17906-63 

1 

3.96142*02 

3.96469*02 

3.74173-01 

9.36406*02 

2.13272-63 

2 

3, 36477*02 

3.57076*02 

6.14139-01 

9.30184*02 

2.03380-63 

3 

3.14310*02 

5.14947*02 

6,62160-01 

9,49995*02 

1.92739-83 

4 

4.21999*02 

4.72631*02 

7.19712-01 

9.49840*02 

i . 81208-63 

5 

4.29690*02 

4.30358*02 

7.88124-01 

9.49714*02 

1,68623-83 

6 

3.87394*02 

3.88077*02 

8,68387-01 

9,49613*02 

i, 54817-63 

7 

3.43106*02 

3.43806*02 

9,69044-01 

9,49340*02 

1,39486-63 

8 

3.02894*02 

3.03309*02 

1,09783*00 

9,49480*02 

1.22231-83 

9 

2.60504*02 

2.60913*02 

1,26720*00 

9,49417*02 

1.03797-83 

10 

2.18940*02 

2.18670*02 

1.49863*00 

9.49379*02 

8.22394-64 

11 

1.76232*02 

1.76687*02 

1,83366*00 

9.49358*02 

5,62784-64 

12 

.1,36724*02 

1.37244*02 

2.32316*00 

9.49350*02 

2.41133-04 

13 

1.12002*02 

1.12136*02 

2.80929*00 

. 9,49349*02 

5.10956-35 


NOTES: 1. 

TEMP ,°K - °R/1 . 8 


2. 

DENSITY, kg/m 3 - 16.02 

X LB/FT 3 

3. 

PRESSURE, N/m 2 = 6.895 

X 10 3 X PSI 

4. 

FLOW, kg/sec - 0.4536 

X LB/SEC 



AIRESEARCH MANUFACTURING COMPANY 

Los Angeles, California 


72-8846 
Page 3-25 






03 

u 

Q- 


TD 

C 


o 

o 


03 


L. 

o 

03 

l_ 
0> 
c 
— 03 
l- cn 
oj a) 
> a: 

a) “o 

L. . — 

3 o 

4-1 (_J 

03 

i— c 
<D cn 
Q. — 
E i/i 

03 0) 

h- Q 


i 

ro 


03 

L_ 

3 

CD 


AIRESEARCH MANUFACTURING COMPANY 


72-8846 
Page 3 ~ 26 






. 002 




2 

Figure 3 - 8 indicates a maximum regenerator pressure drop of 1.26 x 10 N/m 
(1.83 ps i ) . This pressure drop, coupled with the low AP of the remainder of the 
components in the cold gas flow path, will limit leakage rates past the displa- 
cer to very low values. The regenerator represents the major portion of the 
cold gas pressure drop. However, attempts to decrease pressure drops would 
lead to losses in thermal performance which would impose a greater overall 
penalty on the VM. Figure 3~3 presents the gas flow rate at the cold end of 
the regenerator. This plot, along with the gas temperature, is used to 
determine the regenerator losses. 

An ideal regenerator will always supply gas to the cold end of the VM at 
a constant temperature' which is equal to that of the cold expansion volume 
temperature. In the actual case, this temperature varies as a result of the 
finite matrix heat capacity and heat transfer coefficient between the gas and 
matrix. This is evidenced by the temperature plots of Figure 3 - 7 . The losses 
associated with this departure from ideal conditions are estimated by inte- 
grating the excess fluid energy supplied to the cold end over a complete cycle, 
or crankshaft revolution. The loss per cycle is expressed as: 

2tt 

°loss =/ i cp (T - T ref > d9 ( 3 ''> 

Where 

Q, = regenerator loss per revolution 
loss 3 * 

d) = fluid flow rate 

cp = specific heat (function of temperature and pressure) 

T = fluid temperature 
T j. = cold expansion volume temperature 

Figures 3"7 and 3~9. along-with the specific heat characteristics of 
helium, allow a stepwise evaluation of Equation 3"1. The integrated loss as a 
function of crank angle is presented in Figure 3 “ 1 0 . The total loss per cycle 
is 0.1402J (0.0001329 Btu) . At a rotational speed of **00 rpm, this loss trans- 
lates to a total of 0.929 watts. 

Hot Regenerator Design 

1. Introduction 


The hot regenerator serves as a thermal isolator between the sump and the 
high temperature portions of the VM refrigerator. As the gas flows from the 
sump towards the hot end, it is heated to a temperature approaching the maximum 
in the cycle. At the time of flow reversal, an expansion process begins which 
allows heat to be added to the gas at essentially constant temperature. The 
expanding gas is cooled as it flows through the regenerator, thus storing energy 
in the regenerator matrix and reestablishing the temperature profile for the 
next half. Desirable design features are similar to those of cold regenerator, 
with mflior emnhasis nlaced on different items. 
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The hot displacer is considerably larger than the cold, and thus the 
pressure differential across the displacer strongly influences motor drive 
power. Thus regenerator pressure loss becomes an item of prime importance. 
Thermal efficiency assumes a less important role in the hot regenerator than 
in the cold. The thermal losses at the hot end are compensated for by additiona 
heat input on a one to one basis, while at the cold end the thermal power input 
to offset the losses is magnified by the coefficient of performance of the 
refrigerator. Void volume at the higher temperature does not affect performance 
nearly as strongly as at colder temperatures, since the mass of gas stored in 
the high temperature region is smaller, due to the high temperature. 

2 . Method of Analysis 

The hot regenerator is analyzed in the same manner as the cold. A 
description of the analytical methods used in the computer program, along with 
the physical characteristics of candidate packings, is presented in Appendix A 
of this report. 

3 . Initial Iterations 

The experience gained in the design iterations on the cold regenerator 
and from the analysis performed on the GSFC 5 watt VM allowed rapid convergence 
upon the configuration selected for the preliminary design hot regenerator. 

It was originally planned to utilize a hot regenerator packing composed at 
least partially of spherical shot. However, ratioing of the values obtained 
for the cold regenerator indicated that the pressure drop would be too high, 
as would the conduction losses from the hot to the sump region. Order of 
magnitude pressure drops were obtained by assuming that pressure drop varied 
directly with the square of mass velocity (flow per unit free flow area) and 
inversely with fluid density. The estimates of regenerator conduction losses 
were obtained by selecting a reasonable regenerator void volume, applying the 
ground rule of a 30 shot diameter annular height, and thus determining a 
regenerator length. With this type of estimate, it was determined that the 
sum of the thermal conduction losses of the regenerator walls and packing, and 
the idealized heat input was approximately equal to the allowable heat input. 
Since other factors also contribute to the heat input, this was not an accept- 
able solution. In addition, the estimated maximum pressure drops for these 

regenerators ranged from 1.103 x 10^* (1.6) to 4.55 x lO^N/m^ (6.6 ps i ) , which 
would have resulted in excessive drive motor power requirements. 

The results described above led to the decision to utilize screen as the 
hot regenerator packing material. The experience gained from the design of 
the screen packed hot regenerator in the GSFC 5 watt VM refrigerator was applied 
to selection of the core geometry of the preliminary design hot regenerator, 
described in the next section of this report. 

4 . Selected Design 

The preliminary design hot regenerator is annular in shape with an annulus 
height of 0.0038lm (0.150 in.). The total length is 0.1017m (4 in.), which 
is packed entirely with 100 mesh stainless steel screen. The screen has a 
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2 3 2 3 

porosity of 72.5 percent, an area to volume ratio of 10 , 080 m /m (256 in. /in. ) 
and a hydraulic diameter of 0.000288m (0.01132 in.). The frontal area of the 

2 2 

regenerator is 0. 000741m (1.15 in. ). 

The screen packing was selected instead of spherical shot, or a combination 
of the two, in order to minimize pressure drop. The high operating temperature 
of the hot regenerator minimizes the effect of the increased void volume of 
screen packing as compared to spheres. Pressure drop through the screen matrix 
is also predicted more accurately than for the spheres. This is caused by the 
inherent uniformity of the screen material rather than the analytical methods 
used. Slight variations in shot size or departure from a perfectly spherical 
shape may cause deviations in pressure drop. This is an important factor in 
the hot regenerator packing selection because of the strong influence of pres- 
sure drop on drive motor power. Finally, the selection of stainless steel was 
based on the superior heat sink capacity in the operating temperature range of 
interest (see Appendix A). 

Table 3“8 presents the detailed output from the regenerator analysis com- 
puter program for the preliminary design hot regenerator. The output parameters 
listed in the table are metal temperature, gas temperature, density, pressure, 
and flow rate. These parameters are printed as functions of time (angular 
crankshaft displacement) and of position in the regenerator matrix. The angular 
position is referenced to top dead center of the cold displacer. Node 0 repre- 
sents the gas inlet face at the sump end, and Node 11 the outlet at the hot 
end. Nodes 1 through 10 are internal to the matrix. Positive flow rate 
indicates flow from the sump toward the hot end. 

Figures 3“ 1 1 through 3“ 1 3 present plots of key parameters from the data 
of Table 3"8. The hot gas and matrix temperatures are plotted as functions 
of crankshaft position in Figure 3 _ 1 1 • The small difference between the two 
temperatures is indicative of adequate heat transfer. The matrix temperature 
swing of 2.6 R shows that the matrix heat capacity is sufficient. 

2 

Figure 3“12 indicates a maximum pressure drop of 4A80 N/m (0.65 ps i ) . 

This low value of pressure drop was intentionally designed into the regenerator, 
in order to assure that the required drive motor power is low. 

Any time gas is introduced into the hot end of the machine at a temperature 
less than the operating level, additional heat must be added to make up for the 
deficiency. Thus the regenerator inefficiency represents a loss that must be 
compensated for. This loss may be evaluated by use of Figures 3 - 1 1 and 3 _ 1 3 , 
which presents the gas flow rate at the hot end of the regenerator. Equation 
3 - 1 is used to evaluate the regenerator loss, as was done for the cold regene- 
rator. The cumulative losses are shown as a function of crank angle in Figure 
3-14. The total loss per crankshaft revolution is O. 85 IJ (0. 0008059 Btu) , 
which translates to a total loss of 5-67 watts at a speed of 400 rpm. 
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TABLE 3-8 


PRELIMINARY DESIGN 
HOT REGENERATOR 
PERFORMANCE CHARACTERISTICS 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 

9 = 0° 

T I SEC • ) 

= 4,4999‘*-01 

REGULAR 

PRINTOUTS 

DATE 

2301 

i 19 SEP 72 
• COUNT (no. 

T I MF * 16142:11 
OF CALCUi. AT l ONS 

M 

TM(‘i) 

TG(N) 

RG( N) 

PG(N) 

WG(N) 


OFG .R 

OEG.R 

LBM/CF 

PS! 4 

LBM/SPC 

-o 

6.21?7*5 + r2 

6 • 20000*02 

5.53472-01 

V .50556+02 

7.22990-03 

t 

A. 73^lA+C2 

6 . 71583+02 

5.12103-01 

9.50537+02 

7,11036-03 

2 

7. 57747+* 02 

7.54938.02 

4,56951-01 

9.50495+0? 

6.89636-03 

3 

R,4 d 142+^2 

8.43195+02 

4,09197-01 

9,50448+02 

4,70414-03 

4 

9.39fioi*f:2 

9.38748+02 

3,70202-01 

9.50396+02 

6.92987-03 

5 

1.031 l7+*r3 

1.02845+03 

3.37072-01 

9.50338+02 

4.37ri79-n3 

6 

1.12771+^3 

1.12012+03 

3,08033-01 

9.5n274+02 

6,22475-03 

7 

1.21422+03 

1.21174+03 

2. 83085-01 

9.50203+02 

6 , 08959-03 

A 

1.30^64+03 

1.30328+03 

2.62228-01 ‘ 

9.50126+02 

5.96319-03 

9 

1 .39A60 + r*3 

1.39435+03 

2.45518-01 

9.50043+02 

5.84337-03 

10 

1. 48363+03 

1.48162*03 

2,33273-01 

9.49953+02 

5.72788-03 

11 

1.53377 + 0-3 

1.53257+03 

2.27839-01 

9.49907+02 

5.47095-03 


9 = 30° 
TlMEtSEC. 

) * 4.6247^-Tl 

REGULAR 

PRINTOUTS 

OATS 

2384 

s 19 SEP 72 
< COUNT (NO. 

TIME a 16 142 120 
OF CALCULATIONS 

M 

T‘M‘:) 

TG ( N ) 

RO(N) 

PGiNl 

WG(N) 


OEG.R 

OEG.R 

LBM/CF 

PS I A 

LSM/SEC 


*.20V4n. + r2 

6.20000+02 

5.48819-01 

9.77836+02 

6.4Q922-03 

1 

6. 7249*3 + r2 

6.70605+02 

5.27095-01 

9.77820+02 

4,32067-33 

2 

7 , 56 4 + r 2 

7.53437+02 

4.70415-01 

9.77786+02 

6.16201-33 

} 

3*466>G7 + r2 

8.43615+02 

4.21421-01 

9.77747+02 

6.31951-33 

4 

t >.38r26 + r2 

9.35217+02 

3.81194-01 

9.77703+02 

5.89031-33 

5 

1.02977+r3 

1.02698+03 

3.47371-01 

9.77654+02 

5,77238-33 

6 

{ . 121 3^+n3 

1.1187Q+03 

3.17174-01 

9.77600+02 

5.66408-33 

7 

i.?lv>9S*r>3 

1-21037+03 

2. 91505-01 

9.77539+02 

5.56382-33 

R 

1 . 3 0 4 4 4 + f 3 

1.30197+03 

2.70060-01 

9.77472+02 

5.46997-33 

0 

1 . 39 J S44 + r-3 

1.39309+03 

2.52892-01 

5 .77400+02 

5.38091-33 

in 

1 , 48?5°+r 3 

1.48046+03 

2.40333-01 

9.77323+02 

5.29497-33 

8.25257-33 

u 

1 . 53 3l 5 + ” 3 

1-53187+03 

2.34730-01 

9.77282+02 


NOTES: 1. 

TEMP, °K - °R/1 . 8 

2. 

DENSITY, kg/m 3 - 16.02 X LB/FT 3 

3+ 

PRESSURE, N/m 2 - 6.895 X 10 3 X 

4. 

FLOW, kg/sec - 0.4536 X LB/SEC 
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TABLE 3-8 (Continued) 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 



REGULAR 

PRINTOUTS 



9 = 60 ° 



D*TE 

a 19 SEP 72 

TtMF * 16(421: 

T I ME ( SEC . ) 

= 4.74920-01 


2453 

■ COUNT (NO, 

OF CALCULATION 

N 

tm( ■:> 

TS(N) 

RG(N> 

PG(N) 

wG< N) 


DEG.R 

DEG.R 

LBM/CF 

PSIA 

L9M/SEC 

-n 

6.20468*02 

6 . 21000*02 

5.78603-11 

9.95252*0? 

3.641.41-13 

1 

6.7i7i4+n2 

6.70183*02 

5.36519-01 

9.95245*0? 

3.59861-13 

2 

7.55111+02 

7.32783+02 

4,79174-01 

9,95229*0? 

3.52192-03 

3 

R. 45J99 + P-.2 

8.42939+02 

4.29088-01 

9. 95211*0? 

3.45302-13 

4 

9. 3696 ?+i 2 

9.34569+02 

3.88111-01 

9 .9519q*Q? 

3.39158-13 

*5 

1 . 32 a 66 + -:3 

1.02636*03 

3.53375-01 

9.95167*0? 

3.33356-03 

6 

1 . 1 2 - 1 3 ? ♦ 3 

1.11811*03 

3, 22«53-01 

9,9^141*0? 

3.28120-13 

7 

T. 21 19 1 *, 13 

1.20980+03 

2.96847-01 

9 « 9*5112*02 

3.23270-13 

R 

i.30345+fl3 

1.30141+03 

2 , 75n5i -01 

9 .950BQ + 02 

3.18728-13 

9 

1 . 39430. d3 

1.39255+03 

2.57617-01 

9.95045*0? 

3.14414-13 

in 

1 . 481 74*03 

1.47996+03 

2. 4488n-01 

9.9*1008*02 

7.10248-03 

11 

1.93264.13 

1.53157+03 

2.39l9i-01 

9,94988+0? 

3.08192-03 



REGULAR 

PRINTOUTS 






DATS 

= 19 SEP 72 

TIMF * 161421 

t:mE(sec. ) 

s 4.26999-01 


2491 

a COUNT (NO, 

OF CALCUI. AT 1 01 

M 

ru v) 

TG(N> 

Rfi(N> 

PG(N> 

wet V) 

CD 

II 

VO 

0 

0 

DEG.R 

DEG.R 

LBN/CF 

PS I A 

LElM/SPC 

-0 

9.20371*02 

6 .19940+02 

5.10091-01 

9 .97806*02 

-2.71775-04 

1 

S. 71419. -12 

6 . 70450+02 

5.37652-01 

9 .97621 + 02 

-2.59092-04 

2 

7, 5445*., 12 

7.53053+02 

4. 80231-01 

9 .97850*02 

-2.36367-04 

3 

0.44922*12 

8.43278+02 

4.30039-01 

9 .97861 + 02 

-2.15965-14 

4 

9 . 36497+12 

9.34918+0? 

3.88988-01 

9.97667+02 

-1 . 97450-34 

* 5 

1 .02A21.il 

1. 02669 + 03 

3.94195-01 

9.97871+02 

-1 .8Q554-14 

6 

1 . 11984+13 

1.11843+03 

3.23727-01 

9.97873*02 

-1.65035-04 

7 

1 . 21190.13 

1.21010+03 

2.97587-01 

9.97876*0? 

-1 .50661-04 

f\ 

1 . 30303*03 

1.30168+03 

2.75768-01 

9.97877*02 

-1.37196-04 

9 

1.39410+13 

1.39278+03 

2.38324-01 

9.97879+02 

-1.24407-04 

10 

1.48137+13 

1.48009+03 

2.45394-01 

9 . 97BR0 + O2 

-1.12052-14 

11 

1 . 53^42 + 13 

1.53172+03 

2. 399Q3-Q1 

9.97881+02 

-1 .05952-04 


NOTES: 1. 

TEMP, °K = °R/1 

.8 

2. 

DENSITY, kg/m 3 

» 16.02 X LB/FT 3 

3. 

PRESSURE, N/m 2 

« 6.895 X 1 0 3 X PS 

4. 

FLOW, kg/sec ■ 

0.4536 X LB/SEC 
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TABLE 3-8 (Continued) 


NODE 

MATRIX 

GAS 

GAS 


GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 


PRESSURE 

RATE 



REGULAR 

PATNTOUTS 




9 = 120° 




DATE 

• 19 SEP 72 

T1MP a I6i4?t36 

t i me< sec . ) 

8 4,99983-01 



2351 

■ COUNT (NO, 

OF CALCULATIONS) 

N 

T rj l( N) 

TG(N> 

Rr.(N) 


PG<N> 

WG( N) 


D5G.R 

DEG.R ' 

‘ L8M/CF 


PSIA 

IBM/SFC 

-0 

6 . 20653+02 

6 . 22208+02 

5.70403-01 


9.84103*02 

-4,29506-0.3 

1 

6.71807+02 

6.74071+02 

-5i27693-01 


9.84112+07 

-4,23006-03 

2 

7. 85043402 

7.57440+02 

4.7n3i'-01 


9.84132*02 

-4.11365-S3 

3 

3 • 45289402 

8.47604+02 

4.22167-01 


9.64155+02 

-4,00901-03 

4 

9,368384.02 

9.39049+02 

3.82142-01 


9.84181+02 

-3.91409-03 

5 

1.02857+03 

1 . o3o64+Q3 

3,48067-01 


9.84210+0? 

-3.82741-03 

6 

i • I20l6*03 

1.12221+03 

3.18230-01 


9.84242*0? 

-3.74778-03 

7 

1 .21173403 

1.21372+03 

2.92636-01 


9.84278+0? 

-3.67400-03 

ft 

1.30329403 

1.30513+03 

2.71285-01 


9.84317+0? 

-3.60+89-03 

9 

1.39431403 

1.39595+03 

2,54246-01 


9,84359+0? 

-3.53924-03 

10 

1.48151403 

1.48253+03 

2, 41874-01 


9.84405+0? 

-3.47581-03 

11 

1.53299403 

1.53500+03 

2,36219-01 


9.64428+0? 

-3. 44431-03 



REGULAR 

PRINTOUTS 




9 =150° 




DATE 

a 19 SEP 72 

TIME a 16 i 4? 1 43 

TIMECSEC. ) 

» 5. 12418-01 



2601 

a COUNT (NO. 

OF CALCULATIONS) 

M 

TM( N) 

TC t N 1 

RG( N) 


PC ( N ) 

WC ( N ) 


DEG.R 

DEC.R 

LBM/CF 


PSIA 

LBM/SEC 

-0 

6.21459402 

6.23311*02 

5,54740-01 


9.57837*0? 

-6.61087-03 

1 

6.72930+02 

6.75554+02 

5.12937-01 


9.57853+0? 

-6,49170-03 

2 

7 , 562G9+02 

7.58960+02 

4.58015-01 


9.57890+0? 

-6.27824-03 

3 

R. 46384 + 02 

8.49030+02 

4.10514-01 


9.57930+0? 

-6.08637-03 

4 

9.37860+02 

9.40373+02 

3.71690-01 


9.57976+0? 

-5,91229-03 

5 

1.02949+03 

1.03187+03 

3.38576-01 


9.58026+0? 

-5.75330-03 

6 

1.12109+03 

1.12336+03 

3.09564-01 


9.58082+0? 

-5.60727-03 

7 

1.21763+03 

1.23479+03 

2 . 84663-01 


9.58144+0? 

-5.47205-03 

ft 

1. 30409+03 

1.30611+03 

2.63877-01 


9.58210+0? 

-5.34546-03 

9 

1 . 39500*03 

1.39677+03 

2.47282-01 


9.58282+0? 

-5.22535-03 

10 

1 . 48 1 95+03 

1.48306+03 

2.35225-01 


9 , 58339+0? 

-5.10943-03 

It 

1 . 53370+03 

1.53500+03 

2.29724-01 


9.38398+0? 

-5,05227-03 




NOTES: 

1. 

TEMP, °K = °R/1. 

.8 


2. DENSITY, kg/m 3 - 16.02 X LB/FT 3 


3. PRESSURE, N/m 2 - 6.895 X 10 3 X PS I 

4. FLOW, kg/ sec - 0.4536 X LB/SEC 



AIRtSEARCH MANUFACTURING COMPANY 

los Angeles. California 


72-8846 
Page 3~35 



TABLE 3-8 (Continued) 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 



regular 

PRINTOUTS 



6 * 180° 



OAT? 

■ 19 SBP 72 

TINE * 16 i 4? 1 53 

TIMECSEC. J 

1 « 3,24999-01 


2686 

■ COUNT ( NO. 

(IF CALCULATIONS] 

>1 

T *4 ( N ) 

• TGtN-) • 

RGt M) 

PG(Ni 

WG ( N ) 


DEG.R 

DEG.R 

LBM/CF 

PSlA 

L8M/SFC 

-o 

6.22930*02 

6.2452-7*02 

5.37248-01 

9.28427*02 

-7.03001-03 

1 

6.74*34*02 

6.77269*02 

4,96408-01 

9.28446*02 

-6.91612-03 

2 

7.57757*02 

7.60724*02 

4,43322-01 

9.28487+02 

-6.71210-03 

3 

8.47*53+02 

8.50717*02 

3,97448-01 

9.28533*07 

-6,S2868-h3 

4 

0.39742*02 

9.41973*02 

3.59942-01 

9.28585*02 

-6,36224-03 

5 

1 . 03079.03 

1.03340*03 

3.27892-01 

9.28642*02 

-6,21025-03 

A 

1.12733*03 

1.12482*03 

2.99800-01 

9,28706*0? 

-6,b7o68-o3 

7 

1.21362.03 

1.21620*03 

2.75677-01 

9.28776*02 

-5,94l5l-f)3. 

ft 

1.30920*03 

1.30745*03 

2.55521-01 

9.28853*02 

-5,82069-03 

9 

1.39598*03 

l; 397*5*03 

2.39428-01 

9.28936*02 

-5.70617-03 

in 

1.48253*03 

1.48377+03 

2.27751-01 

9 . 29024*07 

-5.59580-03 

u 

1.53422.03 

1.53500*03 

2.2244i-01 

9.29070+07 

-5.54l42-n3 


REGULAR PRINTOUTS 


e * 2 io° 
TIME< SEC. ) 

e 9.37500-01 


OATS 

2770 

« 19 SEP 72 
a COUNT (NO. 

TIME s 16 i 43 i 02 
OF CALCULATIONS) 

N 

T'M'J) 

TGtN) 

RG(N) 

PG ( N ) 

WG ( “J ) 


DEG.R 

DEG.R 

LBM/CF 

PS I 4 

lbm/sec 

-o 

6.23545+02 

6.25378+02 

5.22978-01 

9.03483*07 

-5 . 58501-03 

1 

6 ■ 75859+02 

6.784*3*02 

4.82627-01 

9.03496*07 

-9.49444-03 

2 

7.59723+02 

7.61920*02 

4.31048-01 

9.03527*07 

-9.33221-33 

3 

8,49749*02 

8.51650*02 

3.86904-01 

9.03562*07 

-5.18636-03 

4 

9.40860*02 

9.43040+02 

3.90P8n-01 

9.03601*07 

-5.05401-03 

5 

1.03704*03 

1.03441*03 

3.l89n-oi 

9.03645+07 

-4,93315-03 

6 

1.12362+03 

1.12578*03 

2.91579-01 

9.03693+07 

-4.82219-03 

7 

1.21496+03 

1.21711*03 

2.68097-01 

9.03746+02 

-4.71954-03 

8 

1.30628*03 

1.30831*03 

2.48464-01 

9.03805*0? 

-4.62361-03 

9 

1 . 39893*03 

1.39870*03 

2.32787-01 

9.03868*07 

-4.53277-03 

10 

1 . 48313*03 

1.48422*03 

2.21404-01 

9.03935*07 

-4,44531-03 

11 

1.53457*03 

1.63500*03 

2.16237-01 

9.03970*0? 

-4 . 40226*03 


NOTES: 1. 

TEMP, °K - °R/1 .8 

18/FT 3 

2. 

DENSITY, kg/m 3 - 16.02 X 

3* 

PRESSURE, N/m 2 - 6.895 X 

10 3 X PS 

4. 

FLOW, kg/sec - 0.4536 X 

LB/SEC 
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TABLE 3-8 (Continued) 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 



REGULAR 

PRINTOUTS 



0 = 240° 



DATE 

i 19 SEP 72 

TIME » 16143110 

TI ME f SEC. ) 

s 5,49883-01 


2839 

» COUNT (NO. 

OF CALCULATIONS: 

N 

TM(N> 

TG( N) 

RG(N) 

PGINi 

WGIN) 


OEG.R 

OEG.R 

LBM/CF 

PS I A 

LBM/SPC 

-e 

6.24255+02 

6.25640+02 

5,13767-01 

8.88139+02 

-3,01700-03 

l 

6.76861+02 

6.78778+02 

4.74427-01 

8.88145+02 

-2.96216-03 

2 

7.60763+02 ’ 

7.622*3*02 

4,23728-01 

8.88159+02 

-2,86393-03 

3 

8.50233+02 

6.52136+02 

3.79949-01 

8.88174+0? 

-?, 77560-03 

A 

9.41492+02 

9.43292+02 

3.44146-01 

8.08192+0? 

-2.69545-03 

5 

1.03793+03 

1.03464+03 

3.13495-01 

8.88211+0? 

-2.62226-03 

* 

1.12437+03 

1.12599+03 

2.86607-01 

8.88232+0? 

-2.55508-03 

7 

1. 21577+03 

1.21731+03 

2,63492-01 

8.08256+0? 

-2.49296-03 

8 

i. 30705*03 

1.30849+03 

27*4165-01 

8.88282+0? 

-2.43493-03 

9 

1.39761+03 

1.39885+03 

2.28714-01 

8,08310*0? 

-2.38002-03 

10 

1.45355+03 

1.48431+03 

2.17491-01 

8.08340+0? 

-2.32720-03 

11 

1. 53467*03' 

1.53500*03 

' 2.12391-01 

0.08356+0? 

-2,3oi'2l-h3 



REGULAR 

PRINTOUTS 



e = 270° 



DATE 

s 19 SEP 72 

TIMF s 16 i 43 1 14 

T I ME ( SEC . ) 

• 5,62294-01 


2879 

■ COUNT (NO. 

OF CALCULATIONS 

M 

TM( M ) 

TG(N) 

RRCN) 

PG(Ni 

W G ( N ) 


DEG , R 

DEG , R 

LBM/CF 

PSli 

LBM/SEC 

-o 

6.24334*02 

6.20000+02 

3.17611-01 

8.07017+0? 

1.80261-04 

1 

6.77266+02 

6.79166+02 

4.7337C-01 

8.07013+02 

1,48452-04 

2 

7.60675+02 

7.62631+02 

4.22981-01 

8.87000+02 

9,14750-05 

3 

6.50643+02 

6.52644+02 

3.-79239-01 

8.07008+02 

4.02992-05 

4 

9. 4 1 a6?+n2 

9.43930+02 

3.43470-01 

6 . B70D8+02 

-6,21042-06 

5 

1.03332+03 

1.03539+03 

3.12841-01 

8.87008+02 

-4,86358-05 

6 

1.12*74+03 " 

1 . 12684+03' 

2,85986-01 

8.87008+02 

-8.75760-05 

7 

1.21414+03 

1.21824+03 

2.62909-01 

8.87008+02 

-1.23562-04 

8 

1.30741+03 

1.30950+03 

2.43623-01 

8.87009+02 

-1.57213-04 

9 

1.39793*03 

1.39995*01 

2’, 2B21rf-01 

8.87010+02 

-1.89036-04 

10 

1.48374*03 

1.48481+03 

2.17104-01 

8.87012+0? 

-2 . 19656-04 

11 

1.53474*03 

1.53500+03 

2,12061-01 

8.87013+0? 

-2.34723-04 


NOTES : 1 . TEMP , °K - °R/1 . 8 

2. DENSITY, kg/m 3 - 16.02 X LB/FT 3 

3. PRESSURE, N/m 2 - 6.895 X 10 3 X PS I 

4. FLOW, kg/sec - 0.4536 X LB/SEC 
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TABLE 3-8 (Continued) 


NODE MATRIX 

NO. TEMPERATURE 


9 = 300 ° 

TIME t SEC . > * 3.74960-01 

M TM<M) 

CEG.R 

-C 6.23531+02 

1 4,77067+02 

2 7.60439+02 

3 A. 50357+02 

4 9,41460+02 

5 1,03311+03 

6 1.12<55+o? 

7 1.21597+0? 

8 1.30725+03 

9 1.39760+03 

10 1.48362+03 

11 1.53465+03 


GAS GAS 

TEMPERATURE DENS ITY 


GAS GAS FLOW 

PRESSURE RATE 


REGULAR PRINTOUTS 

04TE « 19 SEP 72 TIME s 16143118 
2916 * COUNT (NO. OP CALCUi ATIONS) 


TO(N) 

DEG.R 

6.20000+02 

6.75660*02 

7.58704+02 

8.4B574+02 

9.39914+02 

1.03146+03 

1.12300+03 

1.21450+03 

1.30588+0? 

1.39652+03 

1.48254+0? 

1.53396+03 


Rr.(N) 

L0M/CF 

5.24427-01 
4,«227n-01 
4.30755-01 
3.86067-01 
3.49467-01 
3.18239-01 
2.90846-01 
2.67292-01 
2,47590-01 
2. 31815-01 
2.2031R-01 
2,15061-01 


PG<Ni 

PSIA 

8.99072+07 
8.99064+07 
8.99048+07 
S. 99029+0? 
8.99008+02 
8 .98986+0? 
8.98961+0? 
8 .98934+0? 
b . 98904+0? 
8 . 98872+0? 
8.98837+0? 
8. 98619+0? 


WG< N) 
L9M/SEC 

3,45505-03 
3.38276-03 
3 . 25325-03 
3.13686-03 
3 . 03l28-rj3 
2,93492-03 
2.84648-03 
?. 7647i-n3 
2.68834-03 
?. 61408-03 
7.54657-03 
?.5l236-fi3 


REGULAR PRINTOUTS 


9 = 330° 

T I ME ( SEC . > * 5.87448-01 

N T'i('l) 

CEC.S 

-n 4.22?76 + r,2 

1 6.76343+^2 

2 7.59518+02 

o 8,49443+02 

4 9.40752+02 

5 ■ 1. 03726+03 

6 1.12375+03 

7 1.21820+03 

8 1.30653+03 

9 1.39712+03 

10 1 ,483C4+o 3 

11 1.53429+03 


DATE « 19 SFP 72 T1MF s 161 4*3 1 26 
2985 » COUNT <nO. OF CALCUi AT l ONS ) 


TG ( N ) 
CEG.R 

6 . 20000+02 
6.74470+02 
7.57030*02 
8 . 46eo9+02 
9.38198+02 
1.02982+03 
1.1.2143+03 
1.21299+03 
1.30442+03 
1.39514+03 
1.48130*03 
1.53319+03 


RO.(N) 

L0M/CF 

5.37147-01 
4. 94867-01 
4.42206-01 
3.9629H-01 
3.58643-01 
3,26566-01 
2.98434-01 
2,74257-01 
2.54031-01 
2.37846-01 
2,26050-01 
2.70628-01 


NOTES : 


PG( N) 
PSIA 

9.21596+0? 
9.21582+0? 
9.21548+09 
9.21511*0? 
9.21469+0? 
9.21423+0? 
9.21 371+0? 
9.21314+0? 
9 . 2i253+0? 
9.211A6+0? 
9.21115+0? 
9.21077*0? 


WGC *J) 
LRM/S-C 

5 , 991 75-03 
8.88841-03 
5.7Q319-03 
5.53672-03 
8.38577-03 
5.24796-03 
5.12189-03 
3.00450-03 
4.89517-03 
4.79164-03 
4 , 69196-03 
4,64289-03 


1 . TEMP , °K a °R/1 . 8 

2. DENSITY, kg/m 3 = 16.02 X LB/FT 3 

3. PRESSURE, N/m 2 - 6.895 X 10 3 X PSI 
i+. FLOW, kg/sec » 0.4536 X LB/SEC 
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TABLE 3-8 (Continued) 


NODE 

MATRIX 

GAS 

GAS 

GAS 

GAS FLOW 

NO. 

TEMPERATURE 

TEMPERATURE 

DENSITY 

PRESSURE 

RATE 

9 = 360° 
T t ME C SEC . ) * 

5,99991-01 

REGULAR 

PRINTOUTS 

DATE 

3068 

» 19 SEP 72 
■ COUNT (NO. 

TIME* » 16143135 
OF CALCULATIONS) 

N 

TMt N ) 

TG(N) 

RGCN) 

PG(N) 

WGCN) 


CEG.R 

OEG.R 

LBM/CF 

PSli 

LBM/SEC 

-o 

4. 21331*02 

6.20000*02 

5.53471-01 

9.50554*07 

7.22987-03 

1 

*.75309*02 

6.73293*02 

5.10820-01 

9.50535*07 

7.11062-03 

2 

7,581.53*02 

7.55403*02 

4,56675-01 

9.50493*02 

6.89674-03 

3 

8.48017*02 

8.45090*02 

4.09247-01 

9.50447*0? 

6.70450-03 

4 

9, 3938**02 

9.36538*02 

3.70282-01 

9.50395*0? 

6,53019-03 

5 

1 . 0339**03 

1.02824*03 

3.37141-01 

5.50337*07 

, 6.37106-03 

A 

\ . 12?5?*03 

1.11993*03 

3 , 06068-01 

9.50273*07 

6.22499-03 

7 

1.21403*03 

1.21155*03 

2.33137-01 

9.50202*0? 

6.08980-03 

ft 

1 . 3054 ?*0 3 

1.30305*03 

2.62276-01 

9.50125*0? 

5.96337-03 

9 

1.39*03*33 

1.39384*03 

2.45597-01 

9.50041*0? 

5.64350-03 

10 

t. 48213*03 

1.4«0l5*03 

2.33449-01 

9,49952*0? 

5,72794-03 

11 

1. 33371*03 

1.53246*03 

2.27848-01 

9.49906*0? 

5.67101*03 


NOTES: 1. TEMP, °K - °R/1.8 

2. DENSITY, kg/m 3 « 16.02 X IB/FT 3 

3. PRESSURE, N/m 2 - 6.895 X 1 0 3 X PS I 

4. FLOW, kg/sec - 0.4536 X LB/SEC 
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FRONTAL AREA 




LOSS = J' WC p (T-T REF )d0 
FRONTAL AREA = 0.00074lm 2 (1 . 1 5 IN. 




Cold End Heat Exchanger Design 


1 . I nt roduct ion 

The cold end heat exchanger transfers the refrigeration heat load from 
the cold finger of the VM refrigerator to the working fluid. The primary 
design criteria are: 

Minimum temperature difference between working fluid and refrigeration 
load --Thi s is an extremely important consideration. Since the refrigera- 
tion temperature is fixed, the larger the temperature difference, the 
lower the working fluid temperature the refrigerator must produce with 
inherent decreased thermodynamic efficiencies. 

Low working fluid pressure drop — The heat exchanger must provide high 
thermal performance and still not impose excessive working fluid pressure 
drop on the system. The heat exchanger pressure drop subtracts directly 
from the cold expansion volume pressure variations, which produce the 
ref r i gerat ion . 

Low void or internal volume- -Void volume at low temperature significantly 
decreases the refrigeration capacity of VM refrigerators. An evaluation 
of this effect is presented elsewhere in this report. 

Flow distribution — Uniform flow within each element of the heat exchanger 
must occur. Potential problems here are twofold: nonuniform flow leads 

to (1) reduced heat exchanger conductance, and (2) fluid elements at 
different temperatures. Subsequent mixing of these elements results is 
entropy increase and reduces the thermodynamic efficiency of the refrig- 
erator . 


Heat exchanger i nterfaces- -The cold end heat exchanger interfaces with 
the cold regenerator, the cold expansion volume, and the refrigeration 
heat load. These interface requirements control the heat exchanger con- 
figuration to some extent. The annular flow passages of the heat 
exchanger interface with the cold regenerator through a perforated plate 
flow distributor and with the cold expansion volume through ports in the 
displaced volume dome. The interface with the refrigeration heat load 
is provided by means of an axial clamp and a high thermal conductivity 
(copper plating) heat path to the heat exchanger outer wall. 

2 . Des i gn I te rat i ons 

The heat exchanger was initially configured in a similar manner as that 
of the GSFC 5 watt VM. The cold finger of the machine was essentially a con- 
stant diameter for the entire length. At the interface between the cold 
regenerator and the cold heat exchanger, the flow was directed radially outward 
to the heat exchanger slots. After traversing the slots in an axial direction, 
the gas was turned radially inward, into a constant area passage, and then 
entered the cold expansion volume. 
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After several design interations, it was found that a heat exchanger 
configuration with twelve rectangular slots 0.00953m (0.375 in.) long, and 

0. 001017 (0.04) x 0.002034m (0.08 in.) in cross section provided the desired 
design characteristics. This heat exchanger required the use of a relatively 
large filler block between the regenerator inner wall and the outer pressure 
vessel wall. The cold end of the refrigerator was subsequently changed because 
of stress considerations. The heat exchanger was reconfigured, but the essen- 
tial design features remained the same. The selected design, which is the 
direct result of the cold finger reconfiguration, is discussed below. 

3. Selected Design 

The preliminary design cold heat exchanger is shown in Figure 3“ 15- Twelve 
axial slots are arranged in an annular configuration around the regenerator 
inner wall. This configuration, with the reduction in outer wall diameter at 
the cold end of the regenerator, eliminated the use of a relatively heavy 
filler block for the heat exchanger, and thus reduced the cantilevered weight 
of the cold end. 

The cyclic flow enters (and exits) the heat exchanger from the cold 
expansion volume via ports at the outer diameter of the displacer bore. The 
working fluid traverses the axial slots and enters the radial distribution 
channels. The flow from the channels is distributed uniformly over the cold 
face of the regenerator by the perforated plate flow distributor. 

The twelve flow channels are 0.00953m (0.375 in.) long, and are rectangular 
in shape. The cross section of each channel is 0.001017 (0 .04) x 0.002034m 
(0.08 in.), with the shorter side in the radial direction. This provides a 
minimum radial height and a maximum prime surface for hea.t transfer. 

The average flow during the cyclic flow process is used for heat transfer 
performance evaluation. The total heat transfer conductance (T]hA) of the cold 
heat exchanger is 0.403 watts/°K (0.224 watts/°R) . At a heat load of 0.25 watts, 
the corresponding film temgerature o drop from the metal heat exchanger wall to 
the working fluid is 0.621°K (1.12 R) . The pressure drop, evaluated at the 

2 

maximum flow rate of the cycle, is 94.5 N/m (0.0137 psi). 

The refrigeration heat load is mounted on a removable copper end cap. 

This end cap is clamped to the cold end of the machine with spring loaded axial 
bolts, and indium foil is employed at the interface. The outer diameter of 
the cold end of the refrigerator is copper plated, in order to allow uniform 
heat distribution over the outer surface of the heat exchanger. Thus a minimum 
temperature drop from the refrigeration load to the heat exchanger wall is 
assured. 

Hot End Heat Exchanger Design 

1 . I nt roduct i on 


The hot-end heat exchanger transfers thermal energy to the working fluid 
of the VM refrigerator, thus providing the power input necessary to drive the 
system. Critical design criteria are similar to those of the cold heat exchanger, 
with emphasis on various parameters changed. 
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Lew working fluid pressure drop — As with the hot regenerator, low pressure 
drop is of prime importance in the hot end heat exchanger. This parameter 
directly affects drive motor power, and thus must be minimized. In 
addition, pressure drop reduces the cold end pressure fluctuations, and 
thus refrigeration capacity. The high fluid temperature, and thus low 
density, in the hot end heat exchanger requires that careful attention be 
paid to pressure drop. High fluid velocities occur for flow rates and 
cross-sectional areas similar to lower temperature portions of the machine. 

Low F i 1m Temperature Drop --The maximum temperature of the heat exchanger 
wall is set by structural considerations. Thus a low film temperature 
drop allows a close working fluid approach to this maximum temperature. 

This in turn increases thermodynamic efficiency, which is a direct function 
of hot end gas temperature. 

Flow distribution — Non-uniform flow in the hot end heat exchanger must be 
avoided for ;fhe same reasons given for the cold end heat exchanger. 

Low void or internal volume --The hot end of the VM refrigerator is the 
least sensitive portion of the machine to void volume. This dependency, 
caused by the high operating temperature, is discussed elsewhere in this 
report . 

Heat exchanger i nterfaces- -The hot heat exchanger provides wprking fluid 
transitions with both the hot displaced volume and the hot regenerator. 

The electrical heater, which provides the energy input to the system is 
bonded directly to the outer surface of the heat exchangei — pressure dome. 
The heat input is transferred directly to the working fluid from the 
heat exchanger walls. 

2 . Design Configuration 

The hot end heat exchanger was patterned after that used on the GSFC 5 watt 
VM , with modifications that simplify the required manufacturing processes. The 
preliminary design configuration is shown in Figure 3“16. 

The working fluid enters and exits the heat exchanger at the hot dis- 
placed volume interface by termination of the heat exchanger inner wall at a 
diameter corresponding to the hot end bearing support. The flow enters a total 
of 64 square passages conta i ned between the heat exchanger inner wall and the 
outer pressure dome. These passages direct the flow initially in a radially 
outward direction. The passages follow the pressure dome contour, until the 
flow emerges in an axial direction into transition slots that form the inter- 
face with the hot regenerator. 

The flow passages extend the full length of the heat exchanger, with no 
changes in cross-sectional area or number of slots. This configuration pro- 
vides adequate heat transfer and a very low pressure drop, which are the desired 
characteristics of the heat exchanger. The constant cross-sectional area slots 
are simpler to produce in manufacture than intersecting ribs and variable 
flow cross-sectional area. 
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Figure 3 - 1 6 . Schematic of Hot End Heat Exchanger 


Heat transfer to the working fluid is accomplished at. both primary and 
secondary heat transfer surfaces. The primary surface consists of the portions 
of the outer pressure dome which form the boundaries of the slots. Heat is 
transferred directly to the working fluid from this surface. The secondary 
heat transfer surface is composed of the sides and inner walls of the flow 
passages. Heat is conducted from the other surface to these secondary surfaces 
and then to the working fluid. The secondary surface is approximately 93 per- 
cent as effective as the primary surface. 

The heat transfer performance calculation is based on the average value 
of working fluid flow rate during the cycle, which is considered a conservative 
estimate. The flow passage length is less than kO hydraulic diameters, which 
is generally too short for fully developed laminar flow to occur. However, 
under the further conservative assumption that fully developed laminar flow 
does occur, the total T]hA of the hot heat exchanger is 6.51 watts/ K (3.62 
watts/°R) . This translates to a film temperature drop of 12.28°K (22.1°R) at 
a total heat load of 80 watts. 

The pressure drop of the hot heat exchanger is calculated at the maximum 
working fluid flow rate during the cycle, and thus represents worst case con- 

2 

ditions. The pressure drop under these conditions is 1179 N/m (0.171 ps i ) , 
which is considered an acceptable value. 

Ambient Sump Heat Exchanger 

1. I ntroduct ion 

The’ sump heat exchanger functions to transfer heat from the working fluid 
of the VM refrigerator for rejection from the system. The design criteria for 
this heat exchanger are similar to those of the cold- and hot-end heat exchange 
with changed emphasis on the various items. The preliminary design criteria 
consist of: 

• Low Working Fluid Pressure Drop --The heat exchanger must provide 
good thermal performance and yet not lead to an excessive pressure 
drop of the working fluid. As with the cold and hot end heat 
exchangers; the pressure drop subtracts from the pressure-volume 
variations in the cold expansion volume thereby reducing the 
refrigeration capacity. In addition, the pressure drop across the 
sump heat exchanger can have a significant affect on the drive motor 
power requirements unless this pressure drop is minimized. 

• Low Void or Internal Vol ume- -Voi d volumes reduce the refrigeration 
capacity of VM refrigerators by decreasing the pressure variations 
or pressure ratio; minimization of the heat exchanger internal void 
volume is therefore important. The relative importance of void 
volume in this heat exchanger as compared to other regions of the 
machine is discussed elsewhere in this report. 
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Minimization of the Film Temperature Drop — The thermodynamic 
efficiency of the refrigerator increases as the temperature of the 
gas in the sump decreases. Minimization of the sump heat exchanger 
film temperature drop allows maximum performance for the fixed heat 
rejection and heat sink temperature. 

Flow Distribution — Uniform flow within the heat exchanger is impor- 
tant for two reasons: (l) non-uniform flow leads to reduced conduc- 

tance of the heat exchanger and (2) non-uniform flow leads to fluid 
elements at different temperatures; subsequent mixing of these 
elements results in an increase in entropy and reduced thermodynamic 
efficiency of the refrigerator. 

Heat Exchanger 1 nterf aces --The sump heat exchanger must interface 
with both the hot and cold regenerators, fluid passages into the 
sump volume, and a cooling collar or clamp which provides the heat 
sink, 

2. Design Configuration 

The configuration of the ambient sump heat exchanger is shown in Figure 
3*17- This configuration is a refinement of the design evolved for the GSFC 
5 watt VM . 

The annular shaped heat exchanger is divided into two sections as shown, 
with both sections being identical in configuration except for Tength. The 
heat transfer surface of each section is formed by brazing an offset copper 
plate fin to the inside surface of the cylindrical section of the sump pressure 
vessel wall. The cylindrical sump filler block fits inside the plate fin 
thereby forming an annular passage forcing flow through the finned surface; At 
the right hand side (Figure 3”17), flow of the working fluid enters and exits 
Section 1 of the heat exchanger as it flows to and from the hot regenerator 
during the cyclic flow process. The average flow rate in this section of the 
heat exchanger is approximately 3 times that in Section 2; this accounts for 
the greater length (larger heat transfer surface) required for this section. 

At the other end of Section 1 of the heat exchanger, the flow enters and exits 
from a flow distribution passage cut into the sump filler block. This distri- 
bution passage or slot is supplied working fluid via ports that connect to the 
active cycle volumes in the crank case and behind the hot displacer as shown. 
The distribution slot is sized to provide uniform flow across the face of the 
heat exchanger. 

On the left hand side of Figure 3 “ 1 7 flow enters and exits Section 2 of 
the heat exchanger as it flows to and from the cold regenerator. This section 
of the heat exchanger is pneuma t i ca 1 1 y connected to the cold regenerator via 
the sump filler block (not shown in Figure 3“17) and the cold-end linear bear- 
ing support. The right hand end of this section of the heat exchanger inter- 
faces and shares the central flow distribution slot with the other section of 
the heat exchanger. 
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COOLING COILS OR HEAT PIPES 


FLOW DISTRIBUTION PASSAGE 


FLOW PATH 
TO COLD REGENERATOR 


SECONDARY 

FLOW DISTRIBUTION 


SECTION 2 PLATE 
FIN SURFACE 


SUMP PRESSURE 
VESSEL WALL 



777/j FLOW PATH 

TO HOT REGENERATOR 


DISTRIBUTION 
FLOW PATH 


SECTION I 
PLATE FIN . 
SURFACE 

SUMP FILLER 
BLOCK 


ALUMINUM SUMP 
COOLING COLLAR 


OFFSET PLATE FIN 
20 FINS PER INCH 


FLOW DISTRIBUTION 
SLOT 


PORT TO BACKSIDE 
OF HOT DISPLACER 

SUMP FILLER 
BLOCK 


8. EXPANDED VIEW OF PLATE. FIN SURFACE (SECTION AA ABOVE) 


Figure 3-17. Sump Heat Exchanger Conf i gurati on 
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The path for heat transfer from both sections of the heat exchanger is 
from the gas to the plate fin surface, from the finned surface through the 
pressure vessel wall and on into the aluminum sump cooling collar. Indium 
foil is placed between the sump pressure vessel and the cooling collar; this 
foil is maintained under a 100 ps i interface pressure to ensure good thermal 
contact. Heat is finally rejected from the system to water cooling coils 
brazed into channels cut in the cooling collar. This collar is designed to allow 
interfacing of the refrigerator with water cooling coils or ammonia heat pipes 
interchangeably. 


3 . Heat Exchanger Characterization 

The rate of heat transfer for each section of sump heat exchanger can be 
expressed as: 


Q - h(Ap + T\ f A f ) AT 


(3-2) 


where 


h = the average heat transfer coefficient 
Ap = basic area of the plate 
= fin effectiveness 
= fin area 

AT = average temperature difference between the working fluid and the 
heat transfer surface 


Referring to Figure 3-18 the following relations can be derived: 


Plate area 


Ap = (1 - N 6) WL (3-3) 

where N = fins per inch 
W = plate width 
L = plate length 


Fin area 



Note: This neglects the fin area exposed between fin and the sump 

filler block--a conservative approach. 
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a. NOMENCLATURE FOR OFFSET FIN 


S-41922 


SUMP PRESSURE VESSEL WALL 



b. INSTALLED PLATE-FIN SURFACE " s-73483 

Figure 3 _ l8. Rectangular Offset Plate-Fin VM Refrigerator 
Sump Heat Exchanger 
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Fin effectiveness 


Tanh (ML ) 


Tl = e 

'f ML 

e 

(3-5) 

where m = VH 

(3-6) 


k = fin material thermal conductivity and the fin length 
is given by 


L = 
e 



( 3 - 7 ) 


F low cross sectional a rea 


A f = | b - 6 (N(b - 6) + 1)| W 


Hydraul i c di ameter 


2(b - 6) (r j 

(b - 6) + |l- ej 


(3-8) 


( 3 - 9 ) 


Performance characteristics unique to a given plate fin surface are 
generally presented as plots of Colburn's j factor and Fanning's friction 
factor f as functions of Reynolds Number. Colburn's j factor is defined by: 


•i ■ & - 2/3 

P 


where C^ = gas heat capacity 
Pr = Prandtl number 


G " Af 


flow rate 


( 3 - 10 ) 


Reynolds number is defined as 


°H G 

u 


(3-11) 


Figure 3 “ 1 9 gives the Colburn j factor for the fin used in the GSFC VM 
refrigerator sump heat exchanger. This surface has 787 fins/m (20 fins/in.), 
an 0.0025A m (0.1 in.) offset length, fin length of 0.000762 m.(0.030 in.) and 
a fin thickness of 0. 0001016 m (O.OOA in.). Figure 3"20 gives the friction 
factor for this surface. The pressure drop is then computed by use of 
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FANNING FRICTION FACTOR, 


0.10 

0.08 


0.06 



Figure 3-20. Fanning Friction Factor vs Reynolds Number for 
Sump Heat Exchanger Heat Transfer Surface 
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p 


( 3 - 12 ) 


AP = 


4fL /V 


H 




where V = gas velocity 
P = gas dens i ty 


g c = gravitational constant 


4. Performance Characteristics 


The sump heat exchanger performance characteristics are summarized in 
Table 3 - 9- The heat transfer performance is based on the average flow in each 
section of the heat exchanger. For the pressure drop, maximum flows were 
assumed. The lengths of the two sections of the heat exchanger have been 
chosen to produce a heat transfer conductance ratio in proportion to the ratio 
of flows in the two sections. The film temperature drop at a heat load of 
80 watts is 2.96 K (5 • 33 R) . 


TABLE 3-9 

SUMP HEAT EXCHANGER DESIGN AND PERFORMANCE SUMMARY 


Parameter 

Section 1 

Section 2 

2 2 

Total Area, m (ft ) 

0.01475 (0.1588) 

0.00794 (0.0855) 

Maximum Flow, kg/sec (lb/sec) 

0.00329 (0.00724) 

0.001237 (0.00272) 

Average Flow, kg/sec (lb/sec) 

0.00208 (0.00458) 

0.000764 (0.00168) 

Fin Effectiveness 

0.837 

0.884 

Fluid Temperature, °K (°R) 

344 (620) 

344 (620) 

Heat Transfer Coefficient, watts/m 2 -°K 
(Btu/ft 2 °R hr) 

1489 (262) 

991 (174.2) 

Conductance (TjhA) , watts/°K (watts/°R) 

19-71 (10.95) 

7-30 (4.06) 

2 

Pressure Drop, N/m (psi) 

429 (0.0623) 

52.7 (0.00765) 


Total for Sections 1 and 2 

Conductance (hA) , watts/°K (watts/°R) 

27-01 

( 15 - 01 ) 
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Cold End Insulat i on 


1 . I ntroduction 


The cold end insulation functions to limit the heat transferred from the 
ambient atmosphere to the cryogenic cold end of the VM refrigerator. The 
insulation system analysis is somewhat unique as compared to the more common 
problem of insulating a constant temperature cryogenic heat sink from the con- 
stant temperature ambient heat source. The cold finger of the VM is composed 
of two regions: (1) a constant temperature portion comprised of the cold end 

heat exchanger and the refrigeration load mounting surface and (2) the cold 
regenerator outer wall, which has a linear temperature gradient imposed on it. 
This linear gradient varies from the cold end temperature to the sump tempera- 
ture, which is above ambient temperature. 

In discussions between NASA/GSFC, AiResearch, and Honeywell Radiation 
Center personnel, further characteristics of the cold end insulation scheme 
were agreed upon. Since the VM refrigerator is a flight prototype machine, it 
is desirable to be able to evacuate the cold end vacuum enclosure and seal it 
off without the requirement for further pumping at a later date. This require- 
ment dictates not only a leak tight enclosure but also the use of materials 
which do not outgas with time under hard vacuum. Thus the use of the cryogenic 
superinsulations (aluminized mylar) is precluded, along with any other organic 
shields supported by fiberglass pads. 

2. Preliminary Analysis 

The analysis is divided into two parts: determination of the radiation 

heat transfer from a constant temperature source (the heat shield) to a sink 
with a linear temperature gradient (the regnerator body and cold end heat 
exchanger) and determination of the heat shield temperature. The output of 
the regenerator analysis computer program (Table 3“7) indicates a linear tem- 
perature gradient along the sink. In addition, axial conduction along the 
metallic radiation shield was expected to result in a constant temperature at 
the shield. The outer wall of the vacuum enclosure is at a constant tempera- 
ture of 294°K (530°R) . In view of the regenerator temperature profile and 
considerations discussed above, the cold-end insulation analysis follows. 

The geometric view factor between two concentric cylinders is developed 
in a number of heat transfer textbooks. However, radiation between two sur- 
faces, where a temperature gradient is imposed on one, is generally not dis- 
cussed. Figure 3 - 2 1 shows the model used to determine the proper heat sink 
temperature for radiation calculations for this case. The radiation heat 
transfer is expressed as: 

. !!l^l ( T f - dx (3-13) 
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TEMPERATURE 



X, LENGTH 

S- 741 42 -A 


T x = t sink = T r MX; WHERE M = SL0PE 


Figure 3-21* Cold End Temperature Profile for 
Heat Leak Calculations 
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whe re 


Q = heat transfer rate 


= emissivity view factor area product 
L = length of cold cylinder 
^*sink = ^ oca ' temperature f(x) 


substitution of the expression for sink temperature results in 


dQ = 


° F 1-2 A 1 


1 


(T 1 - mx) 


dx 


(3-14) 


Expanding and integrating between the limits of x = 0 and x = L results in the 
following expression for the total heat transfer by radiation 


* =aF 1-2 A l 


2T 3 mL - 2T 2 m 2 L 2 + T 1 m 3 l 3 


4. 4 

m L 


(3-15) 


Noting that mL = T^ - T^, the following expression may be obtained 


F 1 - 2 A 1 


[• 8T ? - ' 


2T 3 T - .2T 2 T 2 
12 - 12 


3 4 

• 2T 1 T 2 + T 2 


The radiation heat transfer driving force is defined as 

The expression for the proper average sink temperature wi 
Figure 3“21 may be obtained from Equation 3 _ 1 6 . 


( 3 - 16 ) 


(T ) A - (T . V* 

source sink 

th the conditions of ' 


(r . . Y* - 

\ sink average } 


0.2 


4 3 2 2 3 

T 1 + T ? T 2 + T 1 T 2 + T 1 T 2 


- T, 


(3-17) 


Equation 3-16 permits determination of the radiation heat transfer from a con- 
stant temperature source to a heat sink with a linear temperature gradignt. 

It is interesting to note that if the ambient temperature of 294 K (530 R) and 
the cold end temperature of 62.3 K (112 R) are substituted into Equation 3-16 
(no radiation heat shield) the resulting heat transfer is approximately 75 per- 
cent of that which would occur if the entire regenerator wall were at the 
coldest temperature, 62.3 K (112 R) . 


Verification of the constant temperature of a single radiation shield was 
considered next. The shield is suspended in the vacuum annulus by fiberglass 
pads, and is thus not in contact with any metallic components of the VM . Ignor 
ing the conduction through the fiberglass (a non-conservative assumption) the 
shield receives heat from the outer skin, and rejects it to the regenerator 
wall, by radiation alone. Near the cold end of the machine, the shield is 
cooled by exposure to the low temperature of the cold finger. Near the sump 
end, where the outer wall and regenerator wall are at the same temperature. 
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there is no temperature potential for heat transfer to or from the shield. 

Since the shield is a good heat conductor, the colder end will cool the warm 
end, until overall equilibrium is reached. This complex problem was not solved 
by hand, but by use of a thermal analyzer computer program. 

A thirty-node network was utilized. The ten nodes along the outer wall 
were held at a constant temperature of 2Sk K (530 R) . The ten nodes on the 
regenerator wall were held at the linear temperature gradient imposed by the 
gas flowing through the regenerator. The ten nodes on the radiation heat 
shield were linked to each other by a conduction resistance resulting from a 
shield thickness of 0.010 inches manufactured from aluminum. The shield nodes 
were also linked to the inner and outer walls by radiation resistances which 
were computed using emissivity of all surfaces as 0 . 05 . 

The solution mdicated that the radiation heat shield reached a constant 
temperature of 268 K (A82 R) . The axial temperature variation was less than 
0.566 K (1 R) . The total cold end heat leak calculated for this configuration 
is 68 milliwatts. It should be recalled that the conduction through the fiber- 
glass supports to the radiation shield was neglected. This additional heat 
input to the heat shield would increase the shield temperature and thus the 
overall heat leak. Thus the 68 mw is a theoretical value that is unattainable 
i n actual . pract i ce . 

The preliminary analysis also ignored the portion of the regenerator wall 
which is hotter than ambient. In view of the heat leak calculated above and 
the shortcomings of the thermal model, it was decided to utilize an insulation 
configuration employing two heat shields, and to employ a more detailed thermal 
model . 

3 . Deta i led Ana lysis 


In addition to the preliminary results described above, an analysis per- 
formed by NASA/GSFC indicated that two radiation heat shields would be required 
to adequately insulate the cold end of the GSFC fractional watt VM. Therefore 
a relatively detailed thermal model of the desired configuration was set up 
and solved by use of the AiResearch thermal analyzer computer program. 

The 58 -node thermal network utilized as a model is shown in Figure 3 “ 22 . 
The radiation heat shields were assumed to be made of aluminum with an 
emissivity of 0.05 and a shield thickness of 0.00025^ m (0.010 in.). The 
emissivity of the outer wall, regenerator wall, and cold end was taken as 0.15- 
The areas of the various concentric cylinders were accounted for in calculating 
the emissivity view factors. 

The thermal analyzer computer program was utilized to determine the steady, 
state temperatures of the nodes. The program printout is presented in Table 
3-10 and the resulting temperature distribution is. plotted in Figure 3 - 23 - 
The output data indicates that the radiation shields attain a constant tempera- 
ture. Both shields attain equilibrium at a temperature below that of the outer 
wall and also below the temperature of the warm end of the regenerator wall. 

Thus the inner radiation shield receives heat from two sources; the heat leak 
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OUTER WALL AND BOTH SHIELDS LINKED BY CONDUCTIVE RESISTANCES THRU FIBERGLASS 

OUTER WALL, BOTH SHIELDS, COLD END, AND REGENERATOR WALL LINKED BY RADIATIVE RESISTANCES EXCEPT AT FIBERGLASS 

S-7^35 8-/ 



















TABLE 3-10 



WITH ALUMINUM HEAT SHIELDS 

STEAOY STATE SOLUTION 


NO. OF ITER. 

,= 398 

OTE'TPi .0007 

ACCEL" . 

0000 



O 

u.l 

Q 

O 

7T 

Tfmp. 

HEAT" IN.' 

RM0V 

CPIM ' ' 

KM 

"TEM B C7 

1 

530.00 

. 0000 

-1.000000 

.2300 

loo.ooono 

294.44 

2 

530.00 

.0000 

-l.OOOOOQ 

.2000 

100.00000 

294.44 

3 

530. On 

.0000 

-l.oonooo 

.2000 

103.00030 

294.44 

4 

530.00 

. oono 

-1. 000000 

.2000 

100.00000 

294.44 

5 

530.00 

. oono 

-1. ooooco 

.2000 

100. GUO 30 

294.44 

6 

530.00 

.0000 

-1.000000 

■ .2000 

lao.onono 

294.44 

7 

530.00 

.0000 

-l.oonooo 

.2000 

ino.onono 

294 . 44 

8 

530.00 

.0000 

-1.000000 

.2000 

1CJ0. 00003 

294 .44 

9 

530.00 

.0000 

-1.000000 

.2000 

100. onooo 

294 . 44 

10 

530.00 

. 0000 

-1.000000 

.2300 

loo.onona 

294.44 

11 

530.00 

. .0000 

-1.000000 

.2000 

ino.or.ono 

294.44 

12 

530.00 

.0000' 

"''-1.030000 

".2000 

loo. onooo 

* 294.44 

13 

530.00 

.0000 

-l.oacooo 

.2000 

loa.oncno 

294 .44 

14 

530.00 

.0000 

-1.000000 

.2003 

ino.onouo 

294.44 

15 

5n9 . 55 

.0000 

.010003 

.2000 

103.00300 

283.08 

16 

509.55 

• 0 j n 0 

.310000 

.2000 

loo.ooono 

283.08 

17 

509.55 

.0000 

.010000 

.2000 

ino.coooc 

283.08 

18 

509.55 

.0000 

.010000 

' .2000 

100.00000 

283,08 

19 

509.55 

.0000 

.ninooo 

. .2000 

100.00000 

2R3.08 

20 

509.56 

.oono 

.010000 

.2000 

loo.orono 

283.09 

21 

509. 56 

.oono 

.010000 

.2000 

100.00000 

283.09 

22 

509.55 

. 0000 

.010000 

.2000 

100 . 00030 

203.09 

23 

509,56 

■ . 0000 

.010000 

.2000 

10a. or, one 

283.09 

*' 24 

509.55 

.0000' 

.010003 

.2000' 

100.00000 " 

283.09 

25 

509 , 55 

.0000 

.010000 

.2000 

loo.orono 

283. -08 

26 

509.54 

.oono 

.010000 

.2000 

lnomotiOQ 

283.08 

27 

509.54 

.0000 

■ .010000 

.2000 

100.00000 

283.08 

28 

509.54 

.0000 

.010000 

.2000 

■ lob.orono 

283 . 08 

2V 

471.06 

. 0000 

.010000 

.2000 

loo.onona 

261.70 

30‘ 

471.09 

.0000 

.010000 

"".2000 

100.00003 

261.7? 

31 

471.12 

. 0000 

.010000 

.2000 

100.00003 

261 . 73 

32 

471.17 

. 0000 

.010000 

.2000 

loc.orano 

261.76 

33 

471.23 

. oono 

.010000 

.2000 

loo.ooono 

261.80 

34 

471.3? 

. 0000 

.010000 

. 2000. 

loo.onona 

261. 04 

35 

471.47 

. 0000 

.010000 

.2000 

loo.ooono 

261.93 

36 

471.66 

" .0000 

.010000 

- .2000 

” 100.00000' 

262.03 

37 

471.87 

.ocoo 

.oioono 

.2000 

IOO.OCOPO 

262.15 

38 

472.11 

. 0000 

.010000 

.2000 

100. onooo 

262.28 

39 

472.35 

• 0 0 0 0 

.nicooc 

.2000 

lno.orocQ 

262.4? 

40 

472.31 

.oono 

. 01C003 

.2000 

loo.ooono 

262.51 

41 

472.65 

. 0000 

.010000 

.2000 

ino.onono 

262.53 

42 

472.7? 

.0000 

.010003 

' .2000 

mo.crono 

26?. 6? 

43 

112.00 

. 0000 

-1.000003 

.2000 

100. onooo 

62.2? 

44 

112.00 

. 0000 

-1.000300' 

.2000 

100.000*3 

62.2? 

45 

110.40 

.0000 

- 1 . 0 GO 300 

.2000 

100.0*000 

65.78 

46 

145.Q0 

. 0000 

-1.000003 

. 2 on u 

100.0*000 

8 0.56 

47 

1«3. on 

.0000 

-1.000000 

.2000 

100.030*0 

101.67 

48 

22f>.0n 

' .ClOO 

-l.noccno 

• 2 0 0 u 

loo.ooono 

125.56 

49 

277.00 

.0003 

-l.GQOQCO 

.2000 

loo.orono 

153.89 

50 

328.00 

.0000 

-1.000000 

. 2 U 0 0 

ino.onono 

182.2? 

51 

378.00 

.0000 

-1.000000 

.2000 

ino.onono 

210.00 

52 

430. on 

.0000 

-l.nooono 

.2000 

mo.conoo 

238.89 

53 

400.00 

. oono 

-l.oonooo 

.2000 

loo.onona 

266.67 

54 

523.00 

. ocoo 

-l.nooono 

.2000 

ino.orono 

290.56 

55 

542.00 

. ocoo 

-i.crooono 

.200 0 

ino.onono 

312.2? 

56 

601.00 

.0000 

-i.oooguo 

, 2UU0 

ico. onooo 

333.69 
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Figure 3-23. Temperature Distribution, Cold Finger of Fractional Watt VM 
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from ambient and from the warm end of the cold regenerator. The heat from these 
two sources is radiated to the cold end of the machine, and to the colder por- 
tion of the regenerator wall. 


The heat leak to the cold portion of the regenerator was determined by 
two methods. The first method was summing the heat flow at the individual 
nodes; the second involved the use of the average sink temperature expression 
(Equation 3~ 1 7) - The two methods yield identical results, thus indicating that 
the number of nodes utilized in the thermal network is sufficient. 

The calculated heat leak is summarized in Table 3"11* 


TABLE 3-11 

COLD END INSULATION HEAT LEAK WITH 
ALUMINUM HEAT SHIELDS 


Portion of Cold Finger 

Heat Leak, 
mill iwatts 

Cold End 

21 .9 

Regenerator from Cold End to 

73-3 

Point Where T , , = T , . ,, 
wall shield 


Total 

95.2 


Of the total heat leak, shown in Table 3-11, 60.7 milliwatts is contributed 
by the ambient environment (outer wall) and 34.5 milliwatts is re-radiated from 
that portion of the regenerator wall which is hotter than the inner heat shield. 
The re-radiation portion of the total heat leak will exist no matter how many 
shields are used. 

The heat shields may be modified to decrease the re-radiation slightly, 
but no significant reduction is achieved unless a very complex network of radial 
baffles and shields are utilized. One possible method of decreasing the heat 
leak investigated was reduction of the shield thermal conductivity. This has 
the effect of reducing the inner shield temperature at the cold end of the 
machine, and thus the radiation to the cold end. However, the temperature of 
the shield at the warm end is also increased. This has the effect of increas- 
ing both the temperature potential to the regenerator wall and exposing more 
of the regenerator to the higher temperature. 

In order to evaluate this effect, the shield thermal conductivity was 
changed from that of aluminum to that of stainless steel. This type of shield 
construction would involve an aluminum surface deposited on a low thermal con- 
ductivity core. The total heat leak with this configuration is 93-5 milliwatts. 
The added complexity of construction is not warranted for this decrease in heat 
leak . 
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The heat shields could also be shortened at the hot end of the regenerator, 
in order to decrease the re-radiation. This concept was not investigated in 
detail. The shortened heat shields would have to be isolated from the hot por- 
tion of the regenerator by radial baffles in order to prevent their "seeing" 
the high temperature. The radial baffles would, in turn, have to be thermally 
isolated from the heat shields in order to eliminate conduction. The resulting 
network of baffles and shields would be extremely complex, and very difficult 
to manufacture. 

Another attractive method of decreasing the cold end heat leak is reduc- 
tion of the heat shield emissivity. Aluminum was utilized in this analysis, 
because of its availability, ease of manufacturing operations, and low emis- 
sivity. Other materials (gold, silver, and rhodium, for example) exhibit 
emissivity values considerably lower than that of aluminum. The primary prob- 
lem with silver and gold is oxidation of the highly polished surfaces. Rhodium 
is especially attractive from the standpoint of corrosion resistance. It has 
been used as the reflector material for carbon-arc searchlights, and may be 
readily cleaned by hand without damage to the bright surface. Thus the special 
handling required of highly polished gold and silver surfaces is not required. 

In addition, rhodium is easily plated onto base materials. 

The room temperature emissivity of rhodium is quoted as 0.02 (General 
Electric Heat Transfer data handbook, Section G51 5 • 5 » P9- 24). The thermal 
analyzer computer program was utilized to solve the thermal network of Figure 
3-22, with the emissivities of both heat shields and the outer vacuum enclosure 
equal to 0 . 02 . The resulting temperature distribution .is presented in Table 
3 - 12 .. The temperatures do not differ greatly from those of the aluminum heat 
shield configuration. However, the heat transmitted is considerably lower due 
to the lower emissivity. The calculated heat leak for this configuration is 
summarized in Table 3 “ 1 3 - 

This heat leak represents a more reasonable value than that obtained with 
aluminum heat shields. Thus the rhodium plating will yield an acceptable cold 
end insulation scheme. 

Hot End Insulation 


1 . I ntroduct ion 


The hot end insulation limits the heat transferred from the heater, the 
hot end pressure dome, and the hot regenerator outer wall to the ambient atmo- 
sphere surrounding the VM refrigerator. It is essential to limit this hot end 
loss to a minimum practical value, as heat loss in this manner is a direct loss 
of input power. 

2 . Anal ys i s 


Originally, a multiple radiation shield insulation scheme was visualized 
for the hot end. Accordingly, an analysis of the hot regenerator wall was per- 
formed similar to that for the cold end. The temperature gradient on the hot 
regenerator wall is linear, but the relative roles of the heat source and heat 
sink are reversed from the cold end. The heat sink (the outer vacuum enclosure) 
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TABLE 3-12 

THERMAL ANALYZER SOLUTION OF COLD-END INSULATION NETWORK 




WITH RHODIUM- 

PLATED HEAT 

SHIELDS 



STEADY STATE SOLUTION 





NO. OF ITER. 

s 794 

DTE'ipe .0008 

ACCEL* 

ocoo 



NODE NO" 

rFTi D_ ; 

He AT 1 T ■ 

R7KJV 

CRN 

X . T . 

TE" r C . 

1 

530.00 

.0000 

-l.OOCJQO 

.2000 

130. one. 10 

294 . 44 

2 

530.00 

.0000 

-1.303G0U 

.2000 

ino.oncno 

204.44 

3 

“STEJTTT" 

. 000 O' 

=TTT1 3 a'3DT“ 

.2033 

"inn. oldie 

— 294’. 44“ 

4 

530.00 

. 0303 

-1.GGP3.30 

.2033 

100.01073 

294.44 

5 

530.00 

.0370 

-l.nooooo 

.2300 

100.0*1000 

204 . 44 

6 

53TJ • 30 

70300“ 

=TTTI033TTG“ 

72 70 a 

'in 0; 037 73" 

294344 

7 

530.00 

.0000 

-1.000200 

• 2303 

100.00000 

294.44 

6 

530.00 

.0300 

-1.000003 

.2000 

100.03073 

294.44 

9 

5 3 O'. O' T" 

.0003 

^rrsooona rzana igq . ooooo 

20 4 . 44" 

10 • 

530.00 

. 0303 

-1.3U0300 

.2000 

ICC. OIO'-Q 

29 4 . 4.4 

u 

530.00 

.0000 

-1.030000 

.2000 

loa.oncno 

294.44 

T2 

53 0.0 0 

70700“ 

-1 • G G 0 u 0 Q 

72370 

"ICO’.' 3.7073“ 

294 T44 

13 

530.00 

. 0000 

-l.oooono 

.2000 

100.07013 

294.44 

14 

530.00 

.0300 

-1.G3000G 

.2C00 

170.01003 

294 . 44 


50 5 78'fT- 

: 70303“ 

• CTlOGfTO 

.'2000 

"100.070 70 

mi . or - 

16 

5 0 5 , 8 6 

. 0300 

.010000 

. 200C 

10Q.0O-01D 

281.05 

17 

5-05.88 

. 0000 

•aiaono 

.2310 

mo. 00013 

211 . 03 

18 

5C'5'."55" 

.7303 

.013300 

77303 

■mo.'onor.iC" 

281.03 ' 

19 

505.88 

.0300 

.010300 

.2300 

170.01073 

281.03 

20 

505 . 86 

. 0303 

.013000 

.2000 

100.30313 

211.03 

2l 

505 . '86“" 

. 0303“ 

.'0100DD“ 

72300 
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is at a constant temperature, while the temperature of the heat source (the 
regenerator outer wall) varies linearly. The temperature model utilized in the 
analysis is shown in Figure 3"24. 


TABLE 3-13 

HEAT LEAK OF FRACTIONAL WATT VM COLD FINGER WITH 
RHODIUM PLATED HEAT SHIELDS 


Heat Source 

Heat Leak, 
milliwatts 

Rerad iat i on 

17.0 

Radiation from Outer Wall 

15.0 

Conduction from Outer Wall 

15.0 

Total 

47.0 


The analysis followed that performed for the cold end, with the substitu 
tion of the source temperature into Equation 3 ~ 1 3 • The results indicate that 
the proper source temperature for use in the radiation calculations is repre- 
sented by Equation 3 ~ 1 8 - 


source average 


)* ■ ••■(’! 


L 3 2 2 3 V 

+ T 3 Tz + V2 + V 3 + T 2 


( 3 - 18 ) 


where T^ = maximum source temperature 


T 2 = sink temperature 


This analysis indicates that the heat leak is approximately only 30 per 
cent of that which would occur if the entire regenerator wall were at the 
maximum temperature. 


While this result appears attractive, the constant temperature of the heat 
shield, as evidenced by the cold end insulation analysis, is not. The inner 
heat shield would attain a temperature near that of the constant temperature of 
the hot end dome and heater. This would then result in reradiation from the 
hot shield to the cooler regenerator wall, thus introducing an additional heat 
leak. The re-radiation would manifest itself as an additional, undesirable 
internal heat leak from the hot end to the sump region. Thus, while external 
heat losses would be low, internal losses would occur. This difficulty could 
be overcome by thermally grounding the heat shields at the sump end, thus 
imposing a linear temperature gradient on the shields similar to that on the 
regenerator wall. In addition, radial baffles would be required at frequent 
intervals on the inner shield to prevent the regenerator wall from seeing 
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Figure 3"24. Hot End Temperature Profile for Radiation Heat Leak 
Cal culat ions 
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higher temperature portions of the shield. Thermal grounding of the shields 
would provide additional heat leak paths from the hot end, also increasing the 
total heat leak and thus the thermal power input required. 

An acceptable insulation configuration utilizing heat shields could be 
designed, but would be very complex and difficult to manufacture. Therefore, 
a composite insulation system utilizing fiberglass and Min-K has been utilized, 
similar to that incorporated into the GSFC 5 watt VM. This type of insulation 
system is attractive for the fractional watt VM, since the heater is bonded 
directly to the hot end dome. The heater thus operates at essentially the hot 
end temperature, as compared to that on the 5 watt machine which transfers its 
energy by radiation. The insulation system has been designed so that the tem- 
perature at the interface between the Min-K and the fiberglass is 742 K 
(1335°R) . 

This system is attractive from an assembly standpoint, and the calculated 
heat leak to ambient is 5-9 watts. This is compatible with the overall power 
input requirements for the fractional watt VM. The insulation system is shown 
on the preliminary design outline drawing. 

Pressure Drop and Void Volume Tradeoff Factors 

1. Introduction 


Many components of the VM refrigerator require design tradeoffs between 
the introduction of either void volume or working fluid pressure drop into the 
system. Void volume is defined as any volume within the machine which does 
not vary with crankshaft rotation. Since both of these factors detract from 
the refrigeration capacity of the machine, relative weighting factors for 
their effects are desirable. Unfortunately, minimization of one factor leads 
to large values of the other. Thus, if the pressure drop of a particular flow 
passage is made very low, the void volume becomes large. Conversely, if the 
void volume is minimized, pressure drop becomes relatively large. In addition, 
the effect of void volume depends on the temperature zone of the machine in 
which this volume occurs. 

The first step in the evaluation of th'e tradeoff factors was modification 
of the ideal cycle computer program. 

2 . Modifications to the Ideal Computer Program 


The VM refrigerator ideal cycle computer program, K0270, was modified to 
simplify its use. The program now has the capability to vary any of the seven 
volumes that are input over any desired range, with any desired increment. The 
charge pressure or maximum cycle pressure are now alternate inputs, and which- 
ever one of these is chosen may be varied in the same manner as the volumes. 

The maximum hot and cold end pressure drops are now also inputs. The new 
input and description is presented in Table 3"l4. 

For each parameter that may be varied, there are three inputs; the initial 
value, the increment, and the final value. The final value does not have to 
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TABLE 3-1^ 
INPUT FOR K027O 
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be an even number of increments from the initial value. The program only 
checks to see that the new current value is not greater than the final value. 

If a zero or negative increment is input, the program will change i t to a 
large positive value in order to prevent a large number of insignificantly 
small increments being calculated. Particular attention must be paid to this 
when input is charge pressure (negative values are input). As an example of 

this, it is desired to vary the charge pressure from ^.83 x 10^ (700) to 

5.51 x 10 6 N/m^ (800 psia) in 6.895 x lO^ 4 N/m^ (10 ps i ) increments. The initial 
value input would be -800, the increment is 10, and the final value is -700. 

An attempt has been made to prevent inadvertently receiving very large 
amounts of output. Therefore, all possible combinations of all variables are 
not calculated. Thus if each of the eight variables is allowed to assume ten 

g 

values, eighty cases will be calculated and not 10 . Each variable is incre- 
mented in the order of input until the final value has been exceeded. The 
variable is left at its final value and incrementing of the next variable is 
begun. it is therefore recommended that only one variable be incremented per 
run . 


If the pressure input is negative (charge pressure) the program takes its 
absolute value and calculates the mass of gas in the system from charge pres- 
sure and temperature. If the input is positive (maximum cycle pressure) the 
program determines the crank angle at which maximum pressure occurs. This 
angle is determined by differentiating the expression for pressure as a func- 
tion of crank angle andsetting the derivative equal to zero. The program 
then determines the mass of gas in the system from maximum cycle pressure and 
the expression for pressure as a function of crank angle.. Once the mass of 
gas has been determined, the cycle calculations are made in the same manner. 
The program output has been modified to include the angle at which maximum 
cycle pressure occurs. 

If the pressure drops are input as zero or blank, the calculations are 
not changed from the ideal case. If a positive pressure drop is input, the 
cycle pressures are not modified in the output. The ideal refrigeration and 
heat input are changed. In order to accomplish this, the cycle pressures are 

modified before the JVdV is evaluated. A pressure drop is evaluated at each 
crank angle increment by multiplying the maximum pressure drop by the square 
of the ratio of flow rate and maximum flow rate. If flow is into the cold 
volume, the pressure drop is subtracted from the ideal pressure. If flow is 
out of the cold volume, the pressure drop is added to the ideal pressure. The 
addition or subtraction of pressure drop as a function of flow direction is 
reversed for the hot end. If this modification of the integrals has occurred, 
a message is printed. 

With these modifications to the ideal cycle program completed, a study 
was performed to determine the effects of void volume and pressure drop on VM 
ref r i gerator performance . 
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Figure 3-26. Effect of Sump Void Volume on- Ideal 
Refrigeration Tradeoff Values 
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is more representative of the mass in a given volume. Figure 3“27 presents 
the tradeoff factors of Table 3~ 1 5 as a function of ZT. A linear relationship 
is evident. Thus the tradeoff values could be reduced to a single constant for 
this particular machine, which would then be modified by a temperature- 
compressibility factor. The use of a plot such as Figure 3-27 is considered 
more simple when trade factors at various temperatures are desired. 


TABLE 3-15 

VOID VOLUME TRADEOFF PARAMETERS FOR PRELIMINARY DESIGN 
FRACTIONAL WATT VM REFR'l GERATOR 


Location of 
Void Volume 

Temperature of 
Void Volume, °K (°R) 

Tradeoff Factor, 

AQ /A V, watts/nr (watts/cu in.) 
c 

Cold end 

62.3 (112) 

1.143 x 10 5 (1 .875) • 

Cold regenerator 

204 (366) 

0.380 x 10 5 (0.623) 

Sump 

344 (620) 

0.232 x 10 5 (0.380) 

Hot regenerator 

559 (1077.5) 

0.134 x 10 5 (0.2194) 

Hot end 

853 (1535) 

0.935 x 10 5 (0.1531) 


The effect of pressure drop between the sump and the cold end on ideal 
refrigeration was determined. This was done with the simplified approach 
described earlier under the computer program modifications. The pressure drop 
effect on ideal ref r i gerat i on is constant for maximum AP values from 6895 (1) 

to 3*79 x 10 5 N/m^ (55 psi) and the trade factor is 6.88 x 10 ^ watts/N/m^ 
(0.04738 watts/psi). The corresponding increase in hot end heat input required 
-4 2 

is 2.66 x 10 watts/N/m (1.8356 watts/psi). Again, these values are specific 
to the preliminary design VM. 

Flow Distribution Devices 


1 . I ntroduct ion 


The prime importance of obtaining uniform flow distribution through all 
heat transfer devices, both heat exchangers and regenerators, has been 
repeatedly stressed throughout this report. Test results from the AiResearch 
I R&D VM refrigerator revealed non-uniform flow in the cold end of the refrig- 
erator. The configuration of the cold end of the AiResearch VM refrigerator 
(cold end heat exchanger, displacer, and cold regenerator) is very similar to 
that of the GSFC VM refrigerator. Initial tests on the AiResearch refrigerator 
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Figure 1 - 21 . Refrigeration Tradeoff Values as a 
Function of Void Volume Temperature 
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indicated an unbalance in flow in the cold end heat exchanger and low 
temperature end of the cold regenerator. To overcome this problem, a flow 
distributor was designed and installed in this ref rigerator . As a result of 
the successful testing of this flow distributor, the same basic design was 
incorporated into the GSFC 5 watt VM ref r i gerator . A similar device will be 
employed in the fractional watt machine. This device will be optimized by 
use of the refrigeration tradeoff factors developed in the previous section 
of this report. A generalized optimization procedure is described, and then 
applied to the cold end flow distributor. The procedure will be applied to 
other flow distribution devices within the machine as detailed design 
progresses. 

2 . Optimization Procedure 

A flow distributor optimization procedure was developed utilizing the 
tradeoff factors determined in the preceding section of this report. Factors 
of prime importance affecting flow distributor design are the ratio of axial 
to circumferential pressure drop, the total pressure drop, and the void volume 
contribution to the machine. These factors are interrelated in such a manner 
that the combined effect on refrigeration is determined. This expression is 
then differentiated with respect to the axial dimension and set equal to zero. 
In this manner flow distributor dimensions are determined such that minimum 
loss of refrigeration occurs. 

In general, the dimensions of a flow distributor, such as that at the end 
of a regenerator, are fixed by the geometry of the VM, with the exception of 
the axial length. A sketch of a typical flow distributor (Figure 3 ~ 28) 
i 1 1 us t rates this. 




Figure 3 -28 • Typical Flow Distributor 
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The void volume (ignoring that of the axial holes) may be expressed as a 
function of the variable axial length X: 


V 


D 




Thus, 



(3-19) 


Assuming that the circumferential pressure loss is primarily shock 
(expansion, contraction, and turning of the gas) the pressure loss becomes a 
function of velocity head only. The velocity head varies as the square of the 
flow area. 


A r = (0 - D.) X 

C o i 


(3-20) 


The axial pressure loss of a good flow distributor is generally 10 or more 
times the circumferential pressure loss. If this multiplying factor is denoted 
by N, the total pressure loss of the flow distributor is represented by 
Equation ( 3 ~ 2 1 ) . 


A P = (N + I ) K c 



2 v ' x 2 


(3-21 ) 


whe re 

= the total loss coefficient in the circumferential direction 
w = appropriate mass flow rate 
g^ = gravitational constant 
p = gas dens i ty 

K| = is obtained by substituting Equation ( 3 ~ 20 ) and lumping all 

constant parameters together =: — — — ■ (3-22) 

2 9 C P (°o ' D i> 


■y 
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The loss in refrigeration is expressed as the sum of the pressure drop 
contribution and the void volume contribution. 


Q l - C,AP + C 2 V D 


(3-23) 


= lost refrigeration 
Cj = pressure drop trade factor 
C 2 = void volume trade factor 

substituting Equations ( 3~ 1 9) and ( 3~21 ) in Equation (3~23) » the loss in 
refrigeration due to the flow distributor is 


C K 

q l - + c 2 V 


( 3 - 24 ) 


differentiating and setting the result equal to zero, the value of X is' found 
to be : 


K 2 / 


( 3 - 25 )' 


This value of X represents the minimum loss of refrigeration due to the com- 
bined effects of pressure drop and void volume. The remainder of the design 
procedure consists of calculating the circumferential pressure drop, and 
multiplying it by N to determine the required axial pressure drop. The number 
and size of the axial holes is then selected to give this pressure drop and to 
comply with manufacturing techniques and flow distribution requirements. 

This design method assures optimum flow distributors, and also applies to 
any component for which none of the design variables is prespecified. An 
example of its nonapplicability is the flow distribution slot between the two 
sections of the sump heat exchanger. For this slot, the axial pressure drop 
is fixed as that of the heat exchanger. Once the design value of N is selected, 
the dimensions of the slot are fixed and thus the void volume. The tradeoff 
for this particular slot involves engineering judgement as to what values of N 
are acceptable as a function of void volume. 

An attempt was also made to apply this type of analysis to the hot and 
cold end heat exchangers, but the number of variables involved is too great. 

If the heat transfer conductance is considered as a design requirement, and 
thus a constant, some other variables such as hydraulic radius and ratio of 
friction factor to Colburn j modulus must also be specified. It is felt that 
this procedure would overly constrain the heat exchanger design. A single 
optimum design would not readily emerge from such an analysis. 
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3- Cold End Flow Distributor 


The optimization procedure described above will now be applied to the cold 
end flow distributor. 

The cold regenerator fixes the inner and outer diameter of the annulus 
into which the distributor is installed at 0.01078 m (0.42*4 in.) and 0.026 m 
(1.024 in.) respectively. Thus, allowing 0.000254 m (0.010 in.) thick walls 
on the distributor, D. = 0.01128 m (0.444 in.) and D q = 0.0255 m (1 .004 in.) 

(See Figure 3 “ 28 for nomenclature). Equation ( 3 ~ 1 9 ) thus yields = 0.000411 

m 2 (0.638 in. 2 ) . 

The axial to circumferential pressure drop ratio (N) is chosen as 15, and 
the total loss coefficient in the circumferential direction (K^) is taken as 

5*0. The maximum gas flow rate at the cold end of the cold regenerator (Table 
3-7) is 0.000884 kg/sec (0.001943 lb/sec); this value is reduced by a factor 
of 12, since there are 12 flow slots in the heat exchanger. Substitution of 
these values into Equation (3"22) yields K| = 0.00000234 kgf (0.00000515 lbf) . 

From the preceding section, the void volume trade factor at the cold end of the 

5 3 3 

machine (C,) is 1.143 x 10^ watts/nr (1.875 watts/in. ) and the pressure drop 
* -6 2 

tradeoff factor (C|) is 6.88 x 10 watts/N/m (0.04738 watts/psi). Substitu- 
tion of these values into Equation 3"25 yields 


X = 


2 C, ^ 


C 2 K 2 


1/3 r 


2 x 0.04738 watts x 0.00000515 lbf 
Ibf/i n . 2 0.638 in. 2 x 1.875 watts/in. 3 


1/3 


= 0.000188 m (0.0074 in.) 


Thus the optimum annular flow length for the cold flow distributor is 
slightly greater than 0.000179 m (0.007 in.). The total pressure drop for the 

2 

flow distributor is obtained from Equation 3"21 as 648 N/m (0.0940 ps i ) and 
the refrigeration loss due to this pressure drop is 4.45 milliwatts. The 

“8 3 3 

associated void volume is 7-73 x 10 °rrr (0.00472 in. J ) and the refrigeration 
loss is 8.85 mw. Thus the total refrigeration loss associated with the opti- 
mized flow distributor is approximately 13-3 mw. 

The circumferential pressure drop is obtained by dividing the total pres- 

2 

sure drop by (N + 1). Thus P = 40.5 N/m (0.00587 ps i ) and the axial holes 

C 2 

must be sized to yield a pressure drop of 607 N/m (0.0881 ps i ) . If fifty 
axial holes are chosen as a reasonable number from manufacturing cons i derat ions , 
and a loss coefficient of 1.5 velocity heads is assumed, the resulting hole 
size is 0.00034 m (0.0134 in.) diameter. 
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The flow distributor designed by this method yields minimum refrigeration 
loss. The axial holes are of a reasonable size, and the number of holes does 
not represent manufacturing problems. Other internal flow passages such as 
the sump access ports and ports from the sump heat exchanger to the ambient 
end of the hot displacer will be sized by this method as detail design of the 
VM progresses. 

Design of Cold End Seal 

The cold end seals function to control the rate of leakage which bypasses 
the cold regenerator. This function is critical in the refrigerator design 
since leakage past the seals can result in a significant loss in thermal per- 
formance. Leakage past the seals bypasses the regenerator by flowing through 
the annular space between the cold displacer and cold cylinder walls. At low 
leakage rates, the displacer and cylinder walls effectively regenerate the 
leakage fluid temperatures and the resulting thermal losses are small. As 
leakage rates increase, the walls can no longer function as an effective 
regenerator, and significant losses in overall thermal performance occur. 

1 . Design Configuration 

The basic configuration of the cold end sealing system is shown in Figure 
3"29. Seals incorporated in the cold end sealing system, from the cold end 
toward the sump end of the cold displacer, are as follows: 

(a) A 0.0254 m (1.0 in.) long, close-fit annular seal with a clearance • 
between the inner wall of the regenerator and the seal of 0.0000635 m 
(0.0025 in.) maximum. 

(b) A 5 _ groove labyrinth seal with a tip clearance of 0.0000635 m 
(0.0025 in.), a groove spacing of 0.00127 m (0.05 in.), and a 
nominal tip width of 0.000127 m (0.005 in.) 

(c) The first linear bearing, which acts as an annular seal with a 
clearance of 0.00001017 m (0.0004 in.) and a length of 0.0254 m 
(1 .0 in.) . 

(d) The bearing support member,which acts as an annular seal with a 
clearance of O.OOOO89 m (0.0035 in.) and a length of 0.0429 m 
(1.69 in.). 

(e) The second linear bearing which is identical to the first. 

The arrangement of the machine is such that these sealing elements are 
in series. 
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Figure 3 - 29 . Cold Displacer Sealing Design Configuration 

As with the GSFC 5 watt VM the linear bearings and the bearing support 
are used as part of the sealing system. The mechanical arrangement of the 
machine allows the use of these components as seals without penalty. The only 
disadvantage in their use is that the leakage rate is dependent on bearing 
clearance and increases as the bearing wears. However, the loading and rate 
of wear of these bearings is low, and will provide over two years of operation 
before bearing wear affects the performance (even in the worst case analysis). 

2 . Method of Analysis 

The correlation for analyzing the leakage past both labyrinth and annular 
seals were developed and modified to match the seal test data in Ref. 3"1 • 

The appropriate equations are: 

for Labyrinth seal 

ub, 2 T R7| 

AP. = j (3-26) 

c A p o k l C 2g c 


|cAami) airesearch manufacturing company 


72-88^6 
Page 3-83 



for Annular Seal 


A p i 



48u 

D hV 


Then noting that: 

N N 

■' P T =?- 4P i W 


Where 


(3-27) 


(3-28) 


P.p = the total pressure drop across the series of sealing elements 
i and N = the identity and number of elements respectively 
then further noting that 

di_ = (i). = d). = d). (3-29) 

T L A i 

that is, the flow past each element is the same, the pressure drop can be 
computed as a function of leakage rate from Equation (3 - 28) . 

3 • Performance Characteristics 

a . Leakage Rate 

The configuration of the cold end of the refrigerator is such that seal 
elements can be added in series (or subtracted) as desired by providing flow 
passages Into the active cycle volume at different locations. Figure 3 ~3 0 
gives the leakage rate as a function of pressure drop for the combination of 
all seal elements in series. The design combines all possible seal elements 
to minimize the losses due to leakage. 


Figure 3 “ 30 presents leakage rate data for seal systems with both new 
bearings and bearings after 2 years of wear. As the bearings wear, the 
clearance of the annular seal (the seal which the bearings form with the cold 
displacer) increases and hence the leakage rate increases. The bearing wear 
rate used in the analysis is based on the bearing material tests described in 
reference 3 - 1 • The use of these wear rate data is believed very conservative 
due to the low loading of the linear bearings compared to the wear rate test 
cond i t i ons . 


The dashed vertical line in Figure 3"30 corresponds to the preliminary 
design limit level of pressure drop across the cold end seal. This pressure 
drop level is based on the maximum cold end flow rate at a rotational speed of 
400 rpm and thus represents the worst case conditions. 


f . 1 
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b. 


Thermal Losses 


The cold end leakage is important, in that it influences refrigerator 
thermal performance. A good approximation of the thermal losses associated 
with the cold end leakage is provided by the simple model developed below. 


Figure 3~3 1 depicts the basic elements of the thermal loss model employed. 
Taking an element of length (dX) along the annular flow passage between dis- 
placer and inner regenerator wall, the following differential equation can 
be wr i tten : 


Where 



+ 


hA 

5T T f 

p 


hA 

— — j 

ouC uj 
P 


(3-30) 


T f = temperature of the leakage gas at any location X 

T = temperature of the regenerator and displacer walls at any 

location X 


h = local heat transfer coefficient between leakage gas and the 
surroundi ng walls 


A^ = heat transfer area per unit length along the leakage path 
ub = rate of leakage 



t 

x = 0 

T = T 0 
AT x = 0 

T 0 =T , 

Figure 3 — 31 . Leakage Thermal Loss Model 



w 


t 

= X 

e 

= T a 

= f(x) 
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Then if a linear temperature distribution is assumed along the displacer and 
regenerator walls and ratio of heat transfer coefficient to gas heat capacity 
is taken as a constant, Equation 3"30.can be written as: 


dT f T 

3i r + “ T f 


= or (T -T ) $• + <»T 

v a o X c 


( 3 - 31 ) 


Where the new terms are: 

hA 

a = die” 

P 

T = wall temperature at the sump end of the displacer 

d 

T o = refrigeration temperature or wall temperature at X = o 
X g = length of the displacer 

Solving Equation 3~31 with the boundary conditions of = T q at X = o 
yields: 

- T, = 4~ 2, lx- 1 + - e"" X + T (3-32) 

f X Of a ) o v ' 

e l ; 

if the leakage flow were completely regenerated by the walls, the temperature 
of the fluid would be that of the wall at the end of the displacer; that is: 

T r » T @ X - X 

f a e 

or if flow (leakage) is considered in the reverse direction 

T , @ X = 0 

f o 

In setting up the relation for T^., this latter condition was used as a boundary 

condition assuming flow in the positive X direction (see Figure 3“3l)* Relation- 
ships for flow in the reverse direction are similar but are not developed here. 

In the actual case, due to the cyclic operation of the VM ref r i gera tor , the 
leakage flow does reverse direction. The leakage losses can be estimated, 
however, by considering flow in one direction with an appropriate time span. 

With this cons i derat i on the thermal losses per cycle due to leakage can be 
expressed as : 

\ = fn C P {(Wdb dT < 3 - 33 > 
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Where 


*» 

0 implies integration around the cycle and 
t = time 

Assuming a constant leakage rate, the leakage thermal losses can be expressed 
as : 


Ql= ic p j (T a- T o ) 5r c- e ' aXe )j ^ 3 - 34 ) 

Figure 3“32 gives the losses (Equation 3“34) as a function of leakage 
rate. The thermal losses at leakage rates corresponding to the design limit 
levels of pressure drop and associated cold end leakage rates (see Figure 3 ~ 30 ) 
are shown in Figure 3~32. It is noted these losses are estimated for the 
machine after two years of wear; thermal losses due to leakage are negligible 
for a new machine. In fact, the leakage thermal losses for the machine after 
2 years of wear are expected to be considerably lower than indicated in Figure 
3-32 for two reasons. First, the bearing wear rate used in estimating the 
leakage is considered to yield a conservative result as previously mentioned. 
Second, constant leakage rates at their maximum levels were used in estimating 
the losses by use of Equation 3 ” 3 3 * In the actual case, the leakage rate is 
a periodic function with an average value on the order of 70 percent of the 
maximum. Thus not only are the losses lower due to the decreased leakage flow, 
but the regenerator and displacer walls more effectively regenerate the leakage 
gas, further reducing the losses below the maximum values indicated in Figure 
3-32. 

Design of Hot End Seal 

The hot-end seal functions to control the leakage rate of the working 
fluid which bypasses the hot regenerator. Leakage bypasses the regenerator 
by flowing through the annular space between the hot displacer and the inner 
wall of the regenerator. Excessive leakage results in a loss in thermal per- 
formance, therefore hot end seal design is an important consideration. 

At low leakage rates, the hot displacer and cylinder walls effectively 
regenerate the leakage fluid temperatures and the resulting thermal losses are 
small. As the leakage rates increase, the walls can no longer function as an 
effective regenerator; thus significant losses in overall thermal performance 
result. 

1 . Design Configuration 

The hot end seal configuration (Figure 3 ” 3 3 ) consists of the following 
two elements, from the hot end toward the sump end: (l) a 5-groove labyrinth 

seal with a tip clearance of 0.0000635 m (0.0025 in.), a groove spacing of 
0.00127 m (0.05 in.), and a nominal tip width of 0.000127 m (0.005 in.) and 
(2) a 0.025** m (1.0 in.) long close fit annular seal with clearance between 
the inner wall of the regenerator and the seal of 0.0000635 m (0.0025 in.) 
maximum. 
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Figure 3-32. Cold End Thermal Losses Due to Leakage as 
a Function of Leakage Rate 
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Figure 3 ~3 3 * Hot-End Seal Configuration 
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These sealing elements are similar to those used in the cold-end seal 
except for the diameter. One major difference is the use of linear bearings 
as part of the sealing system for the cold end. The bearings are not used 
in the hot-end seal and thus the performance of this seal is independent of 
wear (or operational time). 

It should be noted that the seal is located in the highest temperature 
region possible within the machine. The leakage rate is an inverse function 
of the temperature, therefore, this location minimizes the leakage rate for a 
given seal configuration. The seal design selected provides very low losses 
as discussed in the following paragraphs. 

2 . Method of Ana lys i s 

The method of analysis is identical to that used for the cold end seals. 

3 . Performance Characteristics 

Figure 3“34 gives the hot-end seal leakage rate as a function of the 
pressure drop across the seal. The design limit pressure drop and the asso- 
ciated leakage rate shown in Figure 3“34 corresponds to the preliminary design 
goal maximum pressure drop across the seal. This pressure drop includes the 
maximum pressure drops across: (l) sump ports to backside of hot displacer, 

(2) Section 1 of the sump heat exchanger, (3) interface between sump heat 
exchanger and hot regenerator, (4) hot regenerator and (5) the hot-end heat 
exchanger. Due to the use of maximum pressure drops--pressure drops, correspond- 
ing to the maximum flow rates at a rotational speed of 400 rpm--the indicated 
design leakage rate is conservative. 

Comparing Figure 3~34 with Figure 3 “30 (Figure 3~30 gives the cold-end 
leakage rate vs pressure drop) a major difference is noted. The hot-end seal 
is completely non-contacting and is not subject to wear. The hot-end leakage 
is therefore independent of operational life. 

Figure 3~35 gives the thermal loss for leakage past the hot-end seal as a 
function of leakage rate. The method of calculating this thermal loss, which 
takes into account the regenerative effect of the displacer and cylinder walls, 
is identical to that used for the cold-end seal thermal losses. The loss 
indicated for the design limit conditions is not significant compared to the 
80 watts maximum power input to the hot-end required to operate the system. 

HEAT PIPE INTERFACE 

This section summarizes the effort expended on Subtask 1.3 entitled Heat 
Pipe Interface Definition. The objective of this effort is to evaluate the 
design of a removable heat pipe interface assembly, connecting the VM sump to 
a space radiator, and to carry the design only to a depth necessary to define 
the interface requirements. 

The design approach is the same as that used in the GSFC 5 watt VM effort, 
and it employs an interface clamp design which is adaptable to use with either 
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Figure 3"34. 
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a water cooling system or an ambient heat pipe assembly. The latter consists 
of one or more heat pipes on each side of the refrigerator crankcase designed 
to allow for the redundancy requirement. 

Design Requirements 

The system parameters influencing the design of the heat pipe assembly 
are listed below. 


Heat rejection (including motor losses) 
Redundancy requi rement 
Sump diameter 

Sump length (allocated for heat pipe assembly) 
Maximum cycle pressure 
Minimum cycle pressure 


70 w 

50-percent 

0.0635 m (2.5-in.) 

0.0556 m 
(2.1875-in.) 

6.895 x 10 6 N/m 2 
(1000 psia) 

6.10 x 10 6 N/m 2 
(885 psia) 


Sump temperature 334°K (140°F) 

The present effort will consist of the following items 


(a) 

Sink temperature selection 

(b) 

Heat pipe 

des i gn 

(c) 

Heat pipe 

assembly design 

(d) 

1 nterf ace 

def i ni tion 

(e) 

1 nterf ace 

temperature drop 


Sink Temperature Selection 

The sink temperature in this case refers to the spacecraft radiator tem- 
perature. The lower the sink temperature the higher the available temperature 
drop between the sump of the refrigerator and the radiator, which simplifies 
the design of the heat pipe assembly. On the other hand, the lower the sink 
temperature the larger the radiator area which is required to reject a given 
heat load. 


A NASA analysis* of the heat rejection radiator was used in conducting 
system tradeoffs to select the sink temperature for the heat pipe assembly. 
The analysis is based on the following assumptions. 


*By Dennis J. Asato, NASA, Goddard Space Flight Center. 
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(1) Isothermal radiator surface 

(2) Environmental factor (5): 0 to 0.6 

(3) | = / (^ , spacecraft orientation) 

Variables in the analysis are defined below 
A = energy absorbed by radiator 
a = absorpt i vi ty 
E = energy emitted by radiator 
e = emi ss i vi ty 

(energy absorbed) 

^ - (energy emitted by isothermal radiator) 

P = power, watts 

2 

S = surface area, ft 

a . 53.2858 x I0~'° - "i” 5 , 

ft 2 °K 4 

T = temperature, °K 

The relationship between the radiator area and the radiator temperature 
is obtained as follows 


E = P + A 
E = ct e T 4 S 

Eliminating A from Equations (3-35) and (3-36) one gets: 
P = E (I - ?) 

Substituting Equation (3-37) into Equation (3-38): 

P = (I - ?) o e T 4 S 


(3-35) 

(3-36) 

(3-37) 

(3-38) 

(3-39) 


Se _ I 

P (l - g) q T 4 


(3-40) 
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The radiator surface area versus the environmental factor £ is shown in 
Figure 3“ 36 for various radiator temperatures. Figure 3"37 depicts the radiator 
area versus the radiator temperature for two values of § (0-2 ar >8 0.45). A 
decrease in the radiator temperature from 322 (120) to 316. 5°K (110°F) results 
in an increase of about 6-percent in the radiator surface area. Hence, a 
preliminary choice of sink temperature range of 316.5 (110) to 322°K (120°F) 
appears to be satisfactory. Hence, the available temperature drop for the heat 
pipe assembly ranges from 4.1 (20) to 16.7°K (30°F) . 

Heat Pipe Design 

Design problems such as working fluid selection, material selection, wick 
design, selection of a heat pipe diameter, power transport capability, tempera- 
ture drop, etc., are discussed in the following paragraphs. 

1 . Working Fluid Selection 

There are four basic criteria that govern the selection of a working fluid 
for the heat pipe. These criteria that vary in relative importance with the 
particular application are: 


n° h fq 

• Large transport factor, N = a - 

^ ^ l 

• High thermal conductivity 

• Moderate vapor pressure 

• Chemical compatibility with wick and container 

The transport factor must be chosen as large as possible for the 

operating temperature range. Figure 3 - 38 presents the transport factor as a 
function of temperature for several working fluids. Water is shown for com- 
parison purposes only; it is not a candidate fluid because of its high freezing 
point. It is clear that ammonia is the most attractive on the basis of trans- 
port factor since is from 1.3 to 10 times that of its closest competitor, 
methanol . 

It is also desirable to have a working fluid with a high thermal conduc- 
tivity, which reduces the conduction temperature drop through the wick. Ther- 
mal conductivities of candidate working fluids for the present application are 
shown in Figure 3“ 3 9 • Again, ammonia is shown to be superior. 

The vapor pressure of the working fluid should not be so high that the 
heat pipe must become a pressure vessel with prohibitively thick walls. It 
should also not be so low that the vapor densities result in high-pressure 
drops in the moving vapor. The vapor pressure- temperature curves for the two 
fluids discussed above are shown in Figure 3~40. Ammonia has the highest 
saturation pressure. Its high vapor pressure at the temperature level of 
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Figure 3-36. Space Radiator Design Chart 
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Figure 3 - 3 7 * Effect of Radiator Area On 
Radiator Temperature 
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Figure 3 _ 39- Thermal Conductivity of Candidate Fluids 
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interest will yield low vapor pressure drops with corresponding small 
saturation temperature drops between evaporator and the condenser. Stated 
another way, it will allow small vapor passage diameters without increasing 
the vapor pressure drop. The vapor pressure is also low enough to allow heat 
pipe construction utilizing reasonable material gages. The vapor pressure of 
methanol is too low for the present application. 

Finally, there must be chemical compatibility between the working fluid, 
the wick, and the container. 

2 . Material Selection 

Aluminum is a very reactive metal. Its resistance to corrosion and 
chemical attack is due to a thin film of aluminum oxide (Al^O^) that forms 

on the surface immediately upon exposure to the atmosphere. The corrosion 
resistance of aluminum to materials other than water is generally very good 
and has resulted in the use of aluminum in the handling of many substances. 

This reactive metal depends upon its alumina surface film for its resistance 
to chemical attack, and any material that degrades or removes this film may 
cause increased corrosion. 

In contrast to most corrosion problems, in the heat pipe application, the 
structural integrity of the tube wall is not the primary consideration. Rather 
it is the quantity of noncondensible gas generated in the thermal decomposition 
or corrosion reaction that must be avoided because. such gas collects in the 
condensing region, and effectively prevents the flow of working fluid vapor into 
this region. 

While the compatibility of ammonia with aluminum has been demonstrated, 
considerable care must be taken to assure that the ammonia is pure. Purities 
less than 99-99 percent have been found to be corrosive in heat pipes, particu- 
larly if the impurities contain water vapor. 

On the basis of the above discussion, it appears that ammonia-aluminum is 
the best combination for the present application. 

3 • Wick Design 

The desirable wick characteristics are dictated by the two functions 
served by the wick, i.e., to deliver and distribute liquid to the evaporator 
region and to serve as a heat transfer path between the heat pipe container 
and the liquid-vapor interface where phase change occurs. The ideal wick con- 
figuration will have : 

• Low resistance to liquid flow 

• High capillary pumping force 

• Low AT at evaporator and condenser 

• Insensitivity to nucleate boiling 
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For a homogeneous wick, properties that enhance one of the above characteristics 
will penalize another. For instance, low resistance to fluid flow means large 
pores while high capillary pumping means small pores or using a small pore wick; 
low resistance to fluid flow means a thick -wick while low AT means a thin wick. 
In the case of a homogeneous wick, therefore, selection of the wick con- 
figuration involves a tradeoff between the various properties. 

One method of obtaining wick configurations that possess all of the 
desirable characteristics to a greater extent has been the use of composite 
wicks. Several configurations have been proposed, and some of these (channels, 
multilayer screens of different mesh sizes, and I-beam wicks) have been evalu- 
ated experimentally. 

Additional factors influencing wick design are chemical compatibility with 
the working fluid and the container (in this application, aluminum), and separa- 
tion of the liquid and vapor flow streams. 

Two candidate wick design concepts are shown in Figure 3"4l • In Concept 1, 
a composite wick is made of two layers of screen of different mesh sizes. The 
loose mesh is used for the liquid return from the condenser to the evaporator 
end to reduce viscous losses in the liquid path. The fine mesh screen is used 
for capillary pumping and to separate the liquid from the high velocity counter- 
current vapor flow. A retaining spring is used to maintain the screen in con- 
tact with the tube wall in applications which require bending of the heat pipe. 

In Concept 2, the main liquid path is isolated from, the heat path, which 
usual l.y results in a smaller temperature drop in the heat pipe. Again in this 
concept, the loose mesh screen is used for the liquid path while the fine mesh 
screen provides the capillary pumping. At the evaporator end liquid is 
delivered to the tube wall through tabs on both ends, and the liquid is dis- 
tributed around the tube through the narrow circumferential grooves. 

While Concept 2 provides lower temperature drop in a heat pipe, Concept 1 
offers more structural integrity, particularly during bending of the heat pipe. 
The advantage of smaller temperature drop is of negligible importance in the 
present application due to the low level of the power transported. At these 
levels, the differences in AT between the two concepts are quite small. 

Since Concept 1 was used in the ambient heat pipe for the 5 _ watt VM 
machine, and since experimental data are available for these heat pipes, it 
was decided to use the same concept for the heat pipes on the fractional watt 
VM machine. 

A . Power Transport Capability 

The design of a heat pipe to meet the specified set of conditions starts 
with satisfying the pressure drop criterion: 


iP c 2 4P l _ + 4P v 4 4P g 


net capi 1 lary 

;> 

1 i qu i d pressure 

4- 

vapor pressure 

± 

body force 

head 


d rop 

T 

drop 


head 
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TUBE WALL 



ONE LAYER OF 30 x 30 
MESH SCREEN 

ONE LAYER OF 200 x 200 
MESH SCREEN 

RETAINING SPRING 


CONCEPT 1 - COMPOSITE SCREEN WICK 



CONCEPT 2 - I-BEAM WICK 

S-77563 


Figure 3-^1- Candidate Wick Designs 
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The heat pipe can operate without drying out the wick as long as the capillary 
pumping pressure term is greater than or equal to the pressure losses in the 
liquid or vapor, and those due to gravitational or acceleration field effects 
(which may help, hinder, or not exist). Each of these pressure terms is con- 
sidered individually in the following paragraphs. 

The capillary head for a square mesh screen has been shown to be idepen- 
dent of the contact angle, and, hence is predicted by the equation 

AP c = y- (3-42) 

c 


where a is the surface tension and r is the 

c 

layers of screen, most of the available data 
correlated by the following equation for the 
3-2) . 


effective pore radius. For single 

were shown to be successfully 
effective pore radius (e.g., Ref. 


d + 6 
r c = 2 


(3-43) 


where d is the wire diameter and 6 is the wire spacing. 

Since in the present concept (composite screen wick) the wire spacing of 
the fine mesh screen is much smaller than the coarse mesh screen it is treated 
as a single layer. 


The liquid pressure drop is predicted by the following correlation 


° \ ( UL a > 

iP t ' 2 P* K h fg a c 9c 


(3-44 ) 


Where 

Q 



and a 

c 


heat flux 

liquid vi scos i ty 

1 i qu i d dens i ty 

latent heat of vaporization 

heat pi pe length 

adiabatic length of heat pipe 

permeability of the screen matrix 

gravitational constant 

total cross sectional area of the wick 



AIRESEARCH MANUFACTURING COMPANY 


72-8846 
Page 3-105 



The vapor pressure drop is predicted by, Reference 3 - 3, 


4ttQ p, (L+L ) 

AP = _ — £l 

v ,2 

p h,. a 9 
v fg v c 

where the subscript v in this case refers 

Finally, the body force head is given 

AP = p, ^ [0. sin 0 + L cos 0] 

g i g c i 


to vapor. 

by, Reference 3 ” 3 » 


(3-45) 


(3-46) 


where 6 is the angle the heat pipe makes with the vertical and 0 ; I s the inner 
diameter of the heat pipe. 

In the present application, the ammonia vapor pressure is quite high at 
the operating temperature which results in high vapor density, and, conse- 
quently, low vapor pressure drop as compared to the liquid pressure drop. 
Hence, the vapor pressure drop AP^ can be neglected. For horizontal operation 

the second term in Equation 3~46 can be neglected and the body force head 
reduces to 


AP = p. D . 

g 1 


Hence, the pressure drop criterion becomes 


2 o 


D. = 


Q n. (L+L ) 
_ max 'L 3 


( d±6) l 9 C i !p ( K h fg a c g c 


(3-47) 


(3-48) 


The 
Equat ion 


heat transport capability of the heat pipe can be obtained from 
3-48 as fol lows : 


Q 


max 


V h fQ 


l 




if&a ) 

a 9 C / 


(3-49) 


All the parameters on the right hand side of Equation 3-49 are known_ 
once the heat pipe diameter is selected. The value of the permeability K was 
determined on the basis of experimental_data obtained on the ambient heat pipe 
for the 5 watt VM refrigerator. These K values compare quite well with those 
published for similar wick designs. 


Equation 3 - 49 was used to predict the maximum power transport capability 
for various heat pipe diameters and lengths. Two heat pipe diameters were 
selected for this parametric study (0.0127 (0.5) and (0.00952 m (0.375-in.)) 
covering the range of interest. The upper bound of the range is determined by 
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the sump length allocated for the heat pipe assembly (0.0556 m ( 2 . 1 8 75 ~ i n . ) ) , 
the desired redundancy which requires a multi-heat pipe approach, and the sump 
diameter which determines the radius of curvature of the heat pipe. The lower 
bound of the diameter range is dictated by the simplicity of fabrication and 
assembly. 

The length of the ambient heat pipe controls the location of the VM 
refrigerator in relation to the space radiator on a particular spacecraft. 

Hence, the effect of the heat pipe length on its maximum power transport 
capability was investigated. 

The length of the condenser end of the heat pipe, and hence the radiator 
length was assumed constant at 0 . 6 1 m (2-ft) . The adiabatic length of the heat 
pipe was changed to vary the overall length of the heat pipe. 

The results of this analysis are shown in Figure 3-42 for a heat pipe in 
the horizontal position. The operating temperature of the heat pipe was assumed 
to be 332°K ( 1 20°F) . 

5. Temperature Drop 

The temperature drop 
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The vapor AT^ is usually quite small and can be neglected. 

For a round heat pipe, with a wick design similar to that shown in Figure 
3 _ 4l Concept 1, the AT can be predicted by the following equation 


AT = Q 
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POWER TRANSPORT CAPABILITY, WATTS 
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Figure 3-42. Effect of Heat Pipe Length on 
Power Transport Capability 
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where 


r Q = heat pipe outside radius 
r. = heat pipe inside radius 
r.. = wick inside radius 

i i 

= evaporator length 
£ c = condenser length 
k = wall thermal conductivity 

K = effective thermal conductivity of the wick-liquid matrix 
Q = heat rate transported by the heat pipe 
In obtaining Equation 3”51» the following assumptions have been made 

(a) Heat input to the heat pipe at the evaporator end is over half the 
circumference only. 

(b) Heat rejection from the heat pipe at the condenser end is over 
the entire circumference. 

(c) The transfer process in the wick is molecular conduction, which is 
substantiated by test data on an identical heat pipe design. 

(d) Temperature drops in the vapor stream, and at the vapor-liquid 
interfaces in both the evaporator and condenser, are negligible. 

All parameters in Equation 3 “ 5 1 are usually known a priori except K. The 
value of K for the present application was determined on basis of the data 
obtained on the identical design used on the 5 - w refrigerator. A sample of 
this data is shown in Figures 3“43 and 3”44. A schematic of the test system 
used to obtain this data is shown in Figure 3 - 45- 

The heat pipe AT for the present design is shown in Figure 3 _ 46 versus 
the power transported for two heat pipe diameters 0.00952 

(0.375) and 0.0127 m (0.5 - in.) and subject to the previous assumptions. The 
increase in AT with the decrease in heat pipe diameter is due mainly to the 
decrease in the thermal resistance at the evaporator end. 

6 . Selection of Size and Number of Heat Pipes 


The 50-percent redundancy requirement necessitates the use of several 
individual heat pipes in the ambient heat pipe assembly. The assembly must be 
designed such that if 50-percent of the heat pipes fail, the heat 
pipe assembly will continue to transport a total of 70-w to the radiator while 
maintaining the sump skin temperature at a maximum value of 334°K (140°F). 
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AVERAGE LOCAL TEMPERATURE 


NOTES: 

1 • TEMP, °K - (°F+460)/1 .8 
2. DISTANCE, m = 0.0254 X IN. 
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Figure 3*43- Effect of Temperature Level on Heat Pipe Performance 
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AVERAGE LOCAL TEMPERATURE 


NOTES: 1. TEMP °K = (°F + 460) /I. 8 

2. DISTANCE, m = 0.0254 x IN. 



Figure . Effect of Temperature Level on Heat Pipe Performance 
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INSULATION 



Figure 3~^5- Schematic of the Test System 
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Figure 3-46. Overall Temperature Drop vs 
Power Transported 
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Figure 3~^2 indicates that an assembly using two heat pipes would be able 
to meet the power transport and redundancy requirements up to a heat pipe 
length of about 1.525 m (5 - f t) . On the other hand, Figure 3 -46 indicates that 
if two heat pipes are used and one fails the resulting heat pipe AT would range 
from 15 (27) to 23.9°K (43°F) . Adding to this value the temperature drops at 
the various interfaces the overall temperature drop (from sump to radiator) 
would significantly exceed the maximum value (16.68 K (30 F) ) selected on the 
basis of the sink temperature study. 

Hence, a four heat pipe assembly appears to be optimum. In case of a 
50-percent failure each heat pipe would transport 35 - w and the heat pipe AT 
would range from 7*22 (13) to 11.67°K (21 F) , depending on the diameter 
selected . 

Finally, a heat pipe diameter of 0.00952 m (0 . 375” in.) would result in a 
lighter heat p.ipe assembly and at the same time can meet the power transport 
and temperature drop requirements. Furthermore the small radius of curvature 
of the heat pipe at the sump (0.0318 m ( 1 . 25“ in.)) favors the smaller heat pipe 
diameter . 

Hence, a heat pipe assembly consisting of four aluminum heat pipes of 

0. 00952 m (0.375-in.) O.D. and using anhydrous ammonia as a working fluid has 
been selected. 

Heat Pipe Assembly Design 

The four heat pipes are mounted in an aluminum block which is made in two 
halves, each containing two heat pipes. The basic design approach of the 
aluminum block is similar to the copper assembly employed on the 5~w refrigera- 
tor. Aluminum was selected for its weight advantage. At the other end of the 
heat pipe assembly the individual heat pipes interface with the space radiator. 

1 . Aluminum Block Design 

The aluminum block design, depicted in Figure 3 -i *7 employs an interface 
clamp design which is adaptable to use with either a water cooling system or 
an ambient heat pipe assembly.. Either the cooling coil or the individual heat 
pipes are soft soldered to the aluminum block, as shown in Figure 3 -i *7, thus 
allowing the interchange of either system. The soft solder provides good 
thermal contact between the heat pipes and the aluminum and, also, allows heat 
to be distributed over most of the pipe circumference. 

The two halves are mounted on the engine sump with indium foil (0. 0000762 m 
(0.003-in.) thickness) at the interface and four bolts located in slots as shown 
in Figure 3 _ ^7. This scheme provides good thermal contact between the block and 
the sump while providing detachab i 1 i ty of either the heat pipe assembly or the 
cooling water assembly. 

The bolts use compression washers under each bolt head and each nut to 
achieve an interface pressure in the range of 3*^77 x 10^ (50) to 6.895 x 10"* 

N/m (100 psi) (between the engine sump and aluminum block). The fluctuations 
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HEAT PIPE SOFT SOLDERED 




Figure 3-47. Aluminum Block Design Details 
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in the engine sump pressure will result in an increase or decrease of the 
deflection experienced by the washers while maintaining the interface pressure 
in the desired range. 

Interface Temperature Drop 

1 . Heat Pipe Operation 


The temperature drop at the various interfaces 

are added to the heat pipe temperature drop to give 

drop in the ambient heat pipe assembly, AT ... 

7 overal 1 


and in the aluminum block 
the overall temperature 
Hence , 


^overall = AT i,, + AT A1 * iT ,, 2 + AT + AT i,j 


(3-52) 


where AT. , 
i J 


temperature drop at the interface between the sump and the 
aluminum block 


AT^ = temperature drop in the aluminum block 

AT. 2 = temperature drop at the interface between the aluminum block 
’ and the heat pipe 

AT = temperature drop in the heat pipe 



temperature drop at the interface between the heat pipe and a 
space radiator 


The thermal resistance at an indium foil interface was based on the data 

in Reference 3 - ^- for interface pressure between 3.^77 x 10"* (50) and 
5 2 

6.895 x 10 N/m (100 ps i ) the interface thermal resistance has the value 
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which gives 

AT. 1 = 0.0556°K (0.1°F) 

The soft solder joint thermal resistance was obtained from the data by 
Yovanovich and Tuarze, Reference 3 - 5- Hence, 


R. 9 - R. = 1 .505 x 10 


-5 °K - m 2 
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(o.wfuai) 
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Assuming two heat pipes to be carrying the entire 70 w results in 


AT. 9 2 = 0.33**°K (0.6° F) 
and AT. 3 = 0.278°K (0.5°F) 

The temperature drop at the radiator interface assumes a condenser length 
of 0 . 6 1 m (2-ft) and that the radiator panel interfaces with one side of the 
heat pipes over an arc length of 0.00025** m (0 . 1 - i n . ) . 

Finally the aluminum block temperature drop, under failure of 2 heat pipes, 
is estimated to be 

AT A1 = 5.36°K (10°F) 

Hence, the overall AT in the ambient heat pipe assembly is estimated to 
have a maximum value of 

^overall ~ 16.4°K (29.5°F) 

2. Water Cooled „ Operation 

The design of the aluminum block is such that water cooling coils may be 
substituted directly for the four heat pipes. The water cooling coils will be 
of the same diameter and geometric configuration as the evaporator end of the 
heat pipes. 

An analysis of the water cooling coils was performed similar to that for 
the 5 watt VM . The analysis indicated that the cooling water must pass through 
the four sections of the water cooling coil in series. The water flow was 

assumed at 0.00315 m /sec (50 gal/hr). The maximum temperature drop from the 
sumg gas to the cooling water with this configuration is approximately 6.67°K 
(12 R) , which is less than the temperature drop during heat pipe operation. 

MECHANICAL DESIGN 

The mechanical design of a long life Vuilleumier refrigerator requires 
detailed attention to every part in order to assure meeting the design require- 
ments. This section presents a description of the various parts of the machine, 
and detailed discussions of the interfacing requirements with the Honeywell 
radiometer and detector assembly. The preliminary design stress analysis is 
also discussed, as are the machine closure techniques. 

Design Description 

The preliminary design generated during Task 1 has drawn heavily on the 
experience gained during the design and development of a similar device for 
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GSFC. Under a former contract to GSFC, an engineering prototype of a 5-watt 
Vuilleumier Cryogenic Refrigerator was developed that met requirements very 
similar to those required under this current contract. In order to enhance the 
reliability and performance characteristics of the fractional watt flight proto- 
type refrigerator, similar or identical design approaches to the 5 watt machine 
have been employed. 

Description of the functional characteristics may best be accomplished by 
treating each detail of the design independently. The basic geometry, function, 
and materials of construction of each detail are described; these components 
are also related to the overall device. 

1 . Hot-End Insulation System 

Thermal containment of the heat supplied to the hot end is accomplished 
by applying Johns Manville Min-K insulation spirally wrapped around the hot- 
cylinder followed by a thin layer of AA-1200 fiberglass to fill out the blanket- 
ing to the outer protective shell. This shell will be fabricated of a corrosion 
resistant alloy. The cylinder will be terminated at one end in an ellipsoidal 
dome and at the opposite end in an O-ring-weld flange. This protective vessel 
will be attached to an evacuation system, thereby allowing the reduction of the 

-A 

pressure within the insulation system to a level of 10 torr or less to provide 
optimum thermal performance.. 

.2. Heater-Heat Exchanger Complex 

The introduction of heat into the VM hot end will be accomplished by a 
Cal-rod type of heating element. These devices are inherently long lived and 
reliable when properly installed. The heater is to be brazed to the elliptical 
dome transition from the cylindrical regenerator to the hot displacer hot bear- 
ing housing. Inside this elliptical surface a channel type forced convection 
heat exchanger has been provided to direct the flow into and out of the hot dis- 
placed volume of the refrigerator. The Cal-rod unit is composed of a corrosion 
resistant outer sheath with the heating element contained in a swaged magnesium- 
oxide insulator. The pressure wall dome is an integral part of the hot end 
pressure vessel and is fabricated of heat treated Inconel 718. The heat 
exchanger liner also will be fabricated from heat treated Inconel 718. 

3 • Hot Displacer Hot Bearing System 

The hot displacer is supported by bearings located at each end. To provide 
the long life sliding bearing required at the hot end of this refrigerator, Boe- 
ing Compact 6-84-1 is to be utilized working against a chrome carbide-nickel chrome 
matrix flame sprayed surface. A similar system was utilized in the 5 _ watt VM 
refrigerator. The Boeing Compact is assembled into a journal bearing with a 
metallic sheath for protection. This is then captured in the Inconel 718 pres- 
sure vessel housing and maintained in axial position by a high rate spring and 
retainer system. The chrome carbide mating surface is provided by flame spray- 
ing directly onto the journal provided by the displacer. This is then finish 
ground to the appropriate dimension. 
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4. 


Hot Regenerator 


The hot regenerator is formed by the installation of precise sheared 100 
mesh stainless steel screen discs. These are bounded on the outside by the 
heat treated Inconel 718 pressure vessel wall and on the inside by the regen- 
erator inner liner. The liner is formed by extending the hot end heat exchanger 
into the cylindrical section thereby forming the inner diameter junction for the 
regenerator. This element simultaneously forms the outer wall for the displacer 
bore and the reference cylinder for the displacer labyrinth seal. The cool end 
of the hot regenerator is terminated by a flow distribution plate which uni- 
formly directs the flow into and out of the regenerator from the sump heat 
exchanger section. This distributor plate also retains the screen disks both 
during assembly and operation of the regenerator. 

5- Hot Displacer 

This component is to be fabricated from heat treated Inconel 718 to pro- 
vide the strength required and also to assure a compatible coefficient of expan- 
sion with that of the regenerator inner liner. The displacer is supported by a 
bearing system at each end. The hot end was described in the hot bearing dis- 
cussion. The sump end will be similarly supported by a Boeing Compact 6-84-1 
s 1 i d i ng • jou rnal bearing mating with a tungsten carbide surface. This surface 
is to be flame sprayed onto the Inconel 718 displacer base material, similar to 
the technique utilized on the hot end, then finish ground to size. 

The displacer itself is composed of two pieces, welded at the central 
support rib. This allows machining of the two pieces to provide a minimum 
weight, maximum strength system. Convective heat transfer within the displacer 
is reduced by the inclusion of fibrous quartz prior to the closure weld. 
Labyrinth seals are to be integrally machined at the hot end of the cylindrical 
section to affect a pneumatic seal for the helium working fluid. 

6. Sump Section 

The region from the hot regenerator to the cold regenerator has been 
designated the sump region. One of the major components contained in this 
region is the ambient heat exchanger. This is formed by brazing an offset 
copper fin to the ID of the Inconel 718 pressure vessel cylindrical section. 

The inner wall of this heat exchanger is formed by the filler block installed 
in the sump area to reduce the dead volume. The filler block is to be fabri- 
cated from high strength aluminum alloy, and will enclose and support the cold 
end bearings for the hot regenerator. Flow passages are to be provided in the 
filler block leading from the cavity behind the hot displacer to the sump heat 
exchanger. The crankshaft area is also ported to the heat exchanger. 

7 • C rankshaf t 

High strength heat treatable alloy will be used to fabricate the multi- 
piece crankshaft. The concept developed for use in the fractional watt VM 
is quite similar to that utilized in the 5 watt refrigerator fabricated under 
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another GSFC contract. For this application, individual cam members and main 
bearing journals are to be pinned together for positive orientation and cen- 
trally retained by a high strength locking screw. This technique was adopted 
in order to allow installation pf the crankshaft — connecting rod complex from 
a central access port. Support of the crankshaft is to be accomplished utiliz- 
ing Boeing compact 6-84-1 sleeves running against tungsten carbide flame sprayed 
on the crankshaft journals and ground to size. Similarly, the connecting rod 
journals will be flame sprayed with tungsten carbide and ground to size. The 
connecting rod bearing mating surfaces will also be composed of solid Boeing 
Compact 6-84-1 journals pressed into the connecting rod ends. 

Axial positioning of the crankshaft for alignment is to be accomplished 
by a system of shimming and finally retained in the proper position by a high 
rate spring system. This latter device is to be incorporated to facilitate 
the installation and removal of the auxiliary break in magnetic drive without 
disturbing the critical location of the crankshaft. 

8. Connecting Rod-Wri stpi ns 

The concept of two piece split connecting rods will also be adopted from 
the 5 watt design. This is to allow final assembly from a central access port. 
Material is to be titanium in order to provide the required strength and also 
light weight for reduced rotating mass. The wrist pins are to be made of solid 
tungsten carbide and these will be running against Boeing Compact 6-84-1 
journals which are an integral part of the connecting rod end. 

9. Cold Displacer . 

The construction of the cold displacer is similar to that of the hot. 

The displacer body is to be fabricated from four sections, welded and finish 
ground to size. The displacer utilizes the cantilever concept of support to 
eliminate the need for the use of low temperature bearings. The bearing sur- 
face is to be flame sprayed tungsten carbide, ground to size, running against 
Boeing Compact 6-84-1. As with the hot displacer, the material is to be 
Inconel 718 in order to provide the expansion compatibility necessary in a 
device which must operate through high thermal excursions with small operatinq 
clearances. The cold end labyrinth seal is integrally machined into the warm 
end of the displacer. The interior of the displacer is packed with fibrous 
quartz. 

1 0. Cold Regenerator 

The use of monel shot 0.000178 and 0.000254 m (0.007 and 0.010 inches) 
in diameter will provide the maximum regeneration with a minimum of axial con- 
duction. Each end of the regenerator will employ 150 mesh monel screens to 
retain the ball packings and to provide radial flow distribution. The inner 
wall of the regenerator is formed by a thin wall Inconel 718 tube which also 
acts as the stationary portion of the labyrinth seal. 

A flow distribution annulus will be provided at the cold end of the 
regenerator and a radial slot distributor is to be utilized at the warm end. 

Both of these devices are provided to assure uniform radial and circumferential 
distribution of the helium fluid and also to prevent channeling. 
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The outer wall of this annular regenerator is formed by the pressure 
vessel wall, which in turn is the inner vessel for the evacuated dewar insula- 
tion system. Here also, Inconel 718 will be used to provide the lightweight 
high strength system necessary for containment of the 1000 ps i working pressure 
of the hel ium. 

1 1 . Cold End Insulation 


The use of two concentric rhodium plated radiation shields will assure 
adequate thermal protection of the low temperature section of the refrigerator. 
The cylindrical shields will be supported between the cold regenerator outer 
cylinder and the dewar outer cylinder by fiberglass pads. Provisions have 

been made to allow evacuation to the required 10 ** torr for optimum thermal 
protection performance. 

1 2 . Cold End DCA Interface and Heat Exchanger 

The cold end heat exchanger is located at the cold end of the cold regen- 
erator, and is composed of axial flow slots. The refrigeration heat load from 
the detector chip assembly (DCA) must be imparted to the gas at this location 
with the minimum temperature drop possible. This will be accomplished by a 
flat surface interface, heavily copper plated for low thermal resistance, and 
ground to an exacting surface flatness. The DCA is then to be held in place 
on this surface by three screws uti 1 i zi ng belleville washers for constant, 
uniform loading of the interface. To enhance this interface thermally, thin 
foil indium will be used to assure uniform, repeatable heat transfer character- 
istics. The final interface requirement between HRC and AiResearch equipment 
is that of the lead matrix. This element will be installed by HRC at the time 
of the DCA integration and will be soldered to the upper flange of the vacuum 
dewar outer vessel. For temporary evacuation and evaluation at AiResearch, 
a cold end closure system will be used, which will be described in a later 
portion of this report. 

1 3 • I nst rumentat i on 


The instrumentation required to monitor the performance of the VM may be 
divided into two categories: (l) Development Test and diagnostic instrumenta- 

tion and (2) flight performance monitoring equipment. In order to provide the 
tools necessary for initial performance analysis it is anticipated that the VM 
will be fitted with a complete complement of diagnostic instruments for labora- 
tory evaluation. A listing of the instruments and their function is presented 
in Table 3 “ 1 6 and also illustrated in Figure 3"48 to show the relative locations. 
For a flight configuration refrigerator, eleven of these devices would be 
retained or installation made compatible with telemetry signal cond i t i on i ng . 

These would include: 


1. Sump temperature 

2. Hot-end temperature 
3 - Cold end temperature 
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TABLE 3-16 

FRACTIONAL WATT VH DEVELOPMENT TEST INSTRUMENTATION 


Des i gnat ion 


Function 

Measurement 


Devi ce 


Location 


T l 

Temperature 

T 2 

Temperature 

T 3 

Temperature 

T 4 

Temperature 

T 5 

Temperature 

T 6 

Temperature 

T 7 

Temperature 

T 8 

Temperature 

T 9 

Temperature 

T ,0 

Temperature 

T ll 

Temperature 

T 1 2 

Temperature 

T I3 

Temperature 

T 1 4 

Temperature 

T l 5 

Temperature 

*,6 

Temperature 

T 1 7 

Temperature 

T 1 8 

Temperature 

T 1 9 

Temperature 

T 20 

Temperature 

T 2 1 

Temperature 

T 22 

Temperature 

T 23 

Temperature 

T 24 

Temperature 

T 25 

Temperature 

T 26 

Temperature 

T 27 

Temperature 


Thermocouple Cu-Con 


Thermocouple Cu-Con 
Thermocouple Ch-Al 


Thermocouple Ch-Al 


Simulated cold end load 
temp 

Simulated cold end load 
temp 

Cold end interface flange 
Cold end heat exchanger 
Cold end heat exchanger 
Cold regenerator, cold end 
Cold regenerator, cold end 

Cold regenerator, mid 
cy 1 i nder 

Cold regenerator, mid 
cy 1 i nder 

Cold regenerator, warm end 
Cold regenerator, warm end 

Sump heat exchanger, cold 
end 

Suipp heat exchanger, cold 
end 

Heat pipe block, cold end 
Heat pipe block, mid block 
Heat pipe block, warm end 

Sump heat exchanger, warm 
end 

Sump heat exchanger, warm 
end 

Hot regenerator, cold end 
Hot regenerator, cold end 

Hot regenerator, mid 
cy 1 i nder 

Hot regenerator, mid 
cy 1 i nder 

Hot regenerator, hot end 
Hot regenerator, hot end 
Hot end heater surface 
Hot end heat exchanger 
Hot beari ng surface 


EI i 

El, 

t 

El, 


Power 

Power 

Power 


Watt meter 
Watt meter 
Watt meter 


Cold end simulated power in 
Drive motor power i n 
Hot end heater power in 


PT, 


Interna 1 Pressure 
Vacuum pressure 
Vacuum pressure 


T ransducer 
Ion gauge 
Ion gauge 


System helium operating 
pressure 

Cold end insulation 
vacuum pressure 

Hot end insulation vacuum 
pressure 


Speed 


Magnet i c monopole 


Engine speed 
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Four heat pipe temperatures 

5. Internal pressure 

6. Motor speed 

7. Hot -end power 

8. Motor power 

14. Auxiliary Bearing Break-In Drive System 

During the break-in period of the NASA/GSFC 5 watt VM, it was found that 
the power required to drive the machine greatly exceeded the original design 
goal. Following the break-in period, the motor power is expected to return to 
the design level, based on the power reductions observed during the first 
1000 hr of break-in. For this reason, it was deemed necessary to provide a 
method of breaking in the fractional watt design without overdesigning the 
airborne motor, which would require excessive power consumption during the 
entire flight. Therefore, a laboratory break-in motor coupling was designed 
(Figure 3“49). This design utilizes the magnetic coupling concept originally 
developed for the 5 watt VM engine, as adapted quite simply to the fractional 
watt VM. To utilize the break-in drive system, the cover on the outboard crank- 
shaft access port is removed from the fully assembled engine. Removal of this 
cover in no way affects the position or slignment of any component within the 
machine since all alignment is maintained by the snap ring located in the bore 
of the main support tube. Removal of the cover exposes the outboard end of the 
crankshaft. An auxiliary drive shaft is attached to the crankshaft with three 
screws. A filler block and a magnetic rotor are mounted on the shaft. A spe- 
cial break-in dome is then placed over the magnetic rotor and the threaded re- 
taining ring installed, thereby resealing the internal volume of the engine. An 
externally driven magnetic rotor is then installed (quite similar to the design 
utilized in the 5 watt VM) which in turn is connected to a laboratory drive 
motor which has approximately ten times the torque capability of the airborne 
motor. 

An additional design feature afforded by this approach is the ability to 
operate the refrigerator with either the laboratory or the airborne motor. 
Ooeration from either motor may be selected without disturbing the test setup 
or opening the sump section of the VM. When the running torque has reached the 
design level, the system may be switched from the laboratory . to the airborne 
motor during operation. In the event that the additional inertia of the break- 
in magnetic coupling precludes starting with the airborne motor alone, startup 
may be accomplished by utilization of both motors. The laboratory motor will 
be disengaged when the refrigerator has achieved nominal running speed. 

Once engine break-in is complete, the entire magnetic drive mechanism is 
removed and the normal "flight" bearing access cover installed. In this way, 
the flight motor may be designed to produce only that torque necessary to drive 
the engine at the optimized minimum power level. 
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Figure 3"49- Break-In Magnetic Coupling 
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1 5. Radiometer and DCA Interfacing Techniques 


Several areas require careful design consideration in integrating the VM 
refrigeration package into the Honeywell radiometer. During the course of the 
Task I activities, several meetings, telephone conversations and communications 
were exchanged between AiResearch and Honeywell Radiation Center personnel. 

The items of information may be grouped into five categories and are discussed 
below by major topic. 

lb. VM-DCA Interface 


Early discussions regarding the Detector Chip Assembly (DCA) centered 
around the needs of Honeywell Radiation Center (HRC) to maintain the detector 
in a precisely located axial and radial position. This leads to the existence 
of two distinct areas to be considered: (1) axial and radial location stati- 

cally, and (2) axial and radial location under operating vibratory conditions. 
This latter facet is to be treated separately and is discussed under the cold 
end vibration testing of the 5 watt VM refrigerator. 

Initial requirements set forth by HRC dictated the axial prepositioning 
of the DCA within ±0.000025** m (±0.0001 in.). AiResearch studied several 
interface configurations for the DCA and the VM but it became obvious that the 
probability of being able to maintain the dimensional location of the DCA 
interface within the limits required by HRC would be quite remote, due to the 
fact that a hard interface point (a ledge or shoulder), ground to a very precise 
dimension, would be required'. The DCA would then have to be clamped into place 
utilizing a heat transfer improvement mechanism such as indium foil. Indium 
foil in itself caused problems in each of the designs considered. First, the 
thickness of indium foil varies from batch to batch, as well as within a given 
sheet, more than the allowable 0.0000508 m (0.0002 in.) required by HRC. 

Second, the relative surface condition of both the VM and the DCA could 
influence the degree of indium foil flow in order to accomplish the required 
thermal interface. Third, the quantity of flow of the indium foil as the 
result of thermal cycling, coupled with constant mechanical load, would produce 
a potentially changing position of the DCA. 

The results of this AiResearch study were relayed to HRC, with a request 
that the requirement for an axial interface be reconsidered. It was agreed 
that no additional effort be expended on the axial concept and that a radial 
scheme would be incorporated. Based on this agreement, several conceptual 
layout sketches were generated. Figures 3“50 through 3"5** illustrate several 
of the approaches considered. Illustrations 3 “ 50 through 3 ~ 5 3 all had one 
common undesirable featu re--name 1 y the indium foil has to be installed on the 
cylindrical surface and a close fitting copper clamping collar installed over 
the cold end of the VM and the indium foil. These techniques would, however, 
provide the axial dimensional prepositioning of ±0.000025** m (±0.0001 in.) 
required by Honeywell. It was felt that the attractiveness of any of these 
approaches was considerably decreased by the problem foreseen with the actual 
installation of the indium foil between the DCA and the VM wherein the foil 
could be crumpled during installation. Therefore, the method would never pro- 
vide a positive interface area that could be repeatedly attained. 
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DCA-VM Radial Interface - Concept 








Figure 3“52. DCA-VM ftadial Interface * Concept 











A method of eliminating this objection is shown in Figure 3“54, but does 
require precise control of the outside diameter of both the DCA and the VM to 
assure the uniform interface pressures necessary for adequate heat transfer 
while simultaneously assuring the axial location of the DCA within the 
0.0000254 m (± 0.0001 in.) tolerance. 

HKC was again contacted and apprised of the problems encountered. It was 
finally agreed to eliminate the axial tolerance requirements in favor of the 
axial interface system. Based on this, several designs were generated at 
AiResearch in order to define the most appropriate approach. Figure 3 ~ 55 repre- 
sents the design adapted for future analysis. This design too was modified 
during the final design layout effort. 

1 7 . Electro Magnetic Interference Compatibility 

Considerable concern was expressed by HRC regarding the EMI compatibility 
of the VM refrigerator. A specification and test method was provided to 
AiResearch, portions of which are included herein for reference: "The EMI 

levels within the dewar IR detector assembly are not to exceed 0.06 microvolt 
RMS in any 15 sec period, nor are spikes exceeding 0.2 microvolt P-P to occur 
more than 1 percent of this period. The frequency spectrum within which these 
measurements are to be made is to extend from 10 Hz to 155 Hz. The EMI levels 
outside the dewar are not to exceed 0.6 microvolts in a 15 sec period within 
the above spectrum nor are spikes exceeding 2.0 micro volts to occur more than 
1 percent of this period. 

The test methods are specified in MIL-STD-462 as applicable. Methods for 
bands below 10K Hz are to be submitted for approval. The probe is to consist 
of a single loop of wire located between the detector mounting position and the 
lead matrix array at right angles to the centerline of the cold cylinder. A 
low noise preamp - 3d b wave optimized for a 20 ohm source is to be connected to 
the loop." 

These requirements were reviewed with the appropriate personnel at 
AiResearch to determine the impact on the engine design. The following points 
resulted from this review: 

(a) The motor enclosure is adequate to eliminate any radiated EMI . 

(b) Leads to the motor must be shielded. 

(c) Cal rod heater is EMI free although leads must be double shielded from 
the GSFC supplied controller. 

(d) Test frequencies go down to 30 Hz per Ml L-STD-461 ; the requested 
requirements would go down to 10 Hz. In their experience, this has 
never been done at AiResearch or local test facilities. It would 
require the acquisition of special test equipment in the form of a 
memory scope for data recording. 



AIRESEARCH MANUF ACTURING COMPANY 

l,,f r.«M< •• I 


72-8846 
Page 3-132 




During the coordination effort with GSFC and HRC, it was concluded that 
any testing for compliance with the EMI requirements would be accomplished at 
HRC with the VM integrated into the Radiometer. This is brought about by the 
unique testing techniques, levels, and instrumentation that is required, and 
which is available only at HRC. Every area of potential EMI problems will be 
evaluated and treated accordingly during the design phase of the VM . 

1 8 . Mounting of the VM in the Radiometer 

During early interface negotiations, HRC presented their desired method 
of integrating the VM into the radiometer. Figure 3 _ 56 illustrates the 
requested method. This represented a drastic deviation from the proposed 
method of integration of the VM since the VM had been designed for two or 
three point mounting. This technique was utilized in order to render the 
lightest system possible and to avoid amplification of Radiometer input dynamic 
loads to the VM. 

Three techniques for incorporating the single flange concept of mounting 
were evaluated and are illustrated in Figures 3 ~ 57 through 3 “59 - Each approach 
required the use of an exoskeleton or significant reinforcement of the mounting 
flange systems in order to provide a system stiff enough to allow the VM to 
function in a normal manner. 

After careful consideration of the numerous alternatives available, it 
was mutually agreed to utilize a two point pickup of the VM at the radiometer. 
This technique provided the least weight with the maximum package structural 
rigidity. Figure 3 “60 illustrates the method of interfacing the second mount 
point that has been adopted. 

1 9 . Cold End Motion 


During the April 1 1 , 1972, interface coordination meeting at GSFC, it was 
determined that it would be necessary to establish the relative motion of the 
VM cold end with respect to the mounting flange interface. This information 
was deemed necessary in order to define whether the Honeywell sensor chip could 
be integrated directly into the AiResearch VM cold end or whether it would be 
necessary to devise an isolated platform between the VM and the sensor. 

GSFC requested the use of their 5 watt VM (which was undergoing performance 
testing at AiResearch) for measuring the relative motion, rather than relying 
solely upon calculated approximations. Several methods of instrumentation were 
evaluated, including the following final candidates: 

(1) Mechanical indicating methods 

(2) Accelerometers 

(3) Proximity probe indicators 
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Figure 3-56. HRC Proposed Radiometer-VM Interface 
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It was determined that the most accurate method would be the use of the 
proximity coils. These devices inherently have the capability of measuring 
position deviation on the order of 25 millionths of an inch and have very high 
response and resolution character i s t i cs . 

Figure 3 - 6 1 illustrates the instruments which were installed on the GSFC 
5 watt VM cold flange. This test setup was to provide two radial modes (90° to 
one another) and one axial pickup. 

In order to install the sensor system, it was necessary to remove the 
vacuum dewar jacket, and the multilayer superinsulation. All thermometry 
instrumentation and load heater leads were maintained intact to eliminate 
possible performance deviations during subsequent performance evaluation 
tests. Prior to installing the proximity probes in the holding device, each 
probe was individually calibrated utilizing a precision micrometer capable of 
accurate readout to 50 millionths of an inch. The calibration test setup is 
shown in Figure 3~62. The curves generated by the X-V plotter are illustrated 
in Figures 3~63 and 3~Sk. Figure 3 ”6 5 depicts the entire instrumentation 
system as it was installed. 

During the running of the vibration tests it was found that the dimensional 
deviation from nominal was so small that only a small portion of the calibration 
curve was useable. Each curve was examined in order to obtain the most linear 
portion and expanded to produce the curve presented in Figures 3- 66 and 3 “67 . 

The engine was. then run at 300 rpm and the resulting motion recorded on an 
oscilloscope via the proximity probes. The oscilloscope traces were photo- 
graphed and are reproduced in Figure 3 ”68 for evaluation. Prior to starting 
the rotation of the VM engine, a zero speed or background noise measurement was 
taken and is also presented in Figure 3“68. From the relative motion noted in 
the two probes perpendicular to the displacer centerline it was determined that 

a deviation of ±1.9 x 10 ^ inches was noted with one while ±6.5 x 10 ^ inches 
was noted with the second. No relative motion could be sensed with the third 
probe (axial motion of the pressure cylinder) in as much as no heat was applied 
to the hot end and therefore there was no pressure variation in the system. 

The results of the 5 watt displacement measurements were analyzed and the 
following conclusions reached. 

(a) The tests were run with the machine operating at 300 rpm but power 
was not applied to the hot end. The tests indicated a relative dis- 
placement of a few millionths of an inch. This is well below the 
Honeywell described allowable displacement (by 2 orders of magnitude) 
and thus problems are not anticipated for the fractional watt 
refrigerator. 

(b) The test results on the 5 watt VM cannot be considered completely 
conclusive with respect to the fractional watt machine since the size 
is different and the tests were not conducted at operating tempera- 
tures and corresponding internal pressure cycling. 
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Figure 3*63. X-Y Plotter Curve: for Probe S/N 4111450 
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Figure 3-67. Expanded Curve for Probe (SN 41X312) 
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(c) Due to (b) above, it is anticipated that the GSFC will require dis- 
placement measurements to be made on the fractional watt machine at 
operating conditions to assure a compatible interface with the 
Honeywell detector assembly. 

Cold End Closure Dome 


During the numerous interface meetings with Honeywell, it was determined 
that the interface for the cold end would exist at the flange provided on the 
VM cold cylinder for mounting of the Honeywell detector. In addition, the 
vacuum dewar is to also terminate in a flange compatible with installation of 
the Honeywell provided lead matrix ring. Both of these interface planes 
require the design and installation of auxiliary mechanisms if the VH is to 
be tested at the AiResearch facilities, prior to shipment to the Honeywell 
Radiation Center. (Integration of HRC equipment into the VM cold end and 
installation of the VM into the radiometer will occur at Honeywell.) For this 
reason, two laboratory test devices have been designed to allow operation at 
AiResearch : 

(1) In order to provide a source of thermal load comparable to that of 
the Honeywell detector and detector mount base, a dummy detector has 
been designed incorporating an electrical resistor to simulate the 
load. This will be installed on the cold heat exchanger at the end 
of the cold cylinder. 

(2) The vacuum dewar will then be completed by the installation o'f a 
temporary closure to allow evacuation of the annulus. This device 
will be held in place by a ring clamp and utilizes an 0-ring seal 
at the interface to the vacuum dewar cylinder. Provisions for an 
evacuation tube will be incorporated along with a pinch-off tube 
cover. With this design, when testing has been completed at 
AiResearch, the dummy load may be removed, the evacuation closure 
installed, and the system evacuated and pinched off for shipment. 

This will protect the cleanliness and efficiency of the superinsula- 
tion during transportation and storage. 

Figure 3“69 depicts the cold end vacuum dewar test closure and the dummy 
load interface mechanism. 

Final Closure Technique 

Early in the design of the fractional watt VM it became clear that it 
would be necessary to devise a method of closure of the various components that 
was unique. The closure system requirements were centered around the require- 
ment for a flight weight system while still maintaining the "prototype" capa- 
bilities of easy access for removal, replacement, or rework of any of the 
internal components. With this approach as a baseline, the two chief areas of 
concern were summarized as follows: 
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NOTE: EVACUATION TUBE SHOWN IN THE PINCHED 

OFF (SHIPPING) CONFIGURATION WITH 
PROTECTIVE CAP INSTALLED 



Figure 3-69. Cold End Vacuum Dewar Test and Shipping Closure 
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(1) The joint must be a "bolt together" design to allow easy access to 
the interior mechanisms while the engine is in the fi tup-break in- 
performance testing phase. 

(2) Once all internal components are correct and the appropriate perform- 
ance has been achieved, the unit is to be permanently sealed and the 
bolting mechanisms removed resulting in the lightest design possible. 
Leakage of the operating fluid is not to exceed approximately 1 x 

-8 

10 scc/sec helium at one atmosphere equivalent differential pres- 
sure. In addition, no contamination of the operating fluid is to 
occur as the result of the retention of any sealing mechanism or as 
the result of excessive outgassing of the sealing system during any 
material joining necessary to accomplish the required seal. 

Figures 3 “70 through 3“8l illustrate several of the concepts generated 
early in the design phase. In general, fault was found with each approach as 
shown with each illustration. The use of a breach-lock mechanism employing a 
double weld joint lead to the solution of the problem. 

After several design iterations, a concept for producing the final closure 
joints was developed. Figures 3"82 and 3“83 illustrate the initial concept as 
it applies to the three primary joints under consideration: the cold end; hot 

end, and crankshaft closures. This method utilizes a threaded ring that allows 
all mechanical and pneumatic loading to be carried by the threads. Two methods 
of sealing the closure may be utilized, dependent upon the particular mode of 
fabrication. Initially a "K" seal or metallic 0-ring could be utilized allow- 
ing easy access to internal components without having to destroy any of the 
closure mechanism. Following the required bearing and engine break in period, 
and after all required entry to the interior components has been accomplished, 
the threaded ring is removed and the sealing mechanism destroyed by punching 
one edge locally. The threaded ring is then replaced and the two small stand- 
ing edges are welded. This procedure results in a sealing system wherein all 
mechanical and pneumatic loads are carried by the threaded ring. The only 
load imparted to the standing edge welds is the internal pressure acting on 
the weld members themselves. It should be noted that the metallic seal was 
punched to assure that the seal was destroyed, thus providing a positive gas 
path to the welded area during helium mass spectrometer leak check, after 
welding is complete. It is currently anticipated that three closure tools 
will be utilized to properly compress the metallic seal so as not to rely upon 
the threaded ring to pull the two components together. 

The principles illustrated in Figures 3 “82 and 3 “83 were then further re- 
fined for closure joints in the sump-crankshaft region (Figure 3-84), the cold 
end of the engine (Figure 3 “ 85 ) and the hot regeneratoi — displacer complex 
(Figure 3 “ 86 ) . Additional design modifications will be incorporated during 
the second and third . tasks of this contract and as the result of the test pro- 
gram to be conducted during the task two activities. In order to fully define 
the test sequence and parameters to be evaluated, the test plan for the final 
closure sealing system has been included as Appendix B of this report. 
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CHANNELS FORHED 
OVER JOINTS 


SAW CUT 
FOR REMOVAL 



<t 

NEGATIVE POINTS: 

1. CHANNELS MUST BE CUT OFF THUS CAUSING METAL CHIPS THAT COULD FALL INTO UNIT 

2. RELIABILITY OF CLOSED CHANNEL TO TAKE FULL LOAD 

3. DIFFICULTY OF FORMING CHANNELS ON COMPLETED ENGINE THAT COULD DAMAGE ENGINE 

Figure 3-70. Formed Channel 



ALL RIVET HOLES & PIN HOLES MUST BE 
WELDED CLOSED AT FINAL WELDING 


NEGATIVE POINTS: 

I . POOR RIVET INSTALLATION ON CURVED SURFACE 

2. RIVETS MUST BE DRILLED OUT, THUS CAUSING METAL CHIPS THAT COULD FALL INTO UNIT 

3. MUST CARRY FULL LOAD 

Figure 3“ 71 • Riveted Assembly 
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ALL BOLT HOLES AND INDEX PIN HOLES MUST BE WELDED 
CLOSED AT FINAL WELDING. 


Figure 3-73* Split Ring Assembled with Bolts 


SPLIT RINGS 


FLANGE SCALLOPED 
TO ALLOW FOR 
PARTIAL WELDING - 
BEFORE REMOVING 
RINGS 



HOT DISPLACER 


Figure 3-7^- 


Removable Split Ring Fixture — Hot Displacer 
Housing Connection to Sump 
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SPLIT RINGS 


BAND CLAMP 



Figure 3~75. Removable Split Ring Fixture — Cold Displacer 
Vessel Connection to sump 






*L HOT DISPLACER. 


Figure 3~76. Bayonet Type Lock— Hot Displacer Housing 
Connect ion to Sump 


s -77552 
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HOT DISPLACER 


Figure 3"77» Removable Split Ring* Fixture Motor Housing. 
Connection to Hot Displacer Housing 


— C HOT DISPLACER 



<fc MOTOR 


S -77579 


Figure 3“ 78 . 


Bayonet Type Lock — Crankshaft Closure Dome 
Connection to Housing 
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<fc MOTOR 


5. HOT DISPLACER 


Figure 3 — 79 • Bayonet Type Lock — Motor Housing Connection 
to Hot Displacer Housing 



•fc COLD DISPLACER 


$-77578 

1 

Figure 3-80. Bayonet Type Lock — Cold Displacer Vessel 
Connection to Sump 
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REMOVABLE SPLIT RINGS 



- £ 

COLD DISPLACER 



S-77574 

Figure 3-81 . Removable Split Ring Fixture 

Cold Displacer Vessel Connection 
to Sump 
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Figure 3-86. Hot End Preliminary Design 




STRESS ANALYSIS 


I ntroduct i on 


During the preliminary design phase (Task I), stress analyses were per- 
formed to establish and/or verify the basic dimensions of the major components. 
The detailed stress analyses will be performed during Task II of this program. 

The preliminary stress analyses were concerned primarily with establish- 
ing the required wall thickness and stiffening ring spacing and geometry for 
the pressure shells, displacers, vacuum jackets, and crankcase housing. The 
relative motion of the cold end displacer and regenerator were also irivesti- 
gated for sinusoidal, random, and self-induced vibration. 

Design Criteria 

Although the preliminary work was performed with an operating pressure of 
7 2 

1.036 x 10' N/m (1500 psia), it became apparent near the end of task 1 that 
the heat leaks and weight penalties imposed by this pressure were too severe. 

A decision was therefore made to reduce the maximum operating pressure to 
6 2 

8.62 x 10° N/m (1250 psia). This analysis is based on the reduced operating 
pressure. The components subjected to working pressure were therefore analyze 
for: 

6 2 

Operating pressure » 8.62 x 10 N/m (1250 psia) 

7 2 

Proof pressure = 1 .67 x Operating Pressure = 1 .44 x 10 N/m 

(2060 psig) 

7 2 

Burst pressure = 2.25 x Operating Pressure = 1.94 x 10 N/m 

(2800 psig) 

The unit was designed to withstand proof pressure without permanent defor 
mation and burst pressure without rupture. 

All components were designed to target factors of safety of 1.15 on yield 
strength, 1.50 on ultimate strength, and 2.0 on elastic instability (buckling) 

For components subjected to high temperatures, the allowable stress was 
that value at which 0.2% creep occurs in 20,000 hr. The temperature profile 
was 922. 2°K (1200°F) at the hot end, 333.3°K ( 1 40° F) at the sump, and 65°K 
(-3^3°F) at the cold end. 

For components subjected to cyclic loads, the allowable stress for steel 

g 

components was based on a life requirement of 10 cycles. 

The ultimate design load factor was 30 g applied separately in any direc- 
tion. This load factor was not superimposed on the vibratory environment 
defined below. 
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The vibratory environment is as follows: 


Sinusoidal Vibration - Thrust Axis ZZ 
Frequency, Hz Level (Q-to-peak) Sweep Rate, Oct. /min 


5-11 

12.2 mm d .a . 
(0.48 in. d .a . ) 

2.0 

11-17 

±2.3 g 

2.0 

17-23 

±6.8 g 

1 .5 

23-200 

±2.3 g 

2.0 


Sinusoidal Vibration - Lateral Axes X-X and Y-Y 


Frequency, Hz 

Level (0-to-peak) 

Sweep Rate, Oct. /min 

5-9-9 

19.0 mm d .a . 
(0.75 in. d .a . ) 

2.0 

7.1-14 

±2.0 g 

2.0 

14-200 

±1 .5 g 

2.0 

Random 

Vibration - Axes ZZ , 

X-X, and Y-Y 

Frequency, Hz 

PSD Level , q 2 /Hz 

Accel., g-rms Duration 


20-300 0.0029 to 0.045 9-1 4 minutes 

increasing from 20 Hz each axis 

at rate of 3 dB/octave 


300-2000 0.045 9*1 4 minutes 

each axis 

Discussion of Results 
1 . Hot End Pressure Shell 

With inside radius of 0.0328 m (1.290 in.), the required cylinder wall 
thickness based on hoop tension is 0.000741 m (0.0292 in.). The wall thickness 
will be increased at the ends to carry the discontinuity stresses. The thick- 
ness of the elliptical dome is 0.00165 m (0.065 in.); required to maintain the 
operating stress level below the creep allowable at 922. 2°K (1200°F). 
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2 . 


Hot End Displacer 


The required cylinder wall thickness, stiffening ring spacing and geometry 

to resist an external pressure of 1.94 x 10 7 N/m 2 (2800 psi) were computed to 
be : 

• Cylinder wall = 0.001017 m (0.040 in.) 

• Ring spacing = 0.03655 m (1.438 in.) 

• Ring width = 0.00203 m (0.080 in.) (along cylinder (£) 

• Ring thickness = 0.00254 m (0.100 in.)(radial) 

3. Cold End Displacer 

The configuration required to resist an external pressure of 2800 psi was: 

• Cylinder wall = 0.000304 m (0.012 in.) 

• Ring spacing = 0.01347 m (0.530 in.) 

• Ring width = 0.000762 m (0.030 in.)(along cylinder (£) 

• Ring thickness = 0.001142 m (0.045 in.)(radia1) 

4 . Cold End Pressure Shell 

The cylinder wall thickness was governed by the requirement that the rela- 
tive motion between the displacer and the regenerator be less than the clearance 
provided. The estimated frequency of the displacer and the regenerator (shell) 
was 452 and 448 Hz, respectively. Since the frequencies are nearly equal, the 
motion of the displacer and shell are in phase and the relative motion is less 
than the clearance provided. The required cylinder wall thickness is 0.000609 m 
(0.024 inches) with a 0.0278 m (1.093 inch) outside diameter. 

The displacement at the end of the pressure shell due to self-induced 

vibration was estimated at 0.508 x 10 ^ m (20 x 10 ^ i nches) (doub le amplitude). 
This is an order of magnitude less than the allowable displacement of 5 . 0 8 x 
“6 -6 

10 m (200 x 10 i nches) (doub 1 e amplitude). The estimate was obtained by 
considering the unbalance of the displacers, connecting rods and crankshaft and 
driving the crankshaft at 400 rpm. 
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5. 


Crankshaft Housing 


The crankshaft housing was analyzed for internal pressure. Particular 
attention was directed to the housing/cross tube intersection and the spherical 
dome/tube intersection. The required wall thickness of the housing was found 
to be 0.001627 m (0.064 inches) and the required wal.l thickness of the cross- 
tube was 0.00292 m (0.125 inches), and the required wall thickness of the tube 
at the spherical dome end was 0.00226 m (O.O 89 inches). 

6 . Crankshaft Assembly 

A preliminary critical speed analysis was performed considering the effects 
of the displacers reciprocating linear unbalance as well as rotary unbalance. 

The first mode critical speed was found to be 8*400 rpm; significantly higher 
than the 400 rpm maximum operating speed. 

Analyses were also performed to verify that the crankshaft preload was 
sufficient to preclude separation of the stacked-parts interfaces and that the 
stress levels in the crankshaft and the connecting rod assemblies are at 
acceptable levels. 

MOTOR DESIGN 

The drive motor for the preliminary design VM refrigerator is a six pole 
squirrel cage induction motor driven by a 21 Hz, 26 volt, three phase power 
supply. Preliminary calculations indicate that the required power supply 
rating is 40 volt-amp (3 phase total) for rated load runnings, and 200 volt- 
amp for short term (less than a second) motor starting. The power supply 
should have both variable voltage and frequency output if running at different 
speeds will be required. The stator and the rotor laminations are made from 
hot rolled, low core loss, high permeability thin silicon steel to increase 
motor efficiency. The motor construction is mechanically rugged, with the con- 
ductor bars and end-rings integrally die case with the lamination stack. No 
brushes or slip rings are required. 

The stator winding utilizes 493°K (220°C) wi re and insulation materials 
to ensure long insulation life. The motor electrical and thermal design pro- 
vides for low loss density to keep the maximum hot spot temperature below 
493 K (220°C) . 

Motor performance, power factor, and efficiency is somewhat reduced by 
the increased air gap necessary to install an approx .• 0.000254 m (0.010 inch) 
thick bore seal. The selection of the bore seal material is based on the 
combination of high electrical resistivity and high mechanical strength. 

Possible materials are Inconel 718, or Ti6-AL-4V Titanium alloy. The estimated 
motor efficiency is approximately 25%, depending on the power supply wave shape, 
and the estimated power factor is 40% at the rated 2.5 watts shaft output load 
condition. The motor is running at approximately 5 % slip at the 400 rpm speed, 
rated load condition, and will maintain ±1% speed at the load for the expected 
temperature variations. 



AIRESEARCH MANUFACTURING COMPANY 


72-8846 
Page 3“ 1 67 




REFERENCES 


3~1 Browning, C. W. and V. L. Potter, "75°K Vuilleumier Cryogenic Refrigerator, 
Final Report for Task II, Analytical and Test Program," AiResearch Report 
72-8497, August 1972. 

3-2 B. D. Marcus, "Theory and Design of Variable Conductance Heat Pipes 
Hydrodynamics and Heat Transfer," Research Report No. 1, Contract No. 
NAS2-5503, TRW. 

3-3 ICICLE Feasibility Study, Final Report, NASA Contract NAS5-12039, RCA- 
Defense Electronic Products, Camden, New Jersey. 

3-4 L. S. Fletcher, P. A. Smuda and D. A. Gyorog, 'Thermal Contact Resistance 
of Selected Low-Conductance Interstitial Materials," AIAA Journal, Vol . 7, 
No. 7, p. 1302. 

3-5 M. M. Yovanovich and M. Tuarze, "Experimental Evidence of Thermal Resis- 
tance at Soldered Joi nts'S Journal of Spacecraft, Vol. 6 , No. 7, p. 856 , 
July 1969 . 



AIRESEARCH MANUFACTURING COMPANY 


72-8846 
Page 3-168 




APPENDIX A 


REGENERATOR ANALYSIS 


General 


The regenerator is one of the most important components of the Vuilleumier 
cycle refrigerator. It is a heat storage device which is used at both the hot 
and cold ends of the refrigerator. At the cold end, it is used to cool the gas 
as it flows from the sump region to the cold expansion volume of the machine. 
The expansion process of the gas reduces the gas temperature further in order 
to provide cooling at the cold temperature. After absorbing heat from the 
refrigeration load, the gas is then exhausted from the expansion volume through 
the regenerator where it regains energy previously stored, to the sump region. 
This reestablishes the temperature profile in the cold regenerator for the in- 
coming gas for the next cycle. 


The process in the hot regenerator is essentially the same as that at the 
cold end. The gas flowing from the sump region to the hot end is heated in the 
regenerator. As the flow reverses, an expansion process begins, which tends to 
cool the gas. This allows heat to be added to the gas at essentially constant 
temperature from the hot end heater. The gas is exhausted from the expansion 
volume through the regenerator where it transfers energy to the regenerator 
matrix and is cooled. This reestablishes the temperature profile in the hot 
regenerator for the incoming gas for the next cycle. 

The design requirements for a good regenerator are that the regenerator 
packing be of a material with a very large heat capacity and that the packing 
have a large heat transfer coefficient and a large heat transfer area. The 
basic tradeoff for the design of a regenerator is between irreversible pressure 
drop and heat transfer potential with a minimum void volume. 


The. analysis of the regenerators for a Vuilleumier refrigerator cannot be 
based on the classical effectiveness parameters which make use of inlet and out- 
let temperatures. The system pressure fluctuates and a considerable amount of 
gas is stored in the regenerator void volume during the flow period. This 
results in the mass flow of gas into the regenerator not being equal to the 
mass flow of gas out of the regenerator at all points in time. 

Computer Program 

In order to analyze the regenerators properly, a computer program making 
use of finite difference techniques has been developed by AiResearch. The 
regenerator is broken into axial nodes; these separate nodes are used for both 
the gas and the matrix. The following basic equations were used to develop the 
difference equations which the computer program solves. 

1 . Gas Nodes 


The continuity equation can be written as 

St Sx u 


(A-l) 
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where p = gas density 


t = time 
u = gas veloc i ty 
x = coordinate 


Setting 


G - o u 


oG _ Sq. 

Now writing the difference form 


the refore 



The momentum equation for the gas is 


du , du 4f 

p aT + pu ^ + 


r 2. dP 

2D, (pU > + §7 9 C " ° 

n 


where f = Fanning friction factor 

D^ = characteristic length 

P = static pressure 

q = 32.2 I bm- f t/ 1 bf-sec^ 

3 c 

Taking the steady-state form and substituting G, 


(A-2) 


(A-3) 


(A-4) 


op _ &(g /gj . I if_\ gJ_g[_ 


Sx 


2D, 


9 C 3x 

Now taking the d i f fe rence ' form of this equation, 

> 

n+ 1 


< p - - p - ( g „ + | 2 /p„ + |) - «-.V> 


_n 
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JI\T 


n n 


I 

+ 2 


fe) 


G . G . 
n+ I n+l 

Pn+I 
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(A-5) 


(A-6) 
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(A-7) 


Now solving for 



The energy equation is 

?(c pT) S(c puT) /UA„ \ 

z + E = [ — L j ( t - T) 

St Sx yi ) 

where c = specific heat at constant volume 
v 

Cp = specific heat at constant pressure 
T = gas temperature 
U = overall heat transfer coefficient 
k. - heat transfer surface area per unit length 
V - void volume in regenerator per unit length 
t - matrix temperature 
This equation can be reduced to 


. T ^ C n GT) 

;;f* — P — 

g- C)X 



(t - T) 


The overall heat transfer coefficient u is defined as 


whe re 


I 



j = Colburn j-factor 

P P randt 1 number 
r 

\ ~ equivalent conduction length 


(A-8) 


(A-9) 


(A-10) 
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Taking the difference form and solving for future gas temperature 



T + 
n 


(AT )(G c T - G c .T 
q n pn n n+ 1 pn+1 n+ 1 


Ax 

n 


vn 


P 

n 


(A-11) 


+ (AT g } | \ (t n " V 

c o 
vn n l 


The stability criterion for this difference equation is 



(A— 12) 


In the case that the time increment At required by the stability criterion 

9 

for the gas is very small compared to that of the matrix, steady-state solutions 
may be selected. In this case, the steady-state energy equation is 


o ( c GT) 

e_ 

Bx 



( t - T) 


(A-13) 


Writing the difference form 



2. Matrix Nodes 


T ) 
n 


The energy equation for the matrix is 



(A-14) 


(A-15) 


(A- 16) 
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where t matrix temperature 

C : matrix speci fie heat 
m 

= matrix mass per unit length 

k = matrix conductivity 
m 1 

A x = matrix effective heat transfer area for conduction 
Wi iting the difference equation and solving for future matrix temperature^ 



whe re 'HA), 

i 

The s tab i 1 i ty 
^ m 


= interface or boundary heat transfer characteristic, 
criterion for the matrix time increment is 



(A- 18) 


3. D i scuss i on 


The computer program is capable of analyzing both transient performance 
and cyclically steady performance. The inputs required by this computer pro- 
gram are the physical characteristics of the regenerator, heat transfer and 
friction loss characteristics, initial conditions, and boundary conditions. 

The physical cha rac t e r i s t i cs are reflected in the matrix areas, volume, length, 
heat capacity, thermal conductivity, and mass. The heat transfer and friction 
loss characteri st ics are read into the computer program in the form of Colburn 
j-factor and Fanning friction factor as a function of Reynolds number. The 
initial conditions must be fully described in terms of pressure and tempera- 
tures of the gas and the matrix. The boundary conditions that are required 
are the time dependent cha racter i s t i cs of pressure, mass flow rate, and gas 
temperature at one end of the . regenera tor, and the return gas temperature at 
the other end (basically defining the thermodynamic process at the other end). 
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Regenerator Cha rac t er i za t i on 


In the analyses of the regenerators, friction coefficients and the Colburn 
j- factor have been taken from the data presented by Kays and London. The data 
for randomly stacked screens are shown in Figure A- 1 and the data for regenera- 
tors packed with spherical shots are* shown in Figure A-2. Standard screen mesh 
data, shown in Table A-l, was used to characterize the screen matrix. Figures 
A-3 and A-4 show the specific heats of Monel, stainless steel, and Inconel as 
a function of temperature. 

I . Area t o Vo 1 ume Ratio 


For screens, the area to volume ratio can be derived to be 


0 = -o 

d 


(A-19) 


where 


- area to volume ratio 
6 - poros i ty 
d - wi re d i amete r 

(A-20) 


For spheres, the area to volume ratio is 


0 


6U - € ) 


where 


d - sphere diamter 


2 . Hyd rau 1 i c Pi ame te r 

Another regenerator cha racter i s t i c which is required is the hydraulic 
diameter. Defining the hydraulic radius as being equal to the flow volume 
divided by the total surface area, the hydraul ic radius of a screen stack is 

found to be 


ed 

H TV V) 


For a sphere, 


The hydrau 1 i c 


the hydraulic radius is equal to 

ed 

6 ( I -e ) 

diameter is equal to 4 times the hydraulic radius. 


(A-21) 


(A-22) 
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REYNOLDS NUMBER, 








TABLE A-1 


U.S. SIEVE SERIES AND 
TYLER EQUIVALENTS 

A.S.T.M. - E-l 1-61 


Sieve Designation 

Sieve Opening 

Nominal Wire Diameter 

S tanda r d 

A1 ternate, 


Inches, 


Inches, 

Tyler Equiv. 
Des i gnat ion , 


No. 

mm 

approx, equiv. 

mm 

approx, equiv. 

mesh 

1 . 4 J mm* 

14 

1 .41 

0.0555 

0.725 

0.0285 

12 

1-19 mm 

16 

1 . 19 

0.0469 

0.650 

0.0256 

14 

1.00 mm* 

18 

1 .00 

0.0394 

0.580 

0.0228 

16 

841 micron 

20 

0.841 

0.0331 

0.510 

0.0201 

20 

707 micron* 

25 

0.707 

0.0278 

0.450 

0.0177 

24 ’ 

595 micron 

30 

0.595 

0.0234 

0.390 

0.0154 

28 

500 micron* 

35 

0.500 

0.0197 

0.340 

0.0134 

32 

420 micron 

40 

0.420 

0.0165 

0.290 

0.01 13 

35 

354 mi.cron* 

45 

0.354 

0.0129 

0.247 

0.0097 

42 

297 micron 

50 

0.297 

0.01 17 

0.215 

0.0085 

48 

250 micron* 

60 

0.250 

0.0098 

0.180 

0.0071 

60 

2 10 mi cron 

70 

0.210 

0.0083 

0.152 

0.0060 

65 

17 7 m i c ron* 

80 

0. 1 77 

0.0070 

0.131 

0.0052 

80 

149 micron 

100 

0.149 

0.0059 

0. 1 10 

0.0043 

100 

125 mi c ron* 

120 

0. 1 25 

0.0049 

0.091 

0.0036 

1 15 

1 05 mi cron 

140 

0. 105 

0.0041 

0.076 

0.0030 

150 

88 micron* 

1 70 

0.088 

0.0035 

0.064 

0.0025 

1 70 

74 micron 

200 

0.074 

0.0029 

0.053 

0.0021 

200 

63 micron* 

230 

0.063 

0.0024 

0.044 

0.0017 

250 

53 micron 

270 

0.053 

0.0021 

0.037 

0.00i5 

270 

44 micron* 

325 

0.044 

' 0.0017 

0.030 . 

0.0012 

325 

37 micron 

400 

0.037 

0.0015 

0.025 

0.0010 

400 


*These sieves correspond to those proposed as an international (ISO) standard. 
It is recommended that wherever possible these sieves be included in all sieve 
analysis data or reports intended for international publication. 
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SPECIFIC HEAT, Btu/lb- 



TEMPERATURE, °R 

$-58038 


Figure A-3 • Monel Specific Heat 
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Ref: AiResearch Material Design Data 

347 STAINLESS STEEL: ASTM STP NO. 227 

INCONEL IOO: G.E. HANDBOOK, SEC. 2 
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Figure A-4. Stainless Steel and Inconel Specific Heats 




APPENDIX B 


FINAL CLOSURE SEALING WELDING 
TEST PLAN 


INTRODUCTION 

The design and function of a closure technique for sealing all entrance 
ports is critical to efficient operation of the fractional watt VM. The 
concept developed under this contract involves a metallic temporary sealing 
mechanism coupled with a permanent weld seal following all fabrication 
operations. An illustration of a candidate design is shown in Figure B-l. 
The attributes of this device may be summarized as follows: 

(a) Preliminary sealing is accomplished by a metallic sealing device 
(a" k" seal or metal "0 M ring) thus allowing repeated sealing 
without altering the components 

(b) All mechanical and pneumatic loads are carried by the threaded 
reta i ner 

(c) Following final entrance into the VM, the two pieces being sealed 
are welded together at the ring. The orientation and position of 
the two pieces is maintained by a typical weld restraining tool 
shown in Figure B-2. This method allows the complete freedom 
needed during a development program but also provides, at minimum 
weight, the reliable leak-tight sealing required of a 2 year flight 
prog ram 

These concepts must be confirmed by a comprehensive test program to 
establish complete compliance with the requirements of the overall project. 
This may be accomplished by subjecting simulated closure mechanisms to the 
conditions and requirements anticipated for the flight model VM. 

TEST SAMPLES 

Figure B-2 illustrates a typical tool necessary to restrain the two 
pressure vessel halves while the threaded ring is removed and replaced prior 
to the final cfosure weld. Figures B-3 and B-4 depict the hot-end closure 
a ssemb ly and the cold-end closure assemb ly respectively. 

In each instance, the exact configuration currently envisioned in the 
actual VM has been incorporated in the test models. In this way, the levels 
of stress that will be encountered in this test program will duplicate 
those that will be found in the actual engine. 

TEST PROCEDURE 

The test procedure for evaluating the acceptabi 1 \ ty of the closure 
concept will be the same for the hot-end and cold-end configurations 
illustrated in Figures B"3 and B-4. The test pieces will be fabricated 
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Figure B-l . Candidate Closure Design 
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HOT DI$PLACER' HOUSING 
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BREAK-IN DRIVE ASSEMBLY TOOL 


l 



Reproduced from 
best available copy. 


Figure B-2. Restraining Tools - Typical 
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Figure B-3- Hot-End Closure Assembly 
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in accordance with applicable engineering drawings and the metallic seals 
obtained from a commercial source. 

The test pair will be assembled, incorporating the metallic seals, and 
after assembly the test specimen will be mechanically loaded with the appro- 
priate assembly tool. Once the restraining ring is torqued into place, the 
assembly tool will be removed and the sealed assembly subjected to a helium 
mass spectrometer leakage test. Once it has been established that the leakage 

rate is 1 x 10 ^ standard cubic centimeters per second helium or less, with 
7 2 

1.034 x 10 N/m (1500 psig) applied to the specimen internally, the sealing 
pair will be taken apart and reclosed. The test setup diagram shown in Figure 
B-5 illustrates the leak detection method. The operations described above will 
be repeated four times to verify that a laboratory-grade seal is established 
that will allow operation of the VM prior to final weld closure. 

After the four closure operations are completed, the test specimen will 
be subjected to the following operations: 

1. Install the appropriate closure tool to restrain the test specimen 
in its original position 

2. Remove the threaded restraining ring 

3- Pierce the metallic sealing member to destroy the seal 

4. Reinstall the threaded restraining ring and torque to required level 

5 . Perform helium mass spectrometer leak tests to determine destruction 
of seal 

6. Perform closure welds 

7- Perform helium mass spectrometer tests of closure wel ds-- leakage not 
to exceed 1 x 10 scc/sec helium 

8. Grind closure welds to allow removal of the threaded restraining 

rings and repeat items 6 and 7 as a weld or component repair operation 

Upon completion of this test sequence, all data gathered will be cataloged 
and reported in the applicable monthly progress report to allow evaluation by 
the cognizant personnel. 
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SPECTROMETER 

Figure B-5. System Test Setup Diagram 
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